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Abstract. In this paper, the scattering and spectral theory of H = A g + V 
is developed, where A g is the Laplacian with respect to a scattering metric 
g on a compact manifold X with boundary and V £ C°°(X) is real; this 
extends our earlier results in the two-dimensional case. Included in this class of 
operators are perturbations of the Laplacian on Euclidean space by potentials 
homogeneous of degree zero near infinity. Much of the particular structure of 
geometric scattering theory can be traced to the occurrence of radial points 
for the underlying classical system. In this case the radial points correspond 
precisely to critical points of the restriction, Vq, of V to dX and under the 
additional assumption that Vq is Morse a functional parameterization of the 
generalized eigenfunctions is obtained. 

The main subtlety of the higher dimensional case arises from additional 
complexity of the radial points. A normal form near such points obtained by 
Guillemin and Schaeffer is extended and refined, allowing a microlocal descrip- 
tion of the null space of H — a to be given for all but a finite set of 'threshold' 
values of the energy; additional complications arise at the discrete set of 'ef- 
fectively resonant' energies. It is shown that each critical point at which the 
value of Vq is less than a is the source of solutions of Hu = uu. The resulting 
description of the generalized eigenspaces is a rather precise, distributional, 
formulation of asymptotic completeness. We also derive the closely related L 2 
and time-dependent forms of asymptotic completeness, including the absence 
of L 2 channels associated with the non-minimal critical points. This phenom- 
enon, observed by Herbst and Skibsted, can be attributed to the fact that 
the eigenfunctions associated to the non-minimal critical points are 'large' at 
infinity; in particular they are too large to lie in the range of the resolvent 
R(cr ± iO) applied to compactly supported functions. 
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1. Introduction 

In this paper, which is a continuation of [5], scattering theory is developed for 
symbolic potentials of order zero. The general setting is the same as in [5], consisting 
of a compact manifold with boundary, A, equipped with a scattering metric, g, and 
a real potential, V € C°°(X). Recall that such a scattering metric on J is a smooth 
metric in the interior of X taking the form 

,„ „. dx 2 h 

(L1) 9 = — + ^ 

near the boundary, where a; is a boundary defining function and h is a smooth 
cotensor which restricts to a metric on {x = 0} = dX. This makes the interior, X°, 
of A a complete manifold which is asymptotically flat and is metrically asymptotic 
to the large end of a cone, since in terms of the singular normal coordinate r — x^ 1 , 
the leading part of the metric at the boundary takes the form dr 2 + r 2 h(y, dy). In 
the compactification of X° to A, dX corresponds to the set of asymptotic directions 
of geodesies. In particular, this setting subsumes the case of the standard metric 
on Euclidean space, or a compactly supported perturbation of it, with a potential 
which is a classical symbol of order zero, hence not decaying at infinity but rather 
with leading term which is asymptotically homogeneous of degree zero. The study 
of the scattering theory for such potentials was initiated by Herbst [12] . 

Let V e C°°(dX) be the restriction of V to dX, and denote by Cv(V) the set 
of critical values of Vq. It is shown in [5] that the operator H = A g + V (where the 
Laplacian is normalized to be positive) is essentially self-adjoint with continuous 
spectrum occupying [minVo,oo). There may be discrete spectrum of finite mul- 
tiplicity in (— minx V, max Vq] with possible accumulation points only at Cv(V). 
To obtain finer results, it is natural to assume, as we do throughout this paper 
unless otherwise noted, that Vb is a Morse function, i.e. has only nondegenerate 
critical points; in particular Cv(V) is a then finite set; by definition this is the set 
of threshold energies, or thresholds. 

From the microlocal point of view scattering theory is largely about the study 
of radial points, i.e. the points in the cotangent bundle where the Hamilton vector 
field is a multiple of the radial vector field (i.e. the vector field A = J^i z id Zi on 
Euclidean space, where [z\, . . . ,z n ) € R"). These correspond in the classical dy- 
namical system to the places where the particle is moving either in purely incoming 
or outgoing sense. In scattering theory for potentials decaying at infinity, there is 
a radial point for each point on the sphere at infinity; thus there is a manifold of 
radial points and the behaviour of the flow in a neighbourhood of these points is 
rather simple, either attracting (at the outgoing radial surface) or repelling (at the 
incoming radial surface) in the transverse direction. Estimates involving commuta- 
tion with the radial vector field A multiplied by suitable powers of \z\ and perhaps 



MICROLOCAL PROPAGATION NEAR RADIAL POINTS 



3 



additional microlocalizing operators, are usually sufficient to control the behaviour 
of generalized eigenfunctions. These are known as Mourre-type estimates and play 
a fundamental role in conventional scattering theory. In the present case, assuming 
Vo is a Morse function, the radial points are isolated and occur in pairs, one pair 
(incoming/outgoing) for each critical point of Vq. The linearized Hamiltonian flow 
at the radial points is rather more complicated since it depends on the Hessian of 
Vo at the critical point, which is arbitrary apart from being nondegenerate. This 
makes the higher dimensional case more intricate than the case dim A — 2 which 
we treated in [5J . Correspondingly one needs more elaborate commutator estimates 
in order to control the behaviour of generalized eigenfunctions. We give a rather 
general and complete analysis of the regularity of solutions of Pu = in a mi- 
crolocal neighbourhood of a radial point of P, using the concept of a test module 
of operators. This is a family of pseudodifferential operators which is a module 
over the zero-order operators, contains P, and is closed under commutation. By 
choosing a test module closely tailored to the Hamilton flow of P near the radial 
point we are able to produce enough positive-commutator estimates to parametrize 
the microlocal solutions of Pu — 0. The construction of appropriate test modules 
(which can be thought of as simply an effective bookkeeping device for keeping 
track of a rather intricate set of commutator estimates) to analyze general radial 
points is the main technical innovation of this paper. 

The general study of radial points was initiated by Guillcmin and Schaeffer [4]. 
This was done in a slightly different context, where P is a standard pseudodifferen- 
tial operator with homogeneous principal symbol and a radial point is one where the 
Hamilton vector field is a multiple of the vector field J2i generating dilations in 
the cotangent space. This setting is completely equivalent to ours, via conjugation 
by a 'local Fourier Transform' (see Section |3T|) . They analyzed the situation in the 
nonresonant case. We refine their analysis by treating the resonant case, which is 
crucial in our application since we have a family of operators parametrized by the 
energy level, and the closure of the set of energies which give rise to resonant radial 
points may have nonempty interior. Moreover, we show that our parametrization 
of microlocal solutions is smooth except at a set of 'effectively resonant' energies 
which is always discrete. 

Recently, Bony, Fujiie, Ramond and Zerzeri have studied the microlocal kernel 
of pseudodifferential operators at a hyperbolic fixed point [2j, corresponding, in 
our setting, to a radial point associated to a local maximum of Vq. Their results 
partially overlap ours, being most closely related to Section 10 of [5] and [BJ. 

1.1. Previous results. The Euclidean setting described above was first studied 
by Herbst [12j . who showed that any finite energy solution of the time dependent 
Schrddinger equation, sou = e~ ltH f with / e L 2 (E"), can concentrate, in an L 2 
sense, asymptotically as t — > oo only in directions which are critical points of Vq. 
This was subsequently refined by Herbst and Skibsted [10] . who showed that such 
concentration can only occur near local minima of Vo- In contrast, solutions of the 
classical flow can concentrate near any critical point of Vq. 

Asymptotic completeness has been studied by Agmon, Cruz and Herbst pQ, by 
Herbst and Skibsted [9], [10], [11] and the present authors in [5]. Agmon, Cruz 
and Herbst showed asymptotic completeness for sufficiently high energies, while 
Herbst and Skibsted extended this to all energies except for an explicitly given 
union of bounded intervals; in the two dimensional case, they showed asymptotic 
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completeness for all energies. These results were obtained by time-dependent meth- 
ods. On the other hand the principal result of [5] involves a precise description of 
the generalized eigenspaces of H 

(1.2) E-°°{a) = {ue C~°°{X)- (H - a)u = 0}; 

note that the space of 'extendible distributions' C~°°(X) is the analogue of tempered 
distributions and reduces to it in case X is the radial compactification of 1™. Thus 
we are studying all tempered eigenfunctions of H. Let us recall these results in more 
detail. 

For any a ^ Cv(V^) the space £ , pp (cr) of L 2 eigenfunctions is finite dimen- 
sional, and reduces to zero except for a in a discrete (possibly empty) subset of 
[minx V, maxVb] \ Cv(V). It is always the case that E pp (a) C C°°(X) consists of 
rapidly decreasing functions. Hence E~^°(<t) C E~°°(<j), the orthocomplement of 
E P p(a), is well defined for a ^ Cv(V). Furthermore, as shown in the Euclidean case 
by Herbst in [TJ], the resolvent, R(a) of H, acting on this orthocomplement, has a 
limit, R{a ± «0), on [minVb,oo) \ Cv(V) from above and below. The subspace of 
'smooth' eigenfunctions is then defined as 

££>» = Sp(<r) (C°°(X) G E pp (a)) C £T°» 

(1.3) V 

Sp(cr) = (R(a + iO) - R(a - i0)) . 

2ni 

In fact 

CWC f]x-^L 2 (X). 

oo 

An alternative characterization of E^ s (a) can be given in terms of the scattering 
wavefront set at the boundary of X. 

The scattering cotangent bundle, SC T*X, of X is naturally isomorphic to the 
cotangent bundle over the interior of X, and indeed globally isomorphic to T*X by 
a non-natural isomorphism; the natural identification exhibits both 'compression' 
and 'rescaling' at the boundary. If (x, y) are local coordinates near a boundary 
point of X, with x a boundary defining function, then linear coordinates (v, fj,) are 
defined on the scattering cotangent bundle by requiring that q E SC T*X be written 
as 

(1.4) q = -v^ + Y'm—, veR, /let"- 1 . 

x * — ' X 

i 

This makes (z/, /x) dual to the basis {—x 2 d x ,xd yi ) of vector fields which form an 
approximately unit length basis, uniformly up to the boundary, for any scattering 
metric. In Euclidean space, v is dual to d r and \ii is dual to the constant-length 
angular derivative r~ 1 d Vi . In the analysis of the microlocal aspects of H — a, in 
part for compatibility with [1], it is convenient pass to an operator 'of first order' 
by multiplying H — a by x^ 1 , i.e. to replace it by 

P = P(a) = x-\H -a). 

The classical dynamical system giving the behaviour of particles, asymptotically 
near dX, moving under the influence of the potential corresponds to 'the bichar- 
acteristic vector field,' see (|2.4p . determined by the boundary symbol, p, of P. This 
vector field is defined on sc Tq X X, which is to say on SC T*X at, and tangent to, the 
boundary sc Tg X X = SC T*X n {x = 0}. It has the property that v is nondecreasing 
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under the flow; we refer to points (y, v, ji) where \x — as incoming if v < and out- 
going if v > 0. What is important in understanding the behaviour of the null space 
of P, i.e. tempered distributions, u, satisfying Pu = 0, is bicharacteristic flow inside 
{p = 0, x = 0}, a submanifold to which it is tangent. The only critical points of the 
flow are at points (y, v, 0) where y is a critical point of P and v — ±^/cr — V(y). 
Thus, the only possible asymptotic escape directions of classical particles under the 
influence of the potential V are the finite number of critical points of Vo- Moreover, 
only the local minima are stable; the others have unstable directions according to 
the number of unstable directions as a critical point of Vb : dX — ► M. 

The classical dynamics of p and the quantum dynamics of P are linked via the 
scattering wavefront set. Let u £ C~°°(X) be a tempered distribution on X (i.e. 
in the dual space of C°°(X; Q)). The part of the scattering wavefront set, WF SC («), 
of u lying over the boundary {x — 0}, which is all that is of interest here, is a 
closed subset of sc Tg X X which measures the linear oscillations (Fourier modes, in 
the case of Euclidean space) present in u asymptotically near boundary points; see 
[T5] for the precise definition. We shall also need to use the scattering wavefront 
set WF S S C («) with respect to the space x s L 2 {X) which measures the microlocal 
regions where u fails to be in x s L 2 (X). There is a propagation theorem for the 
scattering wavefront set in the style of the theorem of Hormander in the standard 
setting; if Pu £ C°°(X), then the scattering wavefront set of u is contained in 
{p = 0} and is invariant under the bicharacteristic flow of P, see [15] . In particular, 
generalized eigenfunctions of u have scattering wavefront set invariant under the 
bicharacteristic flow of P. Note that the elliptic part of this statement is already a 
uniform version of the smoothness of solutions. 

In view of this propagation theorem, it is possible to consider where generalized 
eigenfunctions 'originate', although the direction of propagation is fixed by con- 
vention. Let us say that a generalized eigenfunction originates at a radial point 
q, if q G WF sc (u) and if WF sc (it) is contained in the forward flowout $>+(q) of q; 
thus each point in WF sc (zt) can be reached from q by travelling along curves that 
are everywhere tangent to the flow and with v nondecreasing along the curve, so 
allowing the possibility of passing through radial points, where the flow vanishes, 
on the way. In Part I of this paper we showed, in the two-dimensional case and 
provided the eigenvalue a is a non-threshold value, 

• Every L 2 eigenfunction is in C°°(X). 

• Every nontrivial generalized eigenfunction pairing to zero with the L 2 
eigenspace fails to be in x~ x l 2 L 2 (X). 

• There are generalized eigenfunctions originating at each of the incoming 
radial points in {p = 0}, i.e. at each critical point of Vb with value less 
than a. 

• There are fundamental differences between the behaviour of eigenfunctions 
near a local minimum and at other critical points. The radial point corre- 
sponding to a local minimum is always an isolated point of the scattering 
wavefront set for some non-trivial eigenfunction. For other critical points, 
the scattering wavefront set necessarily propagates and in generic situa- 
tions each nontrivial generalized eigenfunction is singular at some minimal 
radial point. 

• A generalized eigenfunction, u, with an isolated point in its scattering 
wavefront set, necessarily a radial point corresponding to a local minimum 
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of Vo, has a complete asymptotic expansion there. The expansion is deter- 
mined by its leading term, which is a Schwartz function of n — 1 variables. 
The resulting map extends by continuity to an injective map from E™ s (a) 
into (B q L 2 (R n ~ 1 ), where the direct sum is over local minima of Vo with 
value less than the energy a. 
• The space E® ss (cr), consisting of those generalized eigenfunctions which 
are in x~ 1//2 L 2 microlocally near \y = 0}, is a Hilbert space and the 
map above extends to a unitary isomorphism, M + (er), from E® ss (a) to 
(B q L 2 (M n ~ 1 ). A similar map M_(er) can be defined by reversal of sign or 
complex conjugation and the the scattering matrix for P = P(o) at energy 
a may be written 

S(a) = M + (<r)MZ 1 (a). 
In this paper we extend these results to higher dimensions. 

1.2. Results and structure of the paper. We treat this problem by microlo- 
cal methods. Thus, the 'classical' system, consisting of the bicharacteristic vector 
field, plays a dominant role. The main step involves reducing this vector field to 
an appropriate normal form in a neighbourhood of each of its zeroes, which are 
just the radial points. Nondcgencracy of the critical points of Vo implies nondegen- 
eracy of the linearization of the bicharacteristic vector field at the corresponding 
radial points. If there are no resonances, Sternberg's Linearization Theorem, fol- 
lowing an argument of Guillcmin and Schaeffer, allows the bicharacteristic vector 
field to be reduced to its linearization by a contact transformation of sc Tg X X. At 
the quantum level this means that conjugation by a (scattering) Fourier integral 
operator, associated to this contact transformation, microlocally replaces P by an 
operator with principal symbol in normal form. For this normal form we construct 
'test modules' of pseudodifferential operators and analyze the commutators with 
the transformed operator. Modulo lower order terms, the operator itself becomes 
a quadratic combination of elements of the test module. Just as in Part I, we use 
the resulting system of regularity constraints to determine the microlocal structure 
of the eigenfunctions and ultimately show the existence of asymptotic expansions 
for eigenfunctions with some additional regularity. 

However, the problem of resonances cannot be avoided. Even for a fixed operator 
and fixed critical point, the closure of the set of values of a for which resonances 
occur may have non-empty interior. Such resonances prevent the reduction of the 
bicharacteristic vector field to its linearization, and hence of the symbol of P to 
an associated model, although partial reductions arc still possible. In general it 
is necessary to allow many more terms in the model. Fortunately most of these 
terms are not relevant to the construction of the test modules and to the derivation 
of the asymptotic expansions. We distinguish between 'effectively nonrcsonant' 
energies, where the additional resonant terms are such that the definition of the 
test modules, now only to finite order, proceeds much as before and the 'effectively 
resonant' energies, where this is not the case. Ultimately, we analyze the regularity 
of solutions at all (non-threshold) energies. Near effectively nonresonant energies, 
smoothness of families of eigenfunctions may still be readily shown. Effectively 
resonant energies are harder to analyze, but the set of these is shown to be discrete. 
In any case, the space of microlocal eigenfunctions is parameterized at all non- 
threshold energies. At effectively resonant energies the problems arising from the 
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failure of the direct analogue of Sternberg's linearization are overcome by showing 
that, to an appropriate finite order, the operator may be reduced to a non-quadratic 
function of the test module. 

In outline, the discussion proceeds as follows. In sections [2] - [4] we study radial 
points. This is a general microlocal study except that we work under the assumption 
that the symplectic map associated to the linearization of the flow at each radial 
point (see Lemma I2.5j) has no 4-dimensional irreducible invariant subspaces; this 
assumption is always fulfilled in the case of our operator A + V — a. The main result 
is Theorem 13 . 1 1 1 in which the operator is microlocally conjugated to a linear vector 
field plus certain 'error terms'. In the nonresonant case the error terms can be 
made to vanish identically, while in the effectively nonresonant case the error terms 
have a good property with respect to a test module of pseudodiffcrential operators, 
namely they can be expressed as a positive power x e , e > 0, times a power of the 
module. In the effectively resonant case this is no longer possible and we must allow 
'genuinely' resonant terms, but the set of effectively resonant energies is discrete in 
the parameter a in all dimensions. 

We then turn in sections [5j - [7] to studying microlocal eigenfunctions which are 
microlocally outgoing at a given radial point q. The main result here is Theorem 16 .71 
(or Theorem [73] in the effectively resonant case) which gives a parameterization of 
such microlocal eigenfunctions. For a minimal radial point, they are parameterized 
by 6>(R ra_1 ), Schwartz functions of n— 1 variables, for a maximal radial point they 
are parameterized by formal power series in n— 1 variables, and in the intermediate 
case of a saddle point with k positive directions, they are parameterized by formal 
power series in n — 1 — k variables with values in S(R k ). In all cases, the parame- 
terizing data appear explicitly in the asymptotic expansion of the eigenfunction at 
the critical point. 

We next investigate in sections [8] and [9] the manner in which the various radial 
points interact, and prove, in Theorem 19. 2\ a 'microlocal Morse decomposition.' 
This shows that for each non-threshold energy a there are genuine eigenfunctions 
(as opposed to microlocal eigenfunctions) in E^ s (a) associated to each energy- 
permissible critical point. 

Then we turn in sections [TOl and [Til to the spectral decomposition of P and prove 
several versions of asymptotic completeness. First this is established at a fixed, 
non-threshold energy; see Theorem 110.11 which shows that the natural map from 
E® ss (a) to the leading term in its asymptotic expansion (i.e. to its parameterizing 
data) is unitary. Next we prove a form valid uniformly over an interval of the 
spectrum, Theorem 110.101 In section [TTJ a time-dependent formulation is derived, 
as Theorem 111.41 This is based on the behaviour at large times of solutions of the 
time-dependent Schrodinger equation D t u = Pu and is subsequently used to derive 
a result of Herbst and Skibsted's on the absence of L 2 -channels corresponding to 
non-minimal critical points ( Corollary 1 11.7[) . 

1.3. Results used from [5j. Throughout this paper we state the specific location 
of results used from [5] . For the convenience of the reader we summarize here the 
relevant locations. Sections 1-3 of [5] are used as the basic background (and [5J 
Section 3] relies on Section 4 there). The present Section [4] is the analogue of [5j 
Section 5] , although we restate many of the arguments due to the slightly different 
(more general) setting. The basic analytic technique using test modules in Section[5j 
comes from [5j Section 6]. Certain results and methods from Sections 11 and 12 
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of [5] are used here in Sections [9] and [TUl However, the results of the intermediate 
sections 7-10 of [5], while certainly of interest when comparing to the results of 
Section [5] and [7] here, are never used in the present work directly or indirectly. 

In addition, there was an error in the proof of Proposition 6.7 of [5]. While this 
error is minor and is easily remedied, we present the modified proof, together with 
some of the context, here in the Appendix since this proposition lies at the heart 
of the analysis in both papers. 

1.4. Notation. The items listed below without a reference whose definition is not 
immediate from the stated brief description are defined in [15) . 



Notation Description/definition of notation Reference 

Vb restriction of V to dX 

Cv(V) set of critical values of Vq 

SC T*X scattering cotangent bundle over X (|1.4ll 

sc T* ax X restriction of SC T*X to dX (|L4l 

x boundary defining function of X s.t. holds 

y coordinates on dX 

(is, fibre coordinates on SC T*X (jl.4p 

V = (y' \ y", y'") decomposition of y variable (12.12| 
^ = (^', fi", //") dual decomposition of [i variable (12.12| 
r 'ii r j ' r k eigenvalues of the contact map A (|2.12p 
Y" y^/x r " §M 
Y«> y'l'/x^ (E23D 
A (positive) Laplacian with respect to g 

P x^iA + V - a) Sec. El 

H A + V 

R(a) resolvent of H, (H - a)' 1 

R(a ± iO) limit of resolvent on real axis from above/below 

V modified potential Lem. 18.51 
Sp(cr) (generalized) spectral projection of H at energy a (jl.3p 
R(a) resolvent of modified potential (A + V — a)^ 1 

Lg C (X) L 2 space with respect to Riemannian density of g 

H™>°(X) Sobolev space; image of L 2 SC {X) under (1 + A)" m / 2 

H™' l (X) x l H™>°(X) 

ty™'°(X) scattering pseudodiff. ops. of differential order m 

*™' Z (X) x^^iX); maps H™'' l '(X) to H™'- ,n ' l '+ l (X) 

aaA A ) boundary symbol of A G C°° fn. on SC T£ X X 

a d (A) cr 9)0 (A) 

WF sc (m) scattering wavefront set of u; closed subset of sc Tg X X 

WF™''(u) scattering wavefront set with respect to H™> 1 

WF' SC (A) operator scattering wave front set; in its complement 

A is microlocally in ^>*f°(X), i.e. is trivial 

sc H p scattering Hamilton vector field Sec. [5] 

$ + (g) forward flowout from q G sc Tg X X Sec. 1.1 

radial point point in sc Tq X X where p and sc H p vanish Sec. [2] 

RP±(ct) set of radial points of H — a where ±v > 

Min + (er) subset of RP + (er) associated to local minima of Vq 

< partial order on RP + (a) compatible with $ + Def. 18.31 
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£-mic,+ (0, P) microlocal solutions of Pu = in the set O (|4.1I) 

Emic,+ (q, cr) microlocal solutions of (H — a)u — near q (14.411 

E^ ss (cr) space of generalized cr-eigenfunctions of H (|9.ip 

E S (T, a) subset of u G E s css (<r) with WF sc (u) n RP+ (cr) C T $91$ 

Ekin,+ (*) E s (r, Cr), With T = Min + (a) 

test module Sec. [5] 

Isc (0,A4) space of iteratively-regular functions w.r.t. .M (|5.9p 

t rescaled time variable; r = xt Sec. [TT] 

X Sch X x K r dm 



Acknowledgement. We thank an anonymous referee for helpful comments and 
for pointing out some errors in the first version of this manuscript. 

2. RADIAL POINTS 

Let X be a compact n-dimensional manifold with smooth boundary. Recall that 
if (x, y) are local coordinates on X, with x a boundary defining function, then dual 
scattering coordinates (v, /i) on the scattering cotangent bundle are determined. 
The restriction of the scattering cotangent bundle to dX is denoted sc Tg X X and 
has a natural contact structure, the contact form at the boundary being 

(2.1) a = -dv + 22 Hjdyj 

i 

in local coordinates. Recall that a contact structure on a 2n — 1-dimensional man- 
ifold, here sc Tg X X, is given by a nondegenerate one- form, i.e. a one- form a with 
a A (da)™ -1 everywhere non-zero; correspondingly its kernel is a maximally non- 
integrable hyperplanc field on sc Tg X X. One refers to either the line bundle given by 
the span of a, or the hyperplane field given by its kernel, as the contact structure. 

Suppose that P € V E'* ( ; _1 (X) is a scattering pseudodifferential operator of order 
— 1 at the boundary; for example, P = x _1 (A + V — o~). Then the boundary part 
of its principal symbol, p = ag(xP), is a C°° function on sc Tq X X. In this, and 
the next, section we consider radial points of a general real- valued function, p € 
C°°( sc Tg X X), with only occasional references to the particular case, p = \(\ 2 +Vo—<t, 
of direct interest in this paper. Although we discuss radial points in the context 
of boundary points in the scattering calculus this analysis applies directly (and 
could alternatively be done for) radial points in the usual microlocal picture, as 
described in the Introduction. Our objective in this section is to find a change of 
coordinates, preserving the contact structure, in which the form of p is simplified. 
In this section we consider the simplification of p up to second order, in a sense 
made precise below. 

The basic non-degeneracy assumption we make is that 

(2.2) p = implies dp ^ 0; 

this excludes true 'thresholds' which however do occur for our problem, when a is 
a critical value of Vo- It follows directly from (|2.2ll that the boundary part of the 
characteristic variety 

E = {q G sc TQ X X;p{q) = 0} is smooth; 
we shall assume that E is compact, corresponding to the ellipticity of P. 
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Definition 2.1. A radial point for a function p satisfying ()2.2j) is a point q € E 
such that is a (necessarily nonzero) multiple of the contact form a given by 

(|2.1[) . Conversely, if q € E and dp and a are linearly independent at g then we say 
that p is of principal type at g. 

We may extend p to a C°° function on SC T*A, still denoted by p. Over the interior 
sc Tj„J is naturally identified with T*X°, which is a symplectic manifold with 
canonical symplectic form to. Near the boundary, expressed in terms of scattering- 
dual coordinates, 

. I dx \ -> dyi \ \ -» dx sr-^ dyi 

(2.3) u = d\-v-^ + 2^iH— I - {-dv + ^ii l dy l )A — + ^d^^— . 

\ i / i i 

Consider the Hamilton vector field, H x -i p , of x~ x p, which we shall denote sc H p , 
fixed by the identity ui(-, sc H p ) — dp. Then sc H p extends to a vector field on 
SC T*X tangent to its boundary, so sc H p € V b ( sc T*A)Q At the boundary sc H p , 
as an element of Vb( sc T*A), is independent of the extension of p. We denote the 
restriction of sc H p (as a vector field) to sc Tg X X by W, so W is a vector field on 
sc Tg X X. Explicitly in local coordinates 

SC H P = - (d u p)(xd x + n- ■ dfj,) + {xd x p -p + fi- d^p)d v 

(2 ' 4) + E iPihP^i ~ d v,P d ^) + ^ b ( sc T*A); 

i 

since p is smooth up to the boundary, xd x p — at sc Tg X X. Thus, 

(2.5) W = - {d v p)i± -d^ + ifx- d p p - p)d v + ^ I''>; ~ d y 3 P d H ) ■ 

3 

Alternatively W may be described in terms of the contact structure on sc Tq X X. 
Namely W is the Legendre vector field of p, determined by 

(2.6) da(.,W) +70; = dp, a(W) =p 

for some function 7. It follows that W is tangent to S, since dp(W) = r ya(W) = 
7P = at any point at which p vanishes. An equivalent definition of q £ S being a 
radial point is that the vector field W vanishes as q, as follows from (|2.6[) and the 
nondegeneracy of a. 

Definition 2.2. A radial point q G E for a real- valued function p e C°°( sc Tg X X) 
satisfying (|2.2|) is said to be non-degenerate if the vector field W, restricted to 
E = {p = 0}, has a non-degenerate zero at q. Note that this implies that a non- 
degenerate radial point is necessarily isolated in the set of radial points. 

Since the vector field W vanishes at a radial point q, its linearization is well 
defined as a linear map, A 1 on T q sc Tg X X, (later we will use the transpose, A, as a 
map on differentials) 

(2.7) A'v=[V,W](q), 

for any smooth vector field V with V(q) — v; it is independent of the choice of 
extension and can also be written in terms of the Lie derivative 

(2.8) A'v = -C w V{q). 



^Here Vb(M) denotes the space of smooth vector fields on the manifold with boundary M that 
are tangent to dM . 
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Since Wp = jp, A' preserves the subspace T q Y>. Since a is normal to T g S, the 
restriction of da to T q Y> is a symplectic 2-form, u> q . 

Lemma 2.3. At a non-degenerate radial point for p, where dp = Xa, the lineariza- 
tion A' acting on T g E is such that 

S = A'- -AIdesp(2(n-l)) 
is in the Lie algebra of the symplectic group with respect to uj q : 

uj q (Svi,v 2 ) + u) q (vi, Sv 2 ) — 0, V Vi,v 2 S 
Proof. Observe that (|2.6II implies that 
(2.9) L w a = (da)(W,-) + d(a(W)) = ~/a. 

For two vector smooth vector fields Vi, defined near q, 

W(da(V 1 ,V 2 )) = L w (da(V 1 ,V 2 )) 
' 2 ' " ' ! = (L w da)(Vi, V 2 ) + da(L w Vi,V 2 ) + da(Vi,L w V 2 ). 

The left side vanishes at q so using (|2.8p 

(2.11) oj q (A'vi,v 2 ) +u q (v 1 ,A / v 2 ) = XoJ q (vi,v 2 ) V v 1 ,v 2 e T 9 E. 

□ 

It follows from Lemma 12.31 see for example [4] , that A' is decomposable into 
invariant subspaces of dimension 2 and 4, with eigenvalues on the two-dimensional 
subspaces of the form Xr, A(l — r), r < 1/2 real or A(l/2 + ia), A(l/2 — ia), with 
a > 0. 

Note that, by (12. 6|) . d v p(q) = —7(g) = —A, so from (|2.4p . the Hamilton vector 
field sc i/ p is equal to Xxd x modulo vector fields of the form / • W' where W 
is tangent to {x = 0} and f(q) — 0. Therefore if A > 0, then x is increasing 
along bicharacteristics of p in the interior of SC T*X, i.e. the bicharacteristics leave 
the boundary, i.e. 'come in from infinity' if dX is removed, while if A < 0, the 
bicharacteristics approach the boundary, i.e. 'go out to infinity'. Correspondingly 
we make the following definition. 

Definition 2.4. We say that a non-degenerate radial point q for p with dp(q) = 
Xa(q) is outgoing if A < 0, and we say that it is incoming if A > 0. 

For p = \(\ 2 + Vo — cr, we have A = — d„p — —2v. Hence, radial points are 
outgoing for v > and incoming for v < in this case. We next discuss the form 
the linearization takes for p — |£| 2 + Vq — a. 

Lemma 2.5. For the function p = |£| 2 + Vq — a with Vo Morse, the radial points 
are all nondegenerate and the linear operator S associated with each has only two- 
dimensional invariant symplectic subspaces. 

Remark 2.6. In view of the non-occurrence of non-decomposable invariant sub- 
spaces of dimension 4 in this case we will exclude them from further discussion 
below. 

Proof. Choose Riemannian normal coordinates yj on dX, so the metric function 
h satisfies h — \p\ 2 = 0(\y\ 2 ). Since the Hessian of V\qx at a critical point is a 
symmetric matrix, it can be diagonalized by a linear change of coordinates on dX, 
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given by a matrix in SO(n — 1), which thus preserves the form of the metric. It 
follows that for each j, {dyj, dfJ,j) is an invariant subspace of A. □ 

Let I denote the ideal of C°° functions on sc Tg X X vanishing at a given radial 
point, q. The linearization of W then acts on T* ( sc Tg X X) = I/I 2 ; dp(q), or equiv- 
alently a q , is necessarily an eigenvector of A with eigenvalue 0. Similarly, sc H p 
defines a linear map A on T* ( SC T*X) . By (|2.4[) . A preserves the conormal line, 
span da; and the eigenvalue of A corresponding to the eigenvector dx is A. Thus A 
acts on the quotient 

T* ( sc n x X) = T* q ( sc T*X) I spanete, 

and this action clearly reduces to A. 

By Darboux's theorem we may make a local contact diffeomorphism of sc Tg X X 
and arrange that q = (0,0,0). Thus, as a module over C°°( sc Tg X X) in terms of 
multiplication of functions, I is generated by v, yj and the fa, for j — 1, . . . , n — 1. 
Thus in general we have the following possibilities for the two-dimensional invariant 
subspaces of A. 

(i) There are two independent real eigenvectors with eigenvalues in A(R\[0, 1]). 

(ii) There are two independent real eigenvectors with eigenvalues in A(0, 1). 

(iii) There are no real eigenvectors and two complex eigenvectors with eigen- 
values in A(| +i(R\ {0})). 

(iv) There is only one non-zero real eigenvector with eigenvalue ^A. 

Case ([Iv]) was called the 'Hessian threshold' case in Part I. In all cases the sum of 
the two (generalized) eigenvalues is A. 

Lemma 2.7. By making a change of contact coordinates, i.e. a change of co- 
ordinates on sc Tg X X preserving the contact structure, near a radial point q for 
p G C°° : ( sc Tg X X) for which the linearization has neither a Hessian threshold sub- 
space, (jrv|) . nor any non-decomposable ^-dimensional invariant subspace, coordi- 
nates y and /i, decomposed as y = (y',y",y"') and \i — (//, //', p!"), may be intro- 
duced so that 

(i) 

(2-12) (j/',mO = (2/1, ■ ■ •>lfo-i,Mi) ■ ■ ->Ms-i) 

where e'j — dy'j, /j = dfi'j are eigenvectors of A with eigenvalues Xr'j, 
A(l — r'j), j = 1, . . . , s — 1 with r'j < real and negative. 

(ii) = (ys,---,Vm-i,Vs,---,fim-i) where e'j = dy'j , f" = dfj," are 
eigenvectors with eigenvalues Ar", A(l — r'j), j = s, . . . , m — 1 where < 
r'j < 1/2 is real and positive. 

(iii) (y'", jj!") = (y m , ■ ■ ■ , y n -i, Mmi ■ ■ ■ , Mn-i)> where some complex combina- 
tion e'j' , /j", of dy'j' and dfj/j', m < j < n — 1, are eigenvectors with 
eigenvalues Xr'j' and A(l - r'j') with r'j' = 1/2 + ifi'j' , f3'j' > 0. 

Thus if we set e = (e', e" , e'"), f = (/', /", /"') the eigenvectors of A are dv, ej and 
fj, with respective eigenvalues 0, Xrj and A(l — Tj); we will take the coordinates so 
that the Tj are ordered by their real parts. 

Remark 2.8. We emphasize that the change of coordinates here is on the contact 
space, sc Tq X X, and it is, in general, not induced by a change of coordinates on X. 
Analytically it is implemented by a scattering FIO (see Section l3T|) . 
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In coordinates in which the eigenspaces take this form it can be seen directly 
that 

m— 1 n— 1 

(2.13) p = A( - v + ^ rjVjVj + QoiVj^i) + vg\ + g 2 ) 

j=i j=m 

with the Qj elliptic homogeneous polynomials of degree 2, g\ vanishing at least 
linearly and g 2 to third order. 

Remark 2.9. For the function p = |£| 2 + Vq — a with Vq Morse, the eigenvalues 
of A at a radial point q are easily calculated in the coordinates used in the proof 
of Lemma 12.51 Indeed, since the 2-dimensional invariant subspaces decouple, the 
results of [5J Proof of Proposition 1.2] can be used. The eigenvalues corresponding 
to the 2-dimensional subspace in which the eigenvalue of the Hessian is 2a j are thus 




V Q (0) 



where A = —2u(q). 



In fact, below we do not need the full power of Lemma 1^771 Essentially it suffices 
if we arrange that the eigenvectors corresponding to the (in absolute value) larger 
eigenvalues, namely A(l — r'j), if r'- < 0, or A(l — r'j), if r'j G (0, |), are in a model 
form on the two dimensional eigenspaces. The advantage of the weaker conclusion 
is that one has more freedom in choosing the contact change of coordinates. 

Lemma 2.10. (Weaker version of Lemma \2.7\ ) Suppose that ^A is not an eigen- 
value of A. By making a change of contact coordinates, i.e. a change of co- 
ordinates on sc Tg X X preserving the contact structure, near a radial point q for 
p G C x : ( sc Tq X X) for which the linearization has neither a Hessian threshold sub- 
space, (try)) . nor any non-decomposable ^-dimensional invariant subspace, coordi- 
nates y and /i, decomposed as y = [y' ,y" ,y"') and [i = (//' , ji" , //"), may be intro- 
duced so that 

(i) 

(2-14) (y, //) = (y x , . . . , Mi, ■ ■ ■ , Ms-i), 



where some real linear combinations e'j of d/j,'j and dy'j , resp. /j = d/j/j are 
eigenvectors of A with eigenvalues Xr'j, resp. A(l — r'-), j — 1, . . . , s — 1 
with r'j < real and negative. 
(ii) (y", fi") = (y s , . . . , y m -i, p^s, ■ ■ ■ , ^m-i) where some real linear combina- 
tions e'j of dn" and dy'J , resp. f" = dfi" are eigenvectors with eigenvalues 
Xr'- , A(l — r"), j = s, . . . , m — 1 where < r" < 1/2 is real and positive. 
(hi) (?/"', pL 1 ") = (y m , ■ ■ ■ ! Vn-i, t^mi ■ ■ • , f-n-i), where some complex combina- 
tion e'j' , f'j" , of dy'j and dfj,"' , m < j < n — 1, are eigenvectors with 
eigenvalues Xr'J' and A(l - r'j') with r'j' = 1/2 + i(5'j' , f3'j' > 0. 

Again, if we set e = (e', e", e'"), f = (/', /", /"') the eigenvectors of A are dv, ej and 
fj, with respective eigenvalues 0, Xrj and A(l — rj); we will take the coordinates so 
that the rj are ordered by their real parts. In these coordinates a version of (|2. 13|) 
still holds, namely if dj and bj are any functions on sc Tg X X vanishing at (0, 0, 0) 
with differential ej, resp. fj,j = l,...,m—l (so we may take bj = \ij, and we 
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may take aj a R-linear combination of yj and Lij ) then 
m— 1 n— 1 

P = A ( - v + r i a i h i + QjiVhVi) + v 9i + 92) 

7 — 1 7=771 

2.15) 

= A ( ~ v + + C J^I + Y QAVi^i) + v 9i + .92), 
i=i i=l j=m 

where the Cj are real, the Qj are elliptic homogeneous polynomials of degree 2, 171 
vanishes at least linearly and gi to third order. 

As mentioned, Lemma 12.101 is weaker than, hence is an immediate consequence 
of, Lemma 12.71 Although it is by no means essential, this weaker result leaves 
more freedom in choosing the contact map which is useful in making the choice 
rather explicit, if this is desired. In fact, if p = |£| 2 + Vq — a, as in Lemma \2 .51 we 
immediately deduce the following. 

Lemma 2.11. For the function p = |£| +Vq — cr with Vq Morse, the contact map in 
Lemma \2.10\ can be taken as the composition of the contact map on sc Tg X X induced 
by a change of coordinates on X , with the canonical relation of multiplication by a 
function of the form e^l x , cj> g C°°(X). 

Remark 2.12. The canonical relation of multiplication by e % ^l x is given, in local 
coordinates (y, v, fx), by the map 

$0 : (y, v, fx) i-> (y,v + <p(y), ll + d v <p(y)), 

i.e. if we write ^^{y^v, fx) = (y,i>,p,), then p,k = Lit + d Vk <f>(y). Note that while 4> 
is a function on X, the canonical relation only depends on (j)\dx, which is why we 
simply regard <fi as a function on dX and write </>(y) here. 

Proof. As in the proof of Lemma 12.51 we may assume, by a change of coordinates 
on X, that the critical point of Vq over which the radial point q lies is y = 0, 
that h — \fi\ 2 = 0(\y\ 2 ) and that the Hessian of Vb at is diagonal, so for each 
j, (dyj,dfXj) is an invariant subspace of A. Note that in the coordinates (y,v,fx), 
q = (0,1^0,0). With the notation of Remark 12.91 above, if dyj is an eigenvector of 
the Hessian with eigenvalue 2a j then the eigenvectors of A of eigenvalue Xrj , resp. 
A(l — Tj), are ej — ^(1 — r.j)dyj + d/ij, resp. fj = ^rjdyj + dfXj, see Remark 1.3 of 
[5]. In particular, if rj is real, so is fj. 

Now, the contact map <&0 induced by multiplication by e l ^/ x as above acts on 
T* SC T* X X by pull-backs, namely 

y*k dyk + D * dD + J2f 1 *k 

k k 

= ^Vk dyk + D * ( dv + Y ( dy i & dy i ) + Y P-l (d^k + Y d vi dyMv) d Vj ) ■ 

k j k j 

Thus, by the above remark, $ will map q to (0, 0, 0) provided (f)(0) = —vq, d yj <j>(0) = 
for all j. In this case, moreover, the pull-back $^ will map dyk to dyk, dv to dv and 

dLi k to dLi k + J2j 9 yj d yk (j)(y). Correspondingly, by letting <p(y) = -vo + Y^Ji b jVj> 
bj = jTj, ( < i ) 1 )* maps fj to dtij, j — l,...,m — 1. Since the Legendre vector 
field W of ($^)* p is the push-forward of the Legendre vector field W of p under 
<&0, it follows that d/Xj is an eigenvector of the linearization of W with eigenvalue 
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A(l — rj). As $J also maps the 2-dimensional subspaces (dyj,d[ij) (at (0,0,0)) to 
the 2-dimensional subspaces (dyj,dfij) (at q), and the latter are invariant under A, 
so are the former under the linearization of W . This proves the lemma. □ 



3. MlCROLOCAL NORMAL FORM 

Let P £ ^l^~ 1 (X) be an operator with real principal symbol p obeying (|2.2p . as 
in the previous section, and assume that q is a nondegenerate radial point for p. In 
this section we shall reduce p to a normal form, via conjugation with a scattering 
Fourier integral operator. We first pause to define such operators. 

3.1. Scattering Fourier integral operators. Scattering Fourier integral opera- 
tors (FIOs) are defined in terms of conventional FIOs via the local Fourier trans- 
form, as defined in [17j . Let X be a manifold of dimension n with boundary, and 
(x, y) local coordinates where x is a boundary defining function. We can always 
identify a neighbourhood U C dX of yo £ dX with an open set V £ S*™" 1 , which 
we can think of as embedded in R" in the standard way. Correspondingly we may 
identify the interior of a neighbourhood [0, e)^ x U C X of (0, yo) £ X with the an 
asymptotically conic open set (e , oo) x V C R" in R™. If we choose a function 
4> £ C°°(X) supported in [0, e) x x U which is identically 1 in a neighbourhood of 
(0, yo), then the operator T with kernel 

is called a 'local Fourier transform' on X. Here z — [z\, . . . , z n ) £ R™, z ■ y denotes 
the inner product on R™ and du>(y) denotes the standard measure on S'™ -1 (pulled 
back to dX and then to X via the identifications above). Of course, if X is the 
radial compactification of R™ and the identification between U and V is the identity, 
then T really is the Fourier transform premultiplied by the cutoff function <j>. 

It is shown in [17] that T induces a local bijection between sc Tg X X and the 
cosphere bundle of R". In fact, using our identification between U and V C S n ~ 1 
we may represent points in SC T^X as (z, £) where z, = z/\z\ £ V represents a point 
in U and £ represents the point in the fibre given by (y, \x) where v is the parallel 
component of ( relative to 2 and /i is the orthogonal component. The identification 
is then given by the Legendre map 

L(2,Q = (C-z) £ S*R n . 

In other words, T sets up a bijection between scattering wavefront set and 
conventional wavefront set. Moreover, it is shown in [17j that conjugation by T 
maps the scattering pseudodifferential operators A £ ^*£{X) microsupported near 
(yo, Vq, ho) to the conventional pseudodifferential operators microsupported near 
L(yo, fo, jUo) ; with principal symbols related by 

a l {TAT*){L{q))=a{q) 1 

where a is the boundary symbol of A (of order I). 

Definition 3.1. A scattering FIO is an operator E from C co (X) to C~°°(X) such 
that, for any local Fourier transforms T\, T2 on X, TiET\ is a conventional FIO 
onK". 
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A simple example of a scattering FIO is multiplication by an oscillatory factor 
e #(y)/^ Under conjugation by a local Fourier transform this becomes a conven- 
tional FIO given by an oscillatory integral with phase function (z— ^O'C+IClVKC/lCD- 
The scattering resolvent kernel constructed in [7] and [5J, microlocalized to the in- 
terior of the 'propagating Legendrian', is another example. 

It follows then that we can find a scattering FIO quantizing any given contact 
transformation from a neighbourhood of a point q £ sc Tq X X to itself, since we 
may conjugate by a local Fourier transform and reduce the problem to finding a 
conventional FIO quantizing a homogeneous canonical transformation from a conic 
neighbourhood of L{q) £ S*M. n to itself. We can also use the local Fourier transform 
to import Egorov's theorem into the scattering calculus. Namely, if B £ ^ti 1 (^0 is 
a scattering pseudodifferential operator of order —1, with real principal symbol, and 
P £ then also e~ iB Pe iB £ ^l^ 1 (X) is a scattering pseudodifferential 

operator of order —1, whose symbol p' is related to that of P by the time 1 flow of 
the Hamilton vector field of B. This indeed is how we shall conjugate the principal 
symbol p of our operator to normal form. 

3.2. Normal form. In this section we put the principal symbol of P into a normal 
form Pnorm- For later purposes we shall also need the subprincipal symbol of P 
in a normal form, but only only along the 'flow-out', i.e. the unstable manifold, 
of q, which can be done via conjugation by a function; this is accomplished in 
Lemma 16. II (The model form of the subprincipal symbol only plays a role in the 
poly homogeneous, as opposed to just conormal, analysis, which is the reason it is 
postponed to Section [6]) 

For this purpose, we only need to construct the principal symbol u(B) of B as 
in the first subsection. This in turn can be be written as x~ 1 b, b £ C°°( SC T*X), 
so we only need to construct a function b on sc Tg X X such that the pull-back Q*p 
of p by the time 1 flow $ of H x -±i is the desired model form p n orm, where b is 
some extension of b to SC T*X; this property is independent of the chosen extension. 
Thus any B with a(B) — b will conjugate P to an operator with principal symbol 
Pnorm- This construction is accomplished in two steps, following Guillcmin and 
Schaeffer [3] in the non-resonant setting. First we construct the Taylor series of b 
at q = (0, 0, 0), which puts p into a model form modulo terms vanishing to infinite 
order at q. Next, we remove this error along the unstable manifold of q by modifying 
an argument due to Nelson [15] . 

Rather than using powers of X to filter the Taylor series of b, we proceed as in [4] 
and assign degree 1 to y and [i but degree two to v in local coordinates as discussed 
above. Thus, let denote the space of functions 

= u a y a ^C ca { sc T^ x X) 

2a+H + |/3|-2=j 

Note that this is well-defined, independently of our choice of local coordinates, since 
— dv is the contact form a at q, so v is well-defined up to quadratic terms. The 
Poisson bracket preserves this filtration of T in the following sense. If a, b are some 
smooth extensions to SC T*X of elements a £ if, b £ f) J then 

x^c = {x^a^^b} =>c = c|sc T . xX £ f)' i+J . 
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When this holds we write c = {{a, &}}; explicitly, 

f)o f)b 

(3.1) {{a , 6}} = Wa (6) + ^&-— a, 
with W given by ([2~5]) . Thus 

(3.2) {{..•!! :()' •'!)•'• -11"'. 
We then consider the quotient 

so the bracket {{-,-}} descends to 

fl' £1' 

Remark 3.2. These statements remain true with replaced by X J . However, note 
that p = —v in I/I 2 , since dp = —dv at g, but it is not the case that p = —v in g°. 
In fact, p is given by (|3 .3(1 below in g°. 

Using contact coordinates as discussed above, gi may be freely identified with 
the space of homogeneous functions of v, y, [i of degree j + 2 where the degree of 
v is 2. Now let po be the part of p of homogeneity degree two. In order to use 
Lemmas 12.71 and 12.101 we assume throughout the paper from here on that case (iv) 
above Lemma \2. 7| does not apply. Hence from (I2.13P 

m— 1 7i — l 

(3.3) po = A( - v+ ^2 r jVj^3 + ^2 Qj(Vj^j))> P-Po e if 1 - 

j=l j=m 

If we take b G i) 1 , I > 1 and let $ be the time 1 flow of H x -n, then 

(3.4) x$*(x- 1 p)=p + {{p,b}}=p + {{p ,b}}, modulo h i+1 . 

This allows us to remove higher order term in the Taylor series of the symbol 
successively provided we can solve the 'homological equation' 

{{po,b}} = e£f) 1 , modulo f/ +1 . 

Thus we need to consider the range of this linear map; its eigenfunctions are easily 
found from the eigenfunctions of the linearization of W. 

Lemma 3.3. The (equivalence classes of the) monomials p%e a f® with 2a + \a\ + 
\f3\ = I + 2 satisfy 

{{POiPo^f 13 }} = Ra,a.j3Poe a f with eigenvalue 

(3.5) / ^ 'i^ 1 

\ j'=l J'=l 

and give a 6asis o/ eigenvectors for {{po, .}} acting on q 1 . 

Here we identify the differentials ej and fj with linear functions with these dif- 
ferentials. 

Remark 3.4. In fact, the contact coordinates given by Lemma [2.101 suffice for the 
proof of this lemma; the additional information in Lemma 12.71 is not needed. In 
this case, by (|2.15p . 

m — 1 n — 1 

(3.6) po = A( - v + r j e jfj + Qj(Vj,H))' 
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We also remark that we could equally well use the eigenvector basis for {{po, .}} 
acting on q 1 given by v a e a f@ with 2a + \a\ + =1 + 2. This follows from the 
lemma using that v = J2T=i r iVo^3 + YTj= m QjiVjiVj) ~ ^~ 1 Po hi £j°, and y 3 fij as 
well as Qj{yj,^j) are eigenvectors with eigenvalue X(rj + (1 — rj)) = A, and so is 
Po- 

Proof. Taking into account the eigenvalues and eigenvectors of A, all eigenvalues 
and eigenvectors of {{po, .}} can be calculated iteratively using the derivation prop- 
erty of the original Poisson bracket. This implies 

{{p , ab}} = x{x~ 1 p a , x(x~ 1 a)(x~ 1 b)} 

(3.7) = x~ 1 {x~ 1 po, x}ab + x{x~ 1 po, x~ 1 a}b + xa{x~ 1 po, 
= Xab + {{p a , a}}b + a{{p , b}}, 

where each term within {., .} really uses a C°° extensions of the a, b, po to SC T*X, 
followed by evaluation of the bracket and then restriction to sc Tg X X. Since 

{{Po, «}} = x{x~~ po, x~ 1 a} = x{x~ 1 p , x~ 1 }a + {x~ 1 p , a} = — + {%~ Po-, a}, 

on the eigenvectors of {{po, ■}} are the eigenvectors ej and fj of A with eigen- 
values — X + Xrj and — A + A(l — rj). Moreover, in g°, p a is an eigenvector of {{po, ■}} 
with eigenvalue 0. Thus, ej, fj and po satisfy the claim of the lemma. Since the 
other generators of g°, as well as generators of Q 3 , j > 1, can be written as a 
products of the ej, fj and po, the conclusion of the lemma follows by induction. □ 

Definition 3.5. We call the multiindices in the set 

(3.8) / = {(a, a, (3); i? Q ,a,/3 = and 2a + \a\ + > 3} , 
with R a , a ,f3 given by (|3.5[) . resonant. 

Conjugation therefore allows us to remove, by iteration, all terms except those 
with indices in /. Expanding using (|3.3p we deduce the following. 

Proposition 3.6. If P is as above and the leading term of p = ag ^i(P) is given 
by (|3.3p near a given radial point q then there exists a local contact diffeomorphism 
$ near q such that 

m n—1 

(o.y) 



j = l j=m+l (a,a,/3)£l 



modulo 1°° — at q 



with I given by 



Proof. The Taylor series of $ at q can be constructed inductively over the filtration 
f)- 7 as indicated above. At the jth stage, the terms of weighted homogeneity j can 
be removed from p except for those in the null space of {{po, •}}, i.e. the resonant 
terms with R a . a ,p — 0. This leads to (|3.9p in the sense of formal power series. 
However, by use of Borel's Lemma a local contact diffeomorphism can be found 
giving dill). □ 



Now a small extension of Nelson's proof of Sternberg's linearization theorem can 
be used to remove the infinite order vanishing error along the unstable manifold, 
i.e. at v = 0, /i = 0, y" = 0, y'" = 0. 
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Proposition 3.7. Suppose that X and Xq areC°° vector fields onH N with -Xo(0) = 
and X\ = X — Xq vanishing to infinite order at 0. Suppose also that they are both 
linear outside a compact set and equal there to their common linearization, DX(0), 
at which is assumed to have no pure imaginary eigenvalue. Let U(t), Uo(t) be the 
flows generated by X and Xq. If E is a linear submanifold invariant under Xq such 
that 

(3.10) lim U (t)x = V x € E 

t — >oo 

then for all j = 0, 1, 2, . . . and x G E 

(3.11) lim DHU{-t)U (t))x 

t — >oo 

exists, and is continuous in x G E, and 

W-x = lim U(-t)U (t)x, xeE 

t — >oo 

has a C°° extension, G, to M. N which is the identity to infinite order at and such 
that (G~ 1 ) :t X = Xq to infinite order along E in a neighbourhood ofO. 

Remark 3.8. Note that the derivatives in (|3.11| refer to the ambient space M. N , 
and not merely to E. This is useful in producing the Taylor series of G for the last 
part of the conclusion. 

Also, the limit t — ► oo means t — > +oo, as in Nelson's book. 

Proof. We follow the proof of Theorem 8 in [18] . Indeed, if Xq was assumed to be 
linear then Nelson's theorem would apply directly. Dropping this assumption has 
little effect on the proof; the main difference is that a little more work is required 
to show the exponential contraction property, (|3.12[) below. 

Since the real part of every eigenvalue of DX(0) is non-zero, M. N = E + £L 
where E + , resp. E_, is the direct sum of the generalized eigenspaces of DX(Q) with 
eigenvalues with positive, resp. negative, real parts. Since E is invariant under Xq, 
and hence under DX(0), necessarily E C £L. We actually apply the theorem with 
E = E-, but, as in Nelson's discussion, the more general case is useful for the 
inductive argument for the derivatives. 

Let Bj denote a basis of £L consisting of generalized eigenvectors of DX(0) 
with corresponding eigenvalue crj] we shall consider the Cj as differentials of lin- 
ear functions fj on R N . For x G R N , let x{t) = U (t)x, Fj(t) = fj(x(t)). Then 
^t\t=t = {Xofj)(x(t )) where 

X f,(y)=DX(0)f J ( y ) + O(\\y\\ 2 ). 

Moreover, for y G \\y\\ 2 < C\ J2j fj f° r some C\ > 0. So, setting p — J2 fji we 
deduce that 

X QP (y) = J2 2 ^fHv) + 0( P {yf /2 ), 

3 

hence with R(t) = p(x{t)), cq G (super,, 0), there exists 8 > such that for ||i?(f)|| < 
5, 

^-2c R<0, 
dt 

and hence R(t) < e _2cot ||a;|| for t > 0, ||r(a;)|| < 6, x G E— A corresponding- 
estimate also holds outside a compact set, as Xq is given by DX(0) there, so a 
patching argument and (|3.10[) yield the estimate R(t) < Coe _2c °'||a;|| for all x G 
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Since -R(t) 1 / 2 is equivalent to ||.||, we deduce that there are constants C, c > such 
that 

(3.12) |jt7" (t)ar|| < Ce" ct ||x|| V x £ E and t > 0. 

For the remainder of the argument we can follow Nelson's proof even more closely. 
Thus, let k be a Lipschitz constant for X and Xq, and choose m such that cm > k. 
Note that there exists cq > such that for all x £ M. N , 

(3.13) (x) || < collar, X 1 =X-X . 
For h > t 2 > 0, h = t 2 + t, x £ E, 

I = WUi-hpoit^x - U(-t 2 )U (t 2 )x\\ = \\U(-t 2 ) (U(-t)U (t) - Id) U (t 2 )x\\ 

< e Kt *\\(U(-t)U {t) -Id)U (t 2 )x\\ 

by the Lipschitz condition (see [HI Theorem 5] ) . But with X = Xq + X\, by [HI 
Proof of Theorem 6, (5)] 



||l/(-t)^o(*)v-y||< / e KS \\X 1 (U {s)y)\\ds. 
Jo 

Applying this with y = Ua(t 2 )x, we deduce that 

(3.14) I<e Kt2 f e KS \\X 1 (U {s + t 2 )x)\\ds. 

Jo 

Thus, by |3T3)) and (j3T2)) . 

I<e Kt2 f e KS c C m e- cm(s+t ^\\x\\ m ds 
Jo 

<e Kt2 / e KS c C m e- cm ( s+t ^ \\x\\ m ds = 







cm — k 

Letting t 2 —* oo shows that W_x = limt^oo U(—t)Uo(t)x exists, with convergence 
uniform on compact sets, hence W- is continuous in x £ E. Moreover, applying 
the estimate with t 2 = shows that W~(x) — x = C(||a;|| m ). Since m is arbitrary, 
as long as it is sufficiently large, this shows that W- is the identity to infinite order 
at 0, provided it is smooth, as we proceed to show. 

Smoothness can be seen by a similar argument, although we need to put a slight 
twist into Nelson's argument. Namely, first consider the first derivatives, or rather 
the 1-jet. Thus, we work on M. N £(R Ar ). Let (x, £) denote the components with 
respect to this decomposition. These evolve under the flow U'(t), resp. given 

by 

X'( X> = (X(x),DX(x)-0, X (x,0 = (X (x),DX (x)-0, 
where DX(x) and £ are considered as elements of £(R ), and • is composition 
of operators. Note that the second, £(R N ), component of these vector fields is a 
homogeneous degree zero vector field, i.e. it is invariant under push-forward by the 
natural R + -action (by dilations). 

The twist, as compared to Nelson's work, is that we identify £(R ) with M. N , 
which we radially compactify to a (closed) ball B N , which we further embed as 
the closed unit ball in M. N in such a fashion that the smooth structure of the ball 
agrees with the restriction of the smooth structure from M. N . Let l : l w — > 
be this map with range the interior of 5" . Then the push-forward under t of a 
homogeneous degree zero vector field, such as DX(x) ■£ is for each x £ R , extends 
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to a C°° vector field on the closed ball B N , which by homogeneity is tangent to the 
boundary. Furthermore, if l\ = idgiv x i, then and [l\)*X'q extend to C°° 

vector field on R N x B N tangent to the boundary and their difference, 
in addition vanishes to infinite order at {0} x B N . Thus (li)^X' and (ii)*Xq are 
Lipschitz with some Lipschitz constant k': this is automatic over a compact subset 
of M. N x B N , which in fact suffices here, but in fact holds on all of R N x B N since 
outside the inverse image of a compact subset of M. N x B N , X' and X' are linear, 
so in particular their B N component is independent of x. 

To minimize confusion about the 'change of coordinates', we write the coordi- 
nates on M. N x B N as (x, rf) below. With c as in (|3.f 2|) . choose m such that cm > k' . 
Then the infinite order vanishing of at x = yields 

(3-15) IK^JfOMH <c' Q \\x\r 

for all {x,rf). Let U'(t), U' {t) denote the evolution groups generated by (li)^X' 
and (li)*Xq, respectively. Thus, for all real t, 

(3-16) \\U'(t)(x,rj)\\ <e^\\(x, V )\\, 

see [TSl Theorem 5]. So (|3.1 41) still applies, with X\ replaced by {i\)*X' x , k replaced 
by k', etc. Thus, by (j3~12)) and (pUB]) . 

<e^ [ e K ' s ||(( il ),X 1 )(^( S + i 2 )(x,7 ? ))||d s 
Jo 

<e K't 2 f e K's c >c m e- cm{s+t2) \\x\\ m ds 



<e Kt2 / e KS c' n C m e- cm ^ +t2) \\x\\ m ds = ^—— ^ 



Thus, limt^oo U'(— t)U^(t)x exists, with convergence uniform on compact sets, so 
the limit depends continuously on (a;, £) for ie£. 

The higher derivatives can be handled similarly. The resulting Taylor series 
about E can be summed asymptotically, giving G: this part of the argument of 
Nelson is unchanged. □ 



3.3. Effective resonance and nonresonance. Next we apply this general result 
to the symbol p. Following Lemma 12. 7\ when resonances occur we cannot remove 
all error terms even in the sense of formal power series. Consequently we do not 
attempt to get a full normal form in a neighbourhood of the critical point, but only 
along the submanifold 

(3.17) S = {v = 0, y" = 0, y'" = 0, n = 0}, 

which is the unstable manifold for Wq ■ After reduction to normal form, errors which 
are polynomial in the normal directions to S will remain. For later purposes, we 
divide these into two parts. 
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Definition 3.9. With I as in Definition 13. 5( let 
I = I' U I" 

l' ex = {{a,a,J3)eI: a = (a\ a", a'"), (3 = (/?', /?", /?'"), 

(3.18) a = 0, a'" = 0, /?'" = 0, a" = 0, (3" = 0, |/3'| = 1}, 
l" x = {(a,a,(3) € I : a = (a', a", a'"), (3 = ((3',(3",P"'), 

a = 0, a'" = 0, 0" = 0, a' = 0, = 0}. 

An effectively resonant function is a polynomial of the form 

(a,a,/3)£/ el - 

or equivalently 

(a,a,/3)e/„ 

Thus, elements of I C r satisfy (0, a, /3) G J (i.e. are resonant, see Definition 13. 5|) . 
with a = (a'X',0), /3 = (/3',/3",0), and either a" = 0, /3" = 0, \(3'\ = 1, or a 1 = 0, 
/?' = 0. 

Moreover, an effectively resonant function has the form 

(3.19) j2 c*w'Vr'tfY+ E w(eT"(/'y"- 

a',\0'\ = l a", 13" 

For a fixed critical point of a fixed operator P (e.g. P = x _1 (A + V — a) for 
a fixed a), the set I eT is finite. Thus, only a finite number of terms can occur in 
(|3.19p . and hence restricting to polynomials in the definition of effectively resonant 
functions (rather than infinite formal sums) is in fact not a restriction. To see 
this, note that in the expression for R a , a ,p in f|3 . 8[> . we have a — 0, a'" — (3'" = 
and either (i) a" = (3" = and \(3'\ = 1 or (ii) a' = (3 1 = 0. In case (i), if 
(3j = 1 then to have R a , a ,i3 = we need X) a fe r fc = T 'j-> which is only possible for 
< kjl/mhifc \f'k\- in case (h), we need ^ct'-r 1 - + YlP'ji^ ~ r j) — 1' which is 
only possible for \a"\ < 1/ minr^' and \(3"\ < 2. (Actually in case (ii) we must have 

| < 1 in order to satisfy the condition 2a + \a\ + \j3\ > 3 in ([3S]).) 

Definition 3.10. Let Js denote the ideal of C°° functions on sc Tq X X which vanish 
on S and set 

!m-l 
{a",l3")-,Y,r'!a'!+{l-r';)% e{l,2) 
3= « 

An effectively nonresonant function is an element of Js of the form 
(3.21) r em = J2 h s f s + J2 K>,P>^'f*' + Y,W e '"fk 

3=1 (a",P")eI" j,k 

hj G Js, j — 0,1, ... , s, h'^ t0 „ G C°°rT* dx X), {a", 13") G 7", 

7$ G Js, j,k = m,...,n-l. 
Note that J" is finite, hence all sums in the definition arc finite. 
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Theorem 3.11. Using the notation of Lemma \2. 7| for coordinates near a radial 
point of q of p there is a local contact diffeomorphism $ from a neighbourhood of 
(0, 0, ... , 0) to a neighbourhood of q such that $*p = p n orm such that 

n-i 

(3.22) A -1 iWro = -V + ^ ''jlljl 1 , + ^2 QjiVj^j) + 7 'cnr + r eI , 

j j=m 

withr cnr of the form (|3 . 2 1 [) andr eT of the form (| 3 . 1 9 1) / in addition at anon-resonant 
critical point, i.e. if I = 0, then we may take r eia = r et = near q. 

Remark 3.12. If F is an elliptic Fourier integral operator with canonical relation $ 
then P = F~ 1 PF satisfies ag-x{P) — p n0 rm- 

Remark 3.13. As will be seen below, of the two error terms, only r er has any effect 
on the leading asymptotics of microlocal solutions. The construction below shows 
that modulo X°° , r cnr may be chosen to consist of resonant terms only, i.e. to be 
an asymptotic sum of resonant terms. However, this plays no role in the paper; all 
the relevant information is contained in the statement of the theorem. 

Remark 3.14. We do not need the full power of Lemma [2.71 to find $ as in this 
theorem; Lemma 12.101 suffices. Indeed, the terms Cj/J,^ in (|2.15[) can be 

absorbed in r cnr . 

Similarly any term v a yPy a with a + \f3\ > 2 and a ^ 0, or with \(3\ > 3 can be 
included in r er or r cnr . The same is true for any term with \(3\ > 2 such that f3j =/= 
for some j with Re Tj ^ \ . In particular, if Re Tj ^ | for any j, the only terms which 
need to be removed have a + \/3\ < 1. The conjugating Fourier integral operator can 
therefore also be arranged to have such terms only and thus to be of the form e lB , 
with B = Z + (f j x) where Z is a vector field on X tangent to its boundary and / 
is a real valued smooth function on X. Correspondingly, the normal form may be 
achieved by conjugation of P by an oscillatory function, e l ^ x , followed by pull-back 
by a local diffeomorphism of X, i.e. a change of coordinates. However, if Rer, = | 
for some j, some quadratic terms in fj, would also need to be removed for the model 
form, but since they play a role analogous to r or , the arguments of Section O giving 
conormality, are unaffected, and only the polyhomogeneous statements of Section [5] 
would need alterations. However, the contact diffeomorphism (i.e. FIO conjugation) 
approach we present here is both more unified and more concise. 

If p — |C| 2 + Vq — a, the model form of Lemma \2 . 1 01 also only required a change 
of coordinates and multiplication by an oscillatory function (see Lemma |2.11|) . the 
model form of this theorem can be obtained by these two operations, starting from 
the original operator P with symbol p. 

Proof. First we apply Proposition 13.61 Next we need to show that r or as in (|3.19|) 
and r en r as in (|3.19[) can be chosen to have Taylor series at given exactly by the 
error term in (|3 - 9[) . 

So, consider a monomial v a e a f^ with (a, a, (3) £ I. If a'" ^ then j3"' ^ since 
Imr™ > 0, and only the eigenvalues of have negative imaginary parts, and 
conversely. In addition, 2a + |a| + \(3\ > 3 implies that a monomial with a'" ^ or 
13"' ^ has the form v a e & f 'e'f '/£" for some j, k with 2a + \a\ + \f3\ > 1 and 



Re(o + 5^noi+5^(l-r l )A) = 0. 
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Since Re(l — n) > for all I and Ren > for I > s, while rj < for I < 
s — 1, we must have a' 7^ (i.e. a; 7^ for some I < s — 1) and correspondingly 
a + |a"| + \a"'\ + > 0. Due to the latter, v a e & f ^ vanishes on S, so the terms 
with a'" 7^ or /?"' 7^ appear in r onr . 

So we may assume that a'" — (3"' = 0. If a 7^ 0, the monomial is of the form 
v a e ajp v ^ & = a - 1, 25 + |a| + > 1 with 

2^+5^(1-^ = 0. 

Arguing as in the previous paragraph we deduce that the terms with a ^ also 
appear in r Gnr . 

So we may now assume that a = 0, a'" = j3"' = 0. If (3' 7^ 0, the monomial is of 
the form j/V/^/j for some j, and 2a + |a| + |/3| > 2, 

a + 2 r ; 5 ; + - n )Pi = rj < 0. 

We can still conclude that a' 7^ 0, but it is not automatic that a + \a"\ + \(3\ > 0. 
However, if a + \a"\ + \0\ > then v a e a f°fj is again included in r cnr , while if 
a + \a"\ + = 0, then the monomial is included in r er . 

Finally then, we may assume that a = 0, j3' = 0, a'" = /?'" = 0. Since r'j < for 
all j = 1, ...,s - 1 

+ f 1 - ^ E ^ + + a - = i- 

Moreover, the equality holds if and only if a' — 0, in which case this term is included 
in r cr . The terms with a' 7^ can be included in h" & „ ^,,e a fP for some a" < a" , 
/?" < /3", chosen by reducing a" and/or f3" to make 

EW + (i-^e(i,2). 

This can be done since r", 1 — r" 6 (0, 1). 

It follows that A _1 p can be conjugated to the form 

n-l 

(3.23) - v + 2 '•...'/,</'.. + 2 QiiVi'to) + rcnr + rcr + r °°' 

where r cnr ,r cr are as in (|3.21[) . (|3.19[) . with both vanishing if q is non-resonant, 
and Too vanishes to infinite order at (0, 0, 0). Thus, it remains to show that we can 
remove the term in a neighbourhood of the origin. 

To do this we apply Proposition l3.7l Let X' be the Legendre vector field of (|3.23|) . 
and let X[ be the Legendre vector field of , while X' = X' — X[ . Let X be the 
linear vector field with differential equal to DX(0), let x be compactly supported, 
identically 1 near 0, and let X = ~(xX' + (l- X )X), X = -(xX^ + ^-xjX). The 
overall minus sign is due to S being the unstable manifold of X' near the origin, 
hence the stable manifold of — X' Q . Let E be the subspace S of R 2 ™ -1 , defined by 
(|3.17[) . Then Proposition 13.71 is applicable, and G given by it may be chosen as a 
contact diffeomorphism since U(t), Uo(t) are such, see [H Section 3, Theorem 4]. □ 
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3.4. Parameter-dependent normal form. We also need a parameter-dependent 
version of this theorem. Namely, suppose that p depends smoothly on a parameter 
<7, can we make the normal form depend smoothly on a as well? This problem 
can be approached in at least two different ways. One can consider a simply as a 
parameter, so p £ C°°({d sc T*X) x I) = C°°(( BC T§ X X) x I) and then try to carry 
out the reduction to normal form uniformly. Alternatively, one identify p with the 
function p' on the larger space d sc T*(X x 7) arising by the pull-back under the 
natural projection 

p> = ** Pl tt : sc n XxI (X x /) - ( sc n x X) x / 

and then carry out the reduction to a model on the larger space. Whilst the second 
approach may be more natural from a geometric stance, we will adopt the first, 
since it is closer to the point of view of spectral theory of [5] . Clearly the difficulty 
in obtaining a uniform normal form is particularly acute near a value of a at which 
the effectively resonant terms do not vanish. Fortunately in the case of central 
interest here, and in other cases too, the set of points at which such problems arise 
is discrete. 

Lemma 3.15. If P = P(a) = x~ x (A + V — er), q — q(cr) is a radial point of P 
lying over the critical point z = 7t(q) ofVo and I(o~), resp. I er (a), are the sets (|3.8jl . 
resp. (|3.18p , for p{a) then the set TZ Z = lZnt,z U H er ,z, defined by 

^ Tlm,z = {o- £ (Vo(z), +oo) | 3j such that r,- = 

K cr , z = {ae (V (z), +oo) | I el (a) ? 0}, 

i.e., the set of energies a which are either a Hessian threshold (see Lemma \2.7\ ) or 
such that q{o~) has a nontrivial effectively resonant error term (see Definition \3.9\) . 
is discrete in (Vo(7r(g)), +oo). 

Remark 3.16. It follows that if K c (Vq(z), +oo) is compact then K n 7Z Z is finite. 
Thus, to prove properties such as asymptotic completeness, one can ignore all a € K 
which are Hessian thresholds or effectively resonant. 
Note also that by the definition of I er (o-), 

K el , z = {a e (V (z), +oo) | either 3 (0, (a 1 , 0, 0), (/?', 0, 0)) G 1(a) with \/3'\ = 1 

or 3 (0,(0,a",0),(0,/3",0))e/( ( 7)}. 

Proof. Using Remark 12.91 the set TZnt,z of Hessian thresholds is given by {Vb(z) + 
4a j} where aj is an eigenvalue of the Hessian of Vb at z and hence has cardinality 
at most n — 1, so this set is trivially discrete. 

Let K be a compact subset of (Vo(z),+oo). The set K n 7vL CT ,2 of effectively 
resonant energies in K is the the union of zeros of a finite number of analytic 
functions (none of which are identically zero). Indeed, lZ C r,z is given by the union 
of the set of zeros of the countable collection of functions 

tn— 1 s — 1 

- 1 + E < r » + ^ - r »)' -!+(!- »•*) + E 

as k = 1, . . . , s — 1, while a', a", j3" are multiindiccs. But if c > is large enough 
then c _1 > |r"j-(o-)| > c for all j and for all a £ K as K is compact and the rj do 
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not vanish there. Correspondingly, for \a'\ > 

s-1 

-1 + (1 - r k ) + o^-(ct) < -r k - \a'\c < -cT 1 , 
i=i 

and analogously for \a"\ + > |, 

m— 1 

- 1 + E < r » + - ? '») > 1 + CK'I + 1/ 3 "^ > L 

Thus, there are only a finite number of these analytic functions that may vanish in 
K, as claimed. □ 

If q(cr) are the radial points corresponding to z £ Cv(V), and a (fc TZct.z, then we 
will say that q(a) is effectively nonresonant, or that a is an effectively nonresonant 
energy for z. We now prove that, away from effectively resonant energies and 
Hessian thresholds, we have a normal form for p(a) of the form (|3 . 22[) with r or = 
and depending smoothly on a. Thus, for a given critical point z of Vq, consider an 
open interval O C (Vq(z), +oo) \ 1Z Z . Apart from the coefficients hj, h!' a ,, g„, etc., in 
(|3.21|) the only part of the model form depending on a is 

m — 1 

J" {a) = {(a", 13"); E r >K + C 1 - r »W e 

We note that on compact subsets K of O, there is a c > such that r"(a) > c for 
a eK, and then for \a"\ + \(3"\ > 2c" 1 , 

m— 1 

Sa „ „(a) = r»< + (1 - r'!(a))f3'J > 2, 

so if we let 

Jk = U ffe _ff J (a), 

then Jk is a finite set of multiindices. For each multiindex (a", (3") we let 

(3.25) O a „ t0 „ =s-^ ll ((l 7 2)), 

which is thus an open subset of O. 

For the parameter dependent version of the Theorem 13.111 we introduce 

(3.26) S = {(y, v, n, a); v = 0, y" = 0, y'" = 0, (i = 0, a € O}, 
in place of S (j3~T7j) . 

Theorem 3.17. Suppose that p e C°°( SC T^ X X x O), c (Vb(z), +oo) \ TZ Z 
is open, that the symplectic map S induced by the linearization A' of p at q{o~) 
(see Lemma \2.3\) can be smoothly decomposed (as a function of a £ O) into two- 
dimensional invariant symplectic subspaces and that there exists c > such that 
r'j((j) > c for a G O. Then $(cr) and F{a) can be chosen smoothly in a so that 
P„orm(cr) = CTi(P(cr)), P{a) = F (Vp 1 P{o)F (a) , is of the form in Theorem\ML 
with r or = 0, with the sum over J" replaced by a locally finite sum (the sum is over 
Jk over compact subsets K C O), the hj, etc., in (13.21)) depending smoothly on a, 
i.e. they are in C°°(^ c Tq X X x O), vanishing at S as in Theorem \3.11[ and with the 
h'a"p" supported in sc Tq X X x O a "p>> in terms of (13.251) . 
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Remark 3.18. For P — x~ 1 (A + V — a) the conditions of the theorem are satisfied 
for any bounded O — I disjoint from the discrete set of effectively resonant a, since 
in local coordinates (y,fJ.) on E(cr), the eigenspaces of S are independent of a as 
shown in the proof of Lemma [2~5l and the r'j are bounded below by Remark 12.91 

Proof. Since the invariant subspaces depend smoothly on a by assumption, so do 
the eigenvalues of the linearization, and there is smooth family of local contact 
diffeomorphisms, i.e. coordinate changes, under which p(cr) takes the form (|2.13| . 
i.e. 

m — 1 n— 1 

(3.27) p(a) = A(ct)( - v + ^ ^(cr)y^ + °-i( a > VhN) + v 9l + .92) 

3=1 j=m 

the Qj{o~, .), are homogeneous polynomials of degree 2, gi vanishes at least linearly 
and g2 to third order, all depending smoothly on a. 

For the rest of the argument it is convenient to reduce the size of the parameter 
set O as follows. For a G O, let 



0(0-) =[ n vV'((!> 2 )) 

(3.28) X 

\(a",|8»): »„«^W£(-«>,1) / 

an open set (as it is a finite intersection of open sets) that includes er. Thus, 
{0(a) : a G 0} is an open cover of O. Take a locally finite subcover and a 
partition of unity subordinate to it. ft suffices now to show the theorem for each 
element O(ag) of the subcover in place of O, for we can then paste together the 
models p n orm we thus obtain using the partition of unity. Thus, we may assume 
that O — O(a ) for some Gq € O, and prove the theorem with the sum over J" 
replaced by a sum over J"(cto). Hence, on O, for any (a" , j3") either 

a) 8 a »p»(<J Q ) > 1, and then for some (a" J") G J"(<tq), (a",/3") > (a" J") 
(reduce \a"\ + \/3"\ until S & „^„ G (1, 2) - this will happen as r 3 G (0, 1/2)) 

hence S a "/3"(a) > S (cr) > 1 for all a G O by the definition of O(ao), 
or 

b) s a "p"(ao) < 1, and then 5 a "p"(cr) < 1 for all a G O by the definition of 
O(a ). 

In order to make $(ct) smooth in a, we slightly modify the construction of the 
local contact diffeomorphism $i(<r) in Proposition 13.61 so that for any given a we 
do not necessarily remove every term we can (i.e. which are non-resonant for that 
particular a). Namely, we choose the set /' of multiindices (a,a,/3) which we do 
not remove by $i((j) so that V is independent of cr, and such that V contains 
every multiindex which is resonant for some a G O, i.e. /' D U - e o-^(' J )j with I(cr) 
denoting the set of multiindices corresponding to resonant terms for p(a), as in 
Proposition 13.61 With any such choice of J', the local contact diffeomorphism of 
Proposition 13. 6\ $1(0-), can be chosen smoothly in a such that A _1 $*p is of the 
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form 

m n—1 

-v+^2r j (a)y j n j + ^ Q 3 (CT,y J ,/ij)+^c QQ/ 3(cr)^ a e Q / /3 moduloX 00 = f)°° at q, 
j=l j=m+l r 

with c aa f3 depending smoothly on a. 

The requirement I' D U a eoI( a ) means that for (a, a, (3) I', R a ,a,/}(<j) must 
not vanish for a E O. Here we recall that R a ,a,/3(cr) is the eigenvalue of {{po, ■}} 
defined by (|3.5p . namely 

(n— 1 n—1 \ 

a - 1 + a J r i (<0 + J2 PiO- ~ r J j 

Keeping this in mind, we choose /' by defining its complement (I') c to consist 
of multiindices (a, a, j3) with 2a + \a\ + \/3\ > 3 such that either 

(i) a + \/3'\ = 1 and a" = 0, a'" = 0, f3" = 0, 0"' = 0, or 

(ii) | a'" | > 1, p" = 0, or 

(iii) | > 1, a'" = 0, or 

(iv) a = 0, p = 0, |a"'| + = 2, a" = 0, f3" = 0, or 

(v) a = 0, /?' = 0, a'" = /3"' = 0, Sa»p"{a) < 1 (for one, hence all, u e O, as 
remarked above). 

We next show that multiindices in (I') c are indeed non-resonant. In cases (|u])- 
(pUj). ImEa a ^u) =^ since the imaginary part of all terms in (|3.29p (with nonzero 
imaginary part) has the same sign, and there is at least one term with non-zero 
imaginary part, so (a, a, p) is non-resonant. 

In case (jvj) , the non- resonance follows from 

A _1 i?a : Q,/3(cr) < -l+s a "p"{<r) <0, 

since A" 1 R a ,a,f}{^) = — 1 +s a "/3"(cr) + 23<=i a j r ji an d each term in the last sum- 
mation is non-positive. 

In case 0, if a = 1, /?' = then A -1 R a ,a,p(&) = Y^jZi r j a j < since \a'\ > 1 
due to 2a + \a\ + \[3\ > 3. Also in case ([!]), if a = 0, = 1, with say /3i = 1, then 

s-l 

A _1 i? a , a ,/3(cr) = -n + ^ ayrj 
j=i 

which does not vanish since otherwise (a, a, /?) would be effectively resonant - it 
would correspond to one of the terms in the first summation in (|3.19[) . 
Finally, in case (jrv|) . 

s-l 

A -1 Rei? a , Qj/3 (cr) = otjrj < 

3=1 

since a' ^ due to 2a + |a| + \(3\ > 3. 

Thus, all terms corresponding to multiindices in (I') c can be removed from p(a) 
by a local contact diffcomorphism $i(cr) that is C°° in a. So we only need to remark 
that any term corresponding to a multiindex in I' can be absorbed into r enr (er). In 
fact, such a multiindex has either 

1) a+\/3'\ > 2, or 

2) a + |/3'| = 1 and \a"\ + \a'"\ + |/3"| + \(3'"\ > 1, or 
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3) \a"'\ + \(3"'\ > 3 (with neither a'" nor f3"' zero), or 

4) a = 0, p = 0, \a'"\ = 1, = 1, \a"\ + \[3"\ > 1, or 

5) a = 0, (3' = 0, a'" = 0, (3"' = 0, s a »^, > 1. 

The first two cases can be incorporated into the ho or hj terms in (|3.21[) . The third 
and fourth ones can be incorporated into the h'" k term. Finally, in the fifth case, 
any infinite linear combination of these monomials can be written as 

E ^(e'T''(/Y'\ 

(a",/3")eJ"(ao) 

as remarked in (a) after (|3.28[) . 
We thus obtain 

n-l 

A(cr)( - v + ^r J (a)y J p ] + ^ QjiVj^j) + r cnr (cr) + r^), 

j j=m 

with r enl as in (|3.2ip , and vanishes to infinite order at (0,0,0). Finally, we 
can remove the term in a neighbourhood of the origin, smoothly in a, using 
Proposition 13. 71 as in the proof of Theorem 13. Ill thus completing the proof of this 
theorem. □ 



4. MICROLOCAL SOLUTIONS 

In [5l Equation (0.15)] microlocally outgoing solutions were defined using the 
global function v on sc Tg X X. This is increasing along W and plays the role of a 
time function; microlocally incoming and outgoing solution are then determined 
by requiring the wave front set to lie on one side of a level surface of v. In the 
present study of microfocal operators, no such global function is available. However 
there are always microfocal analogues, denoted here by p, defined in appropriate 
neighbourhoods of a critical point. 

Lemma 4.1. There is a neighbourhood 0\ of q in sc Tg X X and a function p G 
C°°(Oi) such that (D\ contains no radial point of P except q, p{q) — 0, and Wp > 
onY.nO with Wp > on E n O x \ {q}. 

Proof. This follows by considering the linearization of W. Namely, if P is conjugated 
to the form (|2.13p . then for outgoing radial points q take p = \y'\ 2 — (\y"\ 2 + \y"'\ 2 + 
|/i| 2 ), defined in a coordinate neighbourhood Oq, for incoming radial points take 
its negative. On E, Wp > c(\y\ 2 + \p\ 2 ) + h for some c > and h € I 3 . As (y,p) 
form a coordinate system on E near q, it follows that Wp > |(|?/| 2 + \p\ 2 ) on a 
neighbourhood O' of q in E. Now let O x C O be such that 6 n E = O' . Then 
Wp(p) = 0, p € 0i, implies p = g, so there are indeed no other radial points in Oi, 
finishing the proof. □ 

Remark 4.2. Below it is convenient to replace Oi by a smaller neighbourhood O of 
q with O C Oi, so p is defined and increasing on a neighbourhood of 0. 

Consider the structure of the dynamics of W in 0. First, p is increasing (i.e. 
'non-decreasing') along integral curves 7 of W, and it is strictly increasing unless 7 
reduces to q. Moreover, W has no non-trivial periodic orbits and 

Lemma 4.3. Let be as m Remark\J^ If 7 : [0, T) -> or 7 : [0, +00) -> is 
a maximally forward-extended bicharacteristic, then either 7 is defined on [0, +00) 
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and limt— >+oo j(t) = q, or 7 is defined on [0, T) and leaves every compact subset K 
of O, i.e. there is Tq < T such that for t > To, l{t) £ K. 

An analogous conclusion holds for maximally backward- extended ^characteris- 
tics. 

Proof. If 7 : [0, +00) — > O then lim^+oo p{j{t)) — p+ exists by the monotonicity 
of p, and any sequence jk '■ [0, 1] — > E, "fk(t) = ^{tk + t), tk — > +00, has a uniformly 
convergent subsequence, which is then an integral curve 7 of W in E with image in 
O, hence in (D\ along which p is constant. The only such bicharacteristic segment 
is the one with image {g}, so lim^ +00 7(f) = q. The claim for 7 defined on [0, T) 
is standard. □ 

As in [5] we make use of open neighbourhoods of the critical points which are 
well-behaved in terms of W. 

Definition 4.4. By a PF-balanced neighbourhood of a non-degenerate radial point 
q we shall mean a neighbourhood, O, of q in sc Tg X X with O C O (in which p is 
defined) such that O contains no other radial point, meets E(er) n O in a MK-convex 
set (that is, each integral curve of W meets E(cr) in a single interval, possibly empty) 
and is such that the closure of each integral curve of W in O meets p — p(q). 

The existence of VF-balanced neighbourhoods follows as in Lemma 1.8]. 
If q is a radial point for P and O a VF-balanced neighbourhood of q we set 

(4.1) E mi c,+ (0, P) = {ue C-°°(X); O n WF sc (Pu) = 0, 

and WF sc (u)DO C {p> p(q)}}, 
with £ ; m i C! _(0, P) defined by reversing the inequality. 

Lemma 4.5. If O 3 q is a W -balanced neighbourhood then every u £ E mlCi ±(0, P) 
satisfies WF sc (u) nOc §±({q})] furthermore, for u £ E miCj ±(0, P) 

WF SC («) no=«^^ WF sc (u). 

Thus, we could have defined -E m ic,±(0, P) by strengthening the restriction on the 
wavefront set to WF sc («)nO C §±({q})- With such a definition there is no need for 
O to be iy-balanced; the only relevant bicharacteristics would be those contained 
in $±({g}). Moreover, with this definition p does not play any role in the definition, 
so it is clearly independent of the choice of p. 

Proof. For the sake of defmiteness consider u £ E m - lc , + (0, P); the other case follows 
similarly. Suppose ( £ O \ {q}. If p(() < p(q), then ( £ WF sc (m) by the definition 
of £' m i Ci+ (0, P), so we may suppose that p(() > p(q). Since q £ $+({q}) we may 
also suppose that £ 7^ q. 

Let 7 : R — » E be the bicharacteristic through ( with 7(0) = £■ As O is W- 
convex, and WF sc (Pu) n O = 0, the analogue here of Hbrmander's theorem on the 
propagation of singularities shows that 

C € WF sc (u) 7 (M) HOC WF SC («). 

As O is W^-balanced, there exists (' £ 7(E) n O such that p((') = p(q). If p(() = 
p(q) — 0, we may assume that £' = £. From this assumption, and the fact that p 
is increasing along the segment of 7 in O, and O is M^-convex, we conclude that 

C'e 7 ((-oo,0])nO. 
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If £' = j (to) for some to G R, then for t < to, p(l(t)) < p(j(to)) — p(q), and 
for sufficiently small \t — to\, 7(t) G O as O is open. Thus, 7(f) G" WF SC («) by the 
definition of £ , m j Cj+ (0, P), and hence we deduce that £ G" WF sc (u). 

On the other hand, if Cf ^ 7(K), then as O is open 7(ifc) G O for a sequence 
tfe — > —00, and as O is W-convex, 7|(_oo,o] C O. Then, again from the propagation 
of singularities and Lemma l4~3l Cf = q. □ 

We may consider E m i c ,±(0,P) as a space of microfunctions, E m i Ci+ (q,P), by 
identifying elements which differ by functions with wavefront set not meeting O: 

Emic,±(<l,P) = E m ^ ± (0,P)/{u e c-°°(xy, WF SC ( U ) n o = 0}. 

The result is then independent of the choice of O, as we show presently. 
If Oi and O2 are two IF-balanced neighbourhoods of q then 

(4.2) Ox C 2 Pmic,±(0 2 ,P) C Ba.icifOl.P). 

Since {w G C _00 (X); WF sc (u) n O = 0} C E mic ,±(0,P) for all O and this linear 
space decreases with O, the inclusions (|4.2j) induce similar maps on the quotients 



Pmic.±(0, P) = Pmic,±(0, P)/{u G C- x (X);WF sc (u) n O = 0}, 
OiC0 2 ^ P mi c,±(0 2 ,P) — ► P mic , ± (0 1; P). 

Lemma 4.6. Provided Oi, for i = 1, 2, are W -balanced neighbourhoods of q, the 
map in (|4.3|) is an isomorphism. 



Proof. We work with £ , m i c ,+ for the sake of definiteness. 

The map in (|4.3|l is injective since any element u of its kernel has a represen- 
tative u G Pmic,+ (0 2 , a) which satisfies q G' WF sc (u), hence WF sc (u) n O2 = by 
Lemma l4~5l so u = in £' m i Ci+ (02, &)■ 

The surjectivity follows from Hormander's existence theorem in the real principal 
type region [13]. First, note that 

R = inf{p(p) : p G $+({«}) n (O \ Ox)} > - p(q) 

since in O, p is increasing along integral curves of W, and strictly increasing away 
from q. Let U be a neighbourhood of 4> + ({q}) n Oi such that U C O, and p > Rq = 
R/2 on f/\ Oi. Let 4 G *- c oo ' (O) be such that WFg C (Id -A) nOifl ®+({q}) = 
and WF' SC (A) C C/. Thus, WF sc (ylu) C U and WF sc (PAw) C U\O u so in particular 
p > Ro on WF sc (P^4u). We have thus found an element, namely u — Au, of the 
equivalence class of u with wave front set in O and such that p > Ro > = p(q) 
on the wave front set of the 'error', Pu. 

The forward bicharacteristic segments from U \ 0\ inside O leave O2 by the 
remark after Lemma I4TT1 since O2 \ 0\ is compact, there is an upper bound T > 
for when this happens. Thus, Hormander's existence theorem allows us to solve 
Pv = Pu on O2 with WF sc (u) a subset of the forward bicharacteristic segments 
emanating from U \ 0%. Then u' = u — v satisfies WF sc (u') C O n {p > = 
p(q)}, WF sc (Pu') n 2 = 0, so v! G P mic ,+ (0 2 ,P), and q £ WF SC («' - u). Thus 
WF SC («' — u) n 0\ = 0, hence u and it' are equivalent in E m i C + (Ox, P). This shows 
surjectivity. □ 

It follows from this Lemma that the quotient space E m i Ct ±(q, P) in (|4.3|) is well- 
defined, as the notation already indicates, and each element is determined by the 
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behaviour microlocally 'at' q. When P is the operator x 1 (A + V — a), then we 
will denote this space 

(4.4) £mic,±(<Z,<x). 

Definition 4.7. By a microlocally outgoing solution to Pu = at a radial point q 
we shall mean either an element of -E m ; Cj +(0, P), where O is a W^-balanced neigh- 
borhood of q, or of £ , m ic,+ (<Z, P)- 

5. Test Modules 

Following [5], we use test modules of pseudodifferential operators to analyze 
the regularity of microlocally incoming solutions near radial points. This involves 
microlocalizing near the critical point with errors which are well placed relative 
to the flow. For readers comparing this discussion to [5], we mention that the 
microlocalizer Q in the following definition corresponds to the microlocalizer Q in 
[5j Equation (6.27)]; the orders in the commutator are different as now we are 
working with P S 1 ^*^ 1 (X). 

Definition 5.1. An element Q G ^*f(X) is a forward microlocalizer in a neigh- 
bourhood O 3 qoia radial point q G sc Tg X X for P G if it is elliptic at q 
and there exist B, F € ^ s i°(O) and G € ^(X) such that 

(5.1) i[Q*Q, P] = (B*B + G)+F and WF' SC (F) n $+({<?}) = 0- 

Using the normal form established earlier we can show that such forward mi- 
crolocalizers exist under our standing assumption that 

the linearization has neither a Hessian threshold subspace, (IrvT) . 

(5.2) 

nor any non-decomposable 4-dimensional invariant subspace. 

Proposition 5.2. A forward microlocalizer exists in any neighbourhood of any 
non-degenerate outgoing radial point q € sc Tg X X for P G , F* ( ;~ 1 (X) at which the 
linearization satisfies (|5.2| . 

Proof. Since the conditions (|5.ip are microlocal and invariant under conjugation 
with an elliptic Fourier integral operator, it suffices to consider the model form in 
Theorem 13.111 which holds under the same conditions (|5.2[) . 
Let R = \n'\ 2 + \y"\ 2 + \y"'\ 2 + | M "| 2 + \fx"'\ 2 , and 

S = {Pnorm = 0, R = 0}, 

so S is the flow-out of q. We shall choose Q G *~ oo ' (X) such that 

(5-3) a d (Q) = q = Xl(|2/T)X2 Wfcnorm), 

where \ii X2, 4> S C^?°(R), xii X2 > are supported near 0, ip is supported near 
0, xii X2 = 1 near and < in [0,oo). Choosing all supports sufficiently small 
ensures that Q G ^ c °°< (O). Note that suppd(x2 R) H S* = 0. On the other hand, 

(5.4) sc ^Xi(E(^) 2 ) = 2 Z>i( SC ^^)Xi(ly'| 2 ) - ^0''y[ + h^dyf), 

j 3 

with hj vanishing quadratically at q. Moreover, on suppx'i (|-| )•> y' is bounded 
away from 0. Since r'j < 0, —J2j r j( y 'j) 2 > on suppxi (|-| 2 )- The error terms 
hj can be estimated in terms of \y'\ 2 , R and p 2 lorm , so, given any C > 0, there 
exists 5 > such that the ~J2jyj( r 'jyj + hj) > if suppxi C (—6,8), R/\y'\ 2 < 
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C and |fWm|/|y'| < C. In particular, taking C = 2, - Y^jV'ji^jyj + h j) > 
on S n suppxi ° (|-| 2 ), for R — Pnorm = on S. Thus (|5.1[) is satisfied (with P 
appropriately specified, microsupported near S), provided that xi is chosen so that 
(— XiXi) 1 ^ 2 is smooth. 

More explicitly, letting \ G C£°(R) be supported in (—1, 1) be identically equal 
to 1 in (—5, 5) with x' < on [0, 00), x > 0> Xi = X2 — 4> = x(-/$)- Indeed, for any 
choice of Se (0,1), |yf > 5/2 on suppxi° |.| 2 , hence R/\y'\ 2 < 2, |p norm |/|y'| < 2 
on suppg n suppxi |-| 2 - With C = 2, choosing <5 G (0, 1) as above, we can write 

a d (i[Q*Q, P}) = - sc H p q 2 = -4A6 2 + /, 

(5 ' 5) (E^^^ + ^)x'l(|y'| 2 )Xl(|y'| 2 )) 1/2 X2(i?)^(Pno rm ), BU P p/n5 = 0, 

i 

which finishes the proof since A < for an outgoing radial point. □ 

A test module in an open set O C sc Tq X X is, by definition, a linear subspace 
M. C ^I'* ( ;~ 1 (X) consisting of operators microsupported in O which contains and is 
a module over (X), is closed under commutators, and is algebraically finitely 
generated. To deduce regularity results we need extra conditions relating the mod- 
ule to the operator P. 

Definition 5.3. If P G ^*^ 1 (X) has real principal symbol near a non-degenerate 
outgoing radial point q then a test module M is said to be P-positive at q if it is 
supported in a M^-balanced neighbourhood of q and 

(i) M is generated by A = Id, A 1} . . . , A N = P over ^*^(X), 

(ii) for 1 < % < N - I, < j < N there exists Cy G ^^(X), such that 

N 

(5.6) i[A i ,xP]=J2xC ij A j 

3=0 

where (Xg(Cij)(q) = 0, for all 7^ j < i, and Recrg(Cjj)(q) > 0. 

As shown in [5] , microlocal regularity of solutions of a pseudodifferential equation 
can be deduced by combining such a P-positive test module with a microlocalizing 
operator as discussed above. We recall and slightly modify this result. 

Proposition 5.4. (Essentially Proposition 6.7 o/[5] ; see Provosition \A. H below for 
a slightly modified statement and a corrected proof). Suppose that P G 1 ^*^~ 1 (X) 
has real principal symbol, q is a non- degenerate outgoing radial point for P, 

(5.7) a dtl (xP-(xP)*)(q) = 0, 

M is a P-positive test module at q, Q,Q' G ^^(A") are forward microlocalizers 
for P at q with WF' SC (Q') being a subset of the elliptic set of Q. Finally suppose 
that for some s < — ^, u G H^1' S (X) satisfies 

(5.8) WF sc (u) HOC "MM) and Pu G C*°°(A). 
Then u G Isc\0',M) where O' is the elliptic set of Q' . 

Proof. As already noted this is essentially Proposition 6.7 of [5] , with a small change 
to the statement and the proof given in Proposition I A. 1 1 below. However, there are 
some small differences to be noted. In [5] (and here in the Appendix), the condition 
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in (|5.6|) was j > i; here we changed to j < i for a more convenient ordering. Since 
the labelling is arbitrary, this does not affect the proof of the Proposition. 

Also, in [5] the proposition was stated for the Oth order operators such as A + V — 
a, which are formally self-adjoint with respect to a scattering metric. This explains 
the appearance oixP both in (|5.7[) and in (|5.6p here, even though in the applications 
below, [Aj,af] could be absorbed in the Cm term. In particular, s < —1/2 in (|5.8j) 
arises from a pairing argument that uses the formal self-adjointness of xP, modulo 
terms that can be estimated by [x s A a , xP], s > 0, A a a product of the Aj. 

The proposition in [5] is proved with (|5.7[) replaced by (xP) = (xP)* , but (|5.7p is 
sufficient for all arguments to go through, since B = (xP) — (xP)* would contribute 
error terms of the form x s A a P> with ag t i(B)(q) = 0, which can thus be handled 
exactly the same way as the Cjj term in (|5.6|) . 

In fact (|5.7|) can always be arranged for any Po £ \l/* ( i _1 (X) with a non-degenerate 
radial point and real principal symbol. Indeed, we only need to conjugate by x k 
giving 

P = x k P x \ k = Z^ABM e H 
2i\ 

satisfies (|5.7I) ; here dp\ q = Xa\ q , with a the contact form. Microlocal solutions 
-Po^o = 0, correspond to microlocal solutions Pu = via u — x k UQ, so u £ H^' S (X) 
is replaced by u G H^' s ~ k (X). □ 

Thus, iterative regularity with respect to the module essentially reduces to show- 
ing that the positive commutator estimates (|5.6|) hold. For each critical point g 
satisfying (|5.2p a suitable (essentially maximal) module is constructed below, so 
microlocally outgoing solutions to Pu — have iterative regularity under the mod- 
ule; that is, that 

(5.9) u G I^(0,M) = {u;M rn u C H^ S (X) for all m}. 

The test modules are elliptic off the forward flow out which is an isotropic 

submanifold of E. Thus, it is natural to expect that u is some sort of an isotropic 
distribution. In fact the flow out (in the model setting just the submanifold S) has 
non-standard homogeneity structure, so these distributions are more reasonably 
called 'anisotropic'. 

First we construct a test module for the model operator when there are no 
resonant terms. Thus, we can assume that the principal symbol is 

m— 1 n— 1 

po = A( - v + TjVito + QiiViiH))- 

3=1 3=171 

Then let M. be the test module generated by Id and operators with principal 
symbols 

(5.10) x' 1 ^, x-^'e'!, x-W>f!, x' 1 ^', x^ 2 ^ and x^po 
over 9%°(X). 

Note that the order of the generators is given by the negative of the normalized 
eigenvalue (i.e. the eigenvalue in Lemma [2771 divided by A) subject to the conditions 
that if the order would be < —1, it is adjusted to —1, and if it would be > 0, it 
is omitted. The latter restrictions conform to our definition of a test module, in 
which all terms of order are included and there are no terms of order less than — 1. 
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These orders can be seen to be optimal (i.e. most negative) by a principal symbol 
calculation) of the commutator with A in which the corresponding eigenvalue arises. 

Lemma 5.5. Suppose P is nonresonant at q. Then the module M generated by 
(|5.10D is closed under commutators and satisfies condition (|5.6D . 



Proof. It suffices to check the commutators of generators to show that M. is closed. 
In view of (|2.4|) (applied with a in place of p), {a,b} = sc H a b, this can be easily 
done. Property (|5.6p follows readily from (|3.ip . Indeed, we have the stronger 
property 

i[Ai,P(a)] = CiAi + Gi, G l e **'°(X), Rec, > 
where A{ is any of the generators of M. listed in (|5.10p . □ 

Remark 5.6. We may take generators of M to be the operators 

D ylj , x'^y'J, x r iD y », x- 1/2 y '?, x 1/2 D y ,„ and 
(5 11) m_1 n_1 

j=l j=m 

Combining this with Proposition [53] proves that, in the nonresonant case, if u is 
a microlocal solution at q, and if WFg C (ii) is a subset of the IF-flowout of q, then 
u e I$(0,M) for all s < -1/2. 

The discrepancy between the 'resonance order' of polynomials in v a e^ f 1 , given 
by a +J2j Pj r j+J2k 7/c(l~ r fe) an d the 'module order' given by the sum of the orders 
of the corresponding module elements is closely related to arguments which allow 
us to regard most resonant terms as 'effectively nonresonant'. To give an explicit 
example, take a resonant term of the form y'iUjiy")^ , corresponding to a term like 
x~ 1 y' i (y"Y [xDyi.) in P. Resonance requires that r- + (1 — r'j) + J2k^k r k = 1 an d 

| > 0. In the module, this corresponds to a product of module elements with an 
additional factor of x t with e > 0, since we can write it 

x%Y[(x-^yi'f*D v ,, e = J2^ r 'k>°- 

k k 

Since, by Proposition ^. 4[ the eigenfunction u remains in x s L 2 (X), for all s < —1/2, 
under application of products of elements of M. , this term applied to u yields a factor 
x t , and therefore it can be treated as an error term in determining the asymptotic 
expansion of u; see the proof of Theorem 16.71 Only the terms with the module 
order equal to the resonance order affect the expansion of u to leading order, and 
it is these we have labelled 'effectively resonant'. 

Next we consider the general resonant case. To do so, we need to enlarge the 
module M. so that certain products of the generators of M. , such as those in the 
resonant terms of Theorem 13.111 are also included in the larger module M. . For a 
simple example, see section 8 of Part I. It is convenient to replace Pq by xD x as 
the last generator of M. listed in (|5.11|) , though this is not necessary; all arguments 
below can be easily modified if this is not done. Let us denote the generators of M. 

by 

A a = Id 1 A 1 =x- Sl B 1 ,...,A N ^ 1 =x- SN - 1 B N - 1 , A N =xD x = x - 1 B N , 

5.12 

Si = - order (At), B t € *~ °°' O). 
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Note that for each i — 1, . . . , N, dagfl(Bi) is an eigenvector of the linearization of 
W; we denote the eigenvalue by ai. Thus, 

Si — min(l, <7j) > for i = 1, . . . , N. 

For any multiindex a £ N N (with N = {1, 2, . . .}) let 

s(a) = min(\J s^, 1), s(a) — s^a,; — s(a) — max(0, SjQij — 1), 

i i 

and let 

tot _ \u\ \ui ,a N 
rt. — J L L 1 A 2 ■ ■ ■ ^-N ■ 

Let ei be the multiindex = (0, . . . , 0, 1, 0, . . . , 0), where the 1 is in the «th slot, if 
i = 1, . . . ,N, and let e = (0, . . . , 0). 

To deal with resonant terms, we define a module M-k generated (over ^f~ c °°' (O)) 
by the operators 

(5.13) x s ^A a eV-°°'~' ( - a HO), \a\<k. 

Note that a = gives Id as one of the generators. Thus, the order of the gener- 
ators in (|5.13p is 'truncated' so that it is always between and —1; in particular 
Mk C , 3>~ 00 ' _1 (0). Since in computations below we will think of , 3>~ oo ' (O) as the 
submodule of M.k consisting of trivial elements, it is convenient to work modulo 
such terms, so we use what is essentially the principal symbol equivalence relation 
on M k where P ~ Q if P - Q G *-°°'°(0). 

While it appears that the ordering in the factors in the product A a matters, this 
is not the case. Indeed, if a is a permutation of {1, . . . , |a|}, and j : {1, . . . , \a\} i— > 
{1, . . . , N} which takes a m -times the value m, m = 1, . . . , N, then 

for this is clear if a interchanges n and n+ 1, as 

X A j(.l) ■ ■ ■ A j{n-l)[ A j{n), A j{n+l)] A j{n+2) ...A, (H) 

G ^-^^W+i-E^^^o) c *-°°' Q (0) 

since s(a) + 1 — Si«j = 1 — s(a) > 0. 
In addition, for Q G *-°°' (O), 

£ s(a) (2A,-(i) . . . ~ x s(a) A,-(i) . . . Aj( m }QAj( m+ x) . . . AjQ a \y 

Similarly, one can shift powers of x from in front of the product to in between 
factors, so in fact the generators can be written equivalently, modulo ^~ c °° (O), as 

(5.14) x s ^B a G *-°°'~ s M(0), \a\ < k, 

where B a = B* 1 ...B% N . 

Moreover, there is an integer J such that M.k = M.J if k > J; indeed this is true 
for any J > 2(r")~ 1 , where r" is the smallest positive eigenvalue of the operator in 
Lemma l2.5l for J > 4 if no eigenvalue lies in (0, i]), since then adding new elements 
to the product simply has the effect of multiplying by an element of '^lf(X). 

In particular, note that the generators in (|5.13[) or (|5.14j) are usually not linearly 
independent: some B aj may be absorbable into a , 3>* ( ; (O) factor without affecting 
s(a). We could easily give a linearly independent (over ^*f(0)) subset of the 
generators, but this is of no importance here. 
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Suppose that P, the normal operator for P(cr) at q, contains resonant terms. 
Then Lemma 15.51 is replaced by 

Lemma 5.7. Let > be a total order on multiindices a satisfying 

(i) | a' | > |a| implies a' > a; 

(ii) | a' | = |a| and J2k s k a 'k > Efc s k a k imply a' > a; 

(hi) \a'\ — \a\ = 1, a' = ej, a = ej, Si = Sj — 1, Oi > o-j imply that a' > a. 

With the corresponding ordering of the generators x~ s ^ A a , the module Mj is a 
test module for P at q satisfying ()5.6() . 

Remark 5.8. (ii) and (iii) could be replaced by (ii)': \a'\ — \a\ and J2 k aka 'k > 
Y], o~kCtk imply a' > a, which would simplify the statement of the lemma. However, 
the proof is slightly simpler with the present statement. Note that (ii)+(iii) is not 
equivalent to (ii)', i.e. the ordering of the generators may be different, but either 
ordering gives (15. 6|) . 

Proof. We first observe that Mj is closed under commutators. Indeed, not only 
is A4 closed under commutators, but the commutators [Ai,Aj] can be written as 

Eilo °l A l with °l e ^so°°'°( X ) and °i = unlcss s i < s i + s j - !■ Expanding 

[x s ^Q a A a ,x s ^Q A% Q ai Qp e *-°°' (O), 

and ignoring momentarily the commutators with powers of x and with Q a and Qp, 
gives a sum of terms of the form 

X sw+W Q a Q p A a ' A ' [Ai , Aj]A a " A?" 

with a = a' + a" + e^, and similarly for (3. Substituting in [Ai, Aj] = E;=o CiAi 
shows that this is an element of the module and is indeed equivalent, modulo 
*-°°'°(0), to 

(C ; :/ (Q)+S( ' 3) ~ S ~ (7<!)) ) x S ^ (l)) A' <(l) , 

(5.15) ;:Si< Si +Sj-l 

7 W = a' + a" +{3' + {3" + ei = a + /3-ei- ej + e h 

provided that 

(5.16) s( 7 w ) < s(a) +s(P). 

Buta(a)+S(/3) > (E^-1) + (E^A-1) = E^+Si+Sj-s^ > E*i7i° - 1 
as Si + Sj — si > 1. Moreover, s(a) + S(f3) > 0, so 

s(a) + 5(/3) > max(J] s fc7 f - 1, 0) = s( 7 «), 

proving (|5. 16[) . 

The commutators 

(5.17) x s ^Q p [x s ^Q a ,A^]A a , x s ^Q a [A a ,x s ^Q ]A' 3 

also lie in Mj. Indeed, [Ai,afQ] = x p - s ' +1 Q' for some Q' £ *^ oc * (O), so they 
are sums of terms of the form x s ^ a ' +s ^' Si+1 Q' A 1 with 7 = a + /3 — ej. Now, 

s( 7 ) < s(a) + s(/3) - Si + 1 
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since s(a)+s(/3) — s» + l > as 1 > Si as well as s(a)+s((3) — Sj + 1 > {J2k s fcCkfc — 1) + 
{J2k s kPk-l)-Si + l = J2 k Sfc7fc -1, so s(a)+s(/3)-Si + l > max(J] fc Sfc7fc-1,0) = 
5(7) indeed, proving that (|5.17|) is in M.j. The commutators 

(5.18) [x?WQ a ,x?WQ fl ]A a A< 1 

can be shown to lie in Mj by a similar argument, this time using 7 = a + /3, and 
1(7) < s(a) + §{f3) + 1. Thus, we conclude that [x~ s{ -^Q a A a , x^QpA 13 } e Mj, 
and hence A4j = M. j+i = ... is closed under commutators. 

Modulo *-°°' (O), a; 5 ^' ))^' may be replaced by a;" 8 ^ )^* . If | 7 (0| > J 
in (|5.15p . then this is written in terms of one of the generators listed in (|5. 14|) 
(or equivalently, modulo \t~°° ,o (0), in (|5.13p ). only after some of the factors in 

B~< ( , which we may always take from B\B$ B 13 , are moved to the front and are 
incorporated in Ci, i.e. they are simply regarded as Oth order operators and C; is 
replaced by Ci = CiBiB 13 B 13 . Notice the principal symbol of Ci always vanishes 
at q in this case. Analogous conclusions hold for the terms in (|5.17|) and (|5.18p . 

On the other hand, if | 7 ^| < J, then x s ^ 7< ^A 7 ' ' is one of the generators in 
(T5TT3T) . and = \a\ + |/3| - 1 if I > 1, and | 7 ( ; )| = \a\ + |/?| - 2 if I = 0. Moreover, 
if J2k s kPk > 1 then 
(5.19) 

^ Sfc7fe ^ = ^Skak + ^SkPk-Si-Sj+si > s fc Q fc + s k /3k - 1 > ^Sfcafc- 

For the terms in (|5.17jl and (|5.18p . if I7I < J, we always get I7I > |a| + \/3\ — 1 since 
7 = a + /3or7 = a + /3 — &i for some i. 

Now we turn to (|5.6p . First, with P replaced by Pq, (|5.6[) is certainly satisfied, 
exactly as in the non-resonant case, since the (Jd,o{B a ) are eigenvectors of the 
linearization of W with eigenvalue given in Section [3] Thus, 

(5.20) i[A a lX - l P ] ]T C\A\ C; e *«°°' (O), 

7 

with UdfiiC'^iq)) = if a 7^ 7 and Re Vd,o{C' a {q)) > 0. So it remains to show 
that it also holds for the resonant terms. If x^^^QpB 13 is a resonant term, then 
s{(3) — 1. Moreover, 

(i) if \(3\ = 1, then x~ 1 QpB l3 = J^^'iv')^ ^y' k f° r some A*' an d some fc; in 
particular it is a summand of r er ; 

(ii) if \(3\ = 2, then either x^QpB 13 = B 3 D y < k for some j > 0, fc, or x~ 1 Q f3 B /3 
is associated to the sum over J" in (|3.21|) : in cither case s kPk > 1. 

We claim that for a resonant term x~ s ^QpB^ , 

(5.21) [a;- s ^B Q ,x- s( ' 3) Q /5 B' 3 ] ~ ^ C^-^S 7 , <? 7 S l-" 00 ' ^), 

7 

and each term on the right hand side has the following property: 

(i) Either a g , (CJ (g) = 0, 

(ii) or | 7 | > j a I, 

(iii) or |7| = \a\, J^k s klk > E fc s fca fc , 

(iv) or I7I = |a| = 1, 7 = e k , a = ej, Sj = s k = 1 and a k > <J r 

Indeed, if |/3| > 3, then either (i) or (ii) holds, depending on whether any factors 
Ak had to be cancelled to write the commutator in terms of the generators in 
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(|5.13|) . If \(3\ = 2, then ^SkPk > 1. Thus, again, either (i) or (ii) holds, or 
|T I = M an d J2k s fc7fe > J2k s k a k by (|5.19|) . so (hi) holds. Finally, if |/3| = 1, then 
x~ 1 QpE>P — Ylu'iv')^ D y > k for some // and some k. Since n < r 2 < . . . < r s _i < 0, 
and the resonance condition is ^f =1 jujrj + (1 — 7"fc) = 1 with |//| + 1 > 3, we 
immediately deduce that /xj = for Z < k. Thus, not only do powers of x commute 
with x~ 1 QpBP , but all A commute with D y ' k and [Aj, (y'Y ] = unless A; = _D y ' 
and /Uj ^ for some j, which in turn implies that j > k, so 1 — > 1 — rj, hence 
(iv) holds. This completes the proof of (|5.21|) . 

By the assumption on the ordering of the multiindices a, we deduce that for all 
resonant terms x~ s 

7 

and either ag i o(C 1 )(q) — 0, or 7 > a. Combining this with (|5.20|) , we deduce that 
M.j satishes (|5.6[) . This establishes the lemma. □ 



Corollary 5.9. Let M — Mj be as in the previous lemma. Suppose that 

(5.22) s <--, ue H™' S (X), Pu £ C°°(X) and WF sc (u) nOc $+({?})• 

Then u e 1^(0, M). 

Regularity with respect to M. can be understood more geometrically as follows. 
Suppose 5 > is sufficiently small so that (x, y', y" , y") define local coordinates on 
the region U given by < x < 8, \y -\ < 6 for all j. Let 

y" y'" 

(5.23) * : -R», <&(x, y' , y" , y"') = (x, y' , Y", Y'"), Yf = Yj" = 

x 3 x L i z 

Thus, $ is a diffeomorphism onto its range <J>(£/°) with 

$-\x,y',Y",Y"') = (x, y ',x r tY",xV 2 Y'"). 



Note that $(Z7°) is not compact; Y" and Y'" are 'global' variables. Thus 
is actually continuous on <J>([/°) since r" > 0. Thus, $ is a blow-up and 
is a somewhat singular blow-down map. In the coordinates (x, y' , Y" , Y'") the 
Riemannian density takes the form 

ax""- 1 dx dy = ax-' n+ ^ >f +(»-m)/2-i dx dy i dY " d Y"\ 

a > 0, a G C°°(X). We thus conclude that (for O small) u e 1^(0, M) if 
and only if for any Q € ^~ c °°' (O) with Schwartz kernel supported in U x U, its 
microlocalization Qu satisfies 

(Yy''(Y'y'''(xD x rD^Df ll Df,,Qu 

(5.24) 

e x s+n/2-T. r 3 /2-(„-m)/4 L 2( x -l ^ dj/ / jy///^ 

for every a, /3, 7" and 7"', i.e. if and only if microlocally u is conormal in (x,y') 
with values in Schwartz functions in (Y", Y'"), with the weight given by s + n/2 — 
£rf/2-(rc-m)/4. 
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We also recall that for conormal functions, the I? and the L°° spaces are very 
close, namely they are included in each other with a loss of x e . Thus, u S lie"* (0,M) 
implies that 

(y»y" (Y'"y'" (xD x ) a D^Df„D^'„Qu 

g a , S +n/2-E^72-(n-m)/4-e i oo^-l dxdy , dY " dY'"), 

for every e > 0. 

6. Effectively nonresonant operators 

We now assume that the normal form p n orm for o~\(P(o~)) at q is such that the 
term r er in Theorem 13.111 vanishes. If this is true, we shall call p n orm effectively 
nonresonant, and a an effectively nonresonant energy for q. The significance of 
the notion of effective resonance in general is that the form of the asymptotics of 
microlocally outgoing solutions of Pu = /, / G C°°(X), is independent of r cnr ; only 
r eI changes this form slightly. Moreover, effective non-resonance is a more typical 
condition than non-resonance. We deal with the effectively nonresonant case in 
this section and treat the effectively resonant case in the following section. In both 
cases, it is convenient to reduce P, and not only its principal symbol, to model 
form. This is accomplished in the following lemma. We recall here our ongoing 
assumption (15. 2|) . 

Lemma 6.1. Let p n0 rm be as in Theorem \3.11\ and P as in Remark \ 3.12\ i.e. 
<jQ -i{P) — Pnorm- Then P can be conjugated by a smooth function to the form 

(6.1) 

(771 — 1 n—1 \ 

xD x + r jVj D y + Qj{^ 1/2 y^x 1/2 D y] ) +R el + b + R\ 
3 = 1 j-m, J 

s — 1 m — 1 

where b is a constant, Qj is a real elliptic homogeneous quadratic polynomial (i.e. 
a harmonic oscillator), Vj and Vo are homogeneous polynomials of degree rj resp. 
1, when yk is assigned degree ru, and R £ x t {AA)' J for some j € N and e > 0. In 
addition, for s < j < m — 1, Vj is actually a polynomial in y s , . . . (i.e. is 

independent of yj, . . . , y m -i) without constant or linear terms, while for j < s—1, 
Vj is a polynomial in yj+i, ■ . ■ , y s -i- 

We call Pnorm o, normal form for P. If p n0 rm is effectively non-resonant then 
Per = 0. 

Remark 6.2. Note that Qj(x~ 1 / 2 yj,x 1 / 2 D yj ) is not completely well-defined since 
Qj is a homogeneous quadratic polynomial, and yj and D yj do not commute. 
However, any two choices for the quantization Qj differ by a constant multiple of 
the commutator [x~ 1 / 2 yj,x 1 / 2 D yj ] = [yj,D yj ], hence by a constant. 

In particular, with the notation of the previous section, Qj (Yj , D v. ) may be 
arranged to be self-adjoint with respect to dYj, by symmetrizing if necessary, which 
changes Qj at most by a constant. 
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Proof. With the notation of Lemma [5.7[ any effectively resonant monomial (defined 
in Definition 13. 9|) gives rise to a term of the form x~ x QpB^ with J2k Sk @ k = 1j 
while the effectively non-resonant terms (defined in Definition 13. 10)1 are of the form 
x _1 QpBP with SkPk > 1- This is indeed the key point in categorizing resonant 
terms as effectively resonant or nonresonant; see the proof of Theorem 16.71 But if 
e = s kPk — 1 > 0, we can rewrite x~ x Q^B^ ~ x € QpA@ (i.e. the difference of the 
two sides is in ^~ C °°' (X)), and QpA 13 £ M.'^'. Since there are only finitely many 
effectively non-resonant terms in (|3.21|) . we deduce that any P with <Ji(P) = Pnorm 
may be written 



(6.2) A- X P = xD x +Y J '■,!/, * I, + Qi( x ~ 1/2 yj> xl/2D vi) + R C r + B + R, 



where i? cr is as in (jlTTjl . R £ x e Mj for some e > 0, and B £ ^*f{X). Note 
that Vj and "Po are polynomials, and the homogeneity claim is the meaning of the 
resonance condition Proposition 13.61 For s < j < m — 1, Vj is independent of 
Uj, . . . , Um-i since < r s < r s +i < . . . < r TO _i; itself cannot appear in Vj due 
to the restriction 2a + \j3\ + \ j\ > 3 in Proposition 13.61 Similarly for j < s — 1, Vj 
is independent of yi, . . . ,yj as r± < < ... < r s _i < 0. This also shows that the 
polynomials 'Pj, j 7^ 0, have no constant or linear terms. 

Let B have symbol b(y, y, /1). Modulo terms in x e J^A\ this can be reduced to the 
symbol b'(0, (y', 0, 0), 0). Finally, by conjugating P norm by a function e l ^ y \ we can 
remove the j/'-dependence of b' . Indeed, the Taylor series of / can be constructed 
iteratively. Let T' denote the ideal of functions of y' that vanish at 0. Conjugating 
P by e l f produces the terms X)j=i r 'jV'j^v ' ■) "> as weu as terms from i? cr , which map 
(T) k -> {T) k+1 . For fc > 1, / 1 — ► T,jZlrjVjdy'J defines a linear map on (X') k , 
k > 1, with all eigenvalues negative since 7^ < for j = 1, . . . , s — 1. Thus, this 
map is invertible, and this shows that b' — b'(0) can be conjugated away in Taylor 
series. Then it is straightforward to check that the infinite order vanishing error 
can also be removed. □ 

Later in this section we show that if p n0 rm is effectively non-resonant, the leading 
asymptotics of microlocally outgoing solutions for (|6.ip and for the completely 
explicit operator 



are the same, if R £ x 1+e A43 for some e > 0, i.e. R is indeed an 'error term'. An 
analogous conclusion holds in the effectively resonant case, with R el included in the 
right hand side of (|6.3p . 

First, however, we study the asymptotics of approximate solutions of P$u = 0. 
The constant b simply introduces a power x~ lb into the asymptotics, as can be seen 
by conjugation of Po by x~ th . Here it is convenient to have the asymptotics for the 
ultimately relevant case, where the operator xP is self-adjoint, stated explicitly, so 
we assume that xPq is formally self-adjoint on L^ C (X), which amounts to 



n-1 



m— 1 n—l 




(6.4) 



Im b 



n 



1 



n — m 



2 



2 



3—1 i= s 
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provided that we have already made Qj self-adjoint as stated in Remark 16.21 Note 
that^ = £"=^Rerf. 

For convenience, we separate the case where q is a source/sink of W, hence of 
the contact vector field of Pq. Recall from the previous section that 

(6.5) Y!' = x~ r "y", y'" = x~ 1,2 y"\ 
and define the exponents 

(6.6) b = b- i — - — , ap, = -^rjffj - ib. 

3=1 

Notice that Re api — > oo as — > oo. 

Proposition 6.3. Suppose that the radial point q is a source/sink ofW, and (16. 4p 
holds. Suppose that u £ 7 (s) (0, M), and P u £ I (s '\0, M) where s < -1/2 < s' . 
Then u takes the form 

(6.7) u = J2 x- il - iK *w k {Y")v k {Y'") + u' 

k 

where the sum is over k £ N, v k (Y) is an L 2 -normalized eigenfunction of the 

harmonic oscillator 

(6.8) 

71 — 1 - 

<L>y.rl>Y;- QA) ,.l)y i = QjiYj^y,) - -{Y,D Yj + Dy^j), - J^, 

im'i/l eigenvalue K k , Wfe are Schwartz functions with each seminorm rapidly decreas- 
ing in k, and u' £ I^ s ~ e '(0, Ai) for every e > 0. 

Conversely, given any rapidly decreasing Schwartz sequence, w k , in Y" , meaning 
one for which all seminorm rapidly decreasing in k, and given any f £ Isc (O, 
tfiere exists u £ n s< - 1/2 Is!\o, M) of the form {6/7]) with WF sc (P Q u - f) n O = 0. 

Remark 6.4. The result is true if we only assume s < s'. However, if s > —1/2, we 
can replace s by s > —1/2, apply the proposition with s in place of s, and then use 
u € Is* (O, Ai) to show that each vanishes. On the other hand, if s' > —1/2, 

the proof of the proposition shows that u 6 7 S c (0,M.) implies u £ Is* ~ e ^(0,.M) 
for every e > 0. 

Proposition 6.5. Suppose that q is a saddle point ofW, and (|6.4|) holds. Suppose 
u £ 1^(0, M), and P a u £ I^ s '\0,M) for some s < s' < oo. Then u takes the 
form 

(6.9) u = J2 xaci " lKk (y'f w f3'AY")MY"') + «' 

where the sum is over k £ N and a finite set of multiindices j3' , v k (Y) and K k are 
as above, wp'.k is a rapidly decreasing Schwartz sequence and v! £ I^ s '-^(0,M) 
for every e > 0. 

Conversely, given any rapidly decreasing sequence of Schwartz functions wp' jk , 

finite in [3' and any f £ Isc \0,Ai) there exists u £ n s< _ 1 / 2 /ic' > (O, M) of the 
form (EU) with WF SC (P « — /) fl O = 0. 
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Remark 6.6. As shown later, x 2 D x gives rise to the terms in Q— Q after the change 
of variables (x,yj) i— > (x, -r^r). If Qj is self-adjoint on L 2 (R,dYj) then Q., has the 
same property. 
Also, with 

n- 1 



2 2^ J 4 ' 

i 

the summand in (ESJ is in 7 s ( f ca ' 3 ' _ ' B_ ^~ e) (0, M) for every e > 0. We 

show below that Im& = B + d, d = -ij]r^ > 0, so the (/?',&) summand is in 

7 BC J j 2 (0,jV1) for every e > 0, and in view of the rapid decay in k, the 

same is true after the k summation. Thus, for u as in (|6.9p . it € lie 2 e ^(0,.M) 
provided s' > d — i.e. decays by a factor x d faster than the microlocal solutions 
at sources/sinks of W. 

Proof of Proposition \6.3[ Suppose that P u = f e I^ s \0,M) for some s' > -1/2. 
Let O 1 be a PF-balanced neighbourhood of q with O 7 c O, and let Q S * S ^. OC ' (A) 
satisfy WF^Q) C O (i.e. Q £ *- DO »°(0)) and WF^Id-Q) n O 7 = 0, with 
Schwartz kernel supported in U x U , 

U = {0 < x < 5, \yj\<6 for all j}. 

(See (|5.23l) for the definition of the diffeomorphism the coordinates Yj, etc.) 
Then, as noted in (|5.24[) . by the definition of 7 S c (0,M), u = Qu satisfies 

(Y"y" (Y"'y'" (xD x ) a Df„Df,'„u £ x s L 2 sc (X) 

for all a, /3", f3"', 7" and 7"'. Here u is a microlocalization of u since WF SC («- Qu) C 
WFs C (Id-Q), so WF sc (ti - Qu) HO' = 9. Moreover, 

P (Qu) = QP Q u + [P , Q]u = Qf+ /', /' G C°°(X), 

since WF sc (m) HO C {<?}, while WF' SC {[P Q , Q}) C WF' SC {Q) n WF^ c (Id -Q) C 0\D T , 
so WF sc (u) n WFg C ([P , Q]) = 0. Thus, with f = Qf + f, 

1 = ^ 

(Y'r"(Y"r"'(xD x ) a Df„Df„f e x s 'L 2 JX), 

for all a, /?", 7" and 7'". 

To prove first part of the proposition, it thus suffices to show that, with the 
notation of (16.71). 



(6.11) u = J2 x~ tl, ~ lKk w k (Y")v k (Y'") + u'. 

k 

Writing the operator P in the coordinates x, Y", Y'" we have 

(6.12) P Q = xD x \ Y + QiW'i d y>>>) + b 

j 

with b = b — i n ~ 4 rn as in (|6.6p . Formal self-adjointness of xPq, i.e. (|6.4p . requires 
that 

T r n—l 1 \ ,, n — m „ 

(6.13) Imb = — - -XX' - ee fl. 
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As already remarked, (I6.10| . which states that / is conormal in x, and Schwartz 
in Y", Y'", and belongs to x s L 2 (dxdy / x n+1 ) , or in terms of the Y coordinates, to 
x s'+n/2-T,r<'/2-(n- m )/4 tf( dxd Y/ x )^ implies (by conormality) that 

f g a , S ' + l/2+B-e L oo 

for every e > 0, where B is defined by (16.13p . More precisely, for all a, ft, 7" and 



(Y")" 1 " '(Y"Y" '(xD x ) a Df„Df',J e a; -'+ 1 /2+s-e i oo 

for every e > 0. Conversely these conditions imply that / satisfies (I6.10P with s' 
replaced by s' — e for every e > 0. 
Writing / in the form 

f(x,Y",Y"') = J2Mx,Y")v k (Y">), 

k 

where is conormal in x, rapidly decreasing as a Schwartz sequence in Y" , a 
particular solution to Pqu = f, is given by 

u = Y,Mx,Y")v k (Y>"), 

(6 - 14) k r - dt 

u k = -i x - ib - iK " / f k (t,Y")? b+iK * 



t 

Since s' + 1/2 > 0, this integral is convergent and u £ I^ s ~ £ - 1 (O, M.) for every e > 0. 

On the other hand, the general solution to Pqu = with u Schwartz in Y" and 
Y'" is given by 

(6.15) u = x- ll - lK »w k {Y")v k (Y"'), 

k 

where w k is rapidly decreasing in k. Since any solution is the sum of the particular 
solution (|6.14p and some homogeneous solution, the first half of the proposition 
follows. 

In fact, the second half also follows by defining 

u = Y J Mxy')MY m )+Y J x ~ li ~ lKk MY>k{Y"') 1 

k k 
with u k as in (|6 . 14[) . Multiplying by a cutoff function (f> £ C°°(X) which is identi- 
cally 1 near (0, 0, ... , 0), it follows that u = <fiu satisfies all requirements. □ 

Proof of Provosition 1 6. 51 We use a similar argument to prove this result. Let O' , 
Q, etc., be as in the previous proof. With u = Qu, as noted in (|5.24p . 

(6.16) (Yy"(Y / y'"(xD x ) a D^Df„Df'„ u e x s L 2 sc (X), 

for all a, ft, 7" and 7"'. One of the main differences with the proof of Proposition l6.3l 
is that microlocalization introduces a non-trivial error, i.e. Pqu is not globally well- 
behaved (not as good as / was microlocally) . However, the error is supported away 
from y' = 0. Indeed, now WF SC («) D O C S, and 

f = P u = Qf + f, f = [P ,Q]u. 

Here WF^ C ([P ,Q]) fl5c {\y'\ > 6 Q } for some S > 0, so /' € I$(0,M) in fact 
satisfies 

(Y"r"(Y"'r"'(xD x rD^Df„Df„f' e x s L 2 sc (X) 
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for all a, j3' ', f3" and (3"', 7" and 7"', with the improved conclusion 
( f>(y')(Yy\Y'y'''(xD x ) a D^D^,Df, 1 f' € C°°(X) 
if <j> is supported in \y'\ < So- Correspondingly, 

(6.17) 4>(y')(Yy\Y'y'\xD x YD^D^,Df,,~f € x s ' Ll(X). 
The operator p) in the coordinates x,y', Y", Y'" now takes the form 

(6.18) P = xD x \y>, Y r> t Y'» + J2 '"'"''V + Qi( Y J"> D y?') + 6. 

with b = b — i "~ m as in (|6.6|) . Again, (|6.17p implies that / is conormal in x, 
smooth in j/', and Schwartz in Y",Y"', and belongs to x s + 1 / 2 + B ~ e j J °° fo r every 
e > 0, where B is defined by (|6.13[) . in the precise sense that for all a, /3, 7" and 
7'", 

4>{y')(Y'y" {Y"y"' {xD^D^DyZj G x s '+ 1 /2+s-e L oo 
for every e > 0. However, now formal self-adjointness of xp) requires that 

(6.19) Im& = B + d, d = -^J2 r 'j >0 > 

3 

so there is a discrepancy of d compared with the previous proposition. Write / in 
the form 

f{x,Y\Y">) = Y J fk{x,y',Y")v k {Y">), 

k 

where f k is rapidly decreasing sequence which is conormal in x, smooth in y' and 
Schwartz in Y" . 

We start by describing solutions of the homogeneous equation Pqu = in U 
which in addition satisfy (|6.16p . Decomposing u in terms of the v k , and factoring 
out a power of x for convenience, i.e. writing u = J2 k x~ lb ~ lK,k u k (x, y', Y")vk(Y"'), 
we see that the coefficients u k satisfy 

(xd x \ y > t Y",Y»> + ^2 r 'jy'j9y' j )u k = 0. 

j 

Since u is smooth in the interior of U, Pqu = amounts to demanding that u k be 
constant along each integral curve segment of the vector field xd x + ^2j r 'jy'jd y '., 
with the value of u depending smoothly on the choice of the integral curve. (We 
remark that U is convex for this vector field; \y'\ is increasing as x — ► 0.) Thus, 
u k (x,y',Y") = u k {Y',Y") with u k smooth in Y' and Schwartz in Y" . Here Y- = 

y'j/x r i ; note that t'a < 0. Expanding u k in Taylor series around Y' = to order N, 
we see that 

u k (x,y',Y")- x~^ r '^(y'fwp'AY") 

\0'\<N-1 

is a finite sum of terms of the form x~ r i®i (y 1 ) 13 ' u k ^i (Y', Y") with u kt p> smooth 
(Schwartz in Y'"), where the sum runs over f}' with \j3'\ = N. Thus, given any 
s" (e.g. a" = s'), we can choose N sufficiently large so this difference lies in 
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Isc \o,M), which means it is ignorable for our purposes. Thus, the general solu- 
tion to Pqu = in U which satisfies (|6.16|) is given by 

(6.20) u = J2 xa *'~ lKk (y'f^>AY")vk(Y'"), 

modulo any Jj* \0,M) (where the sum is understood as a finite one, due to the 
remark above), where the seminorms of wpi k are rapidly decreasing in k for each 
&• 

In expressing a particular solution u of Poii — f in terms of /, we need to 
integrate along integral curves of the vector field xd x + J2j r 'jV'j^v' < an d since r'^ < 0, 
\y'\ ► oo as x — * along such curves (unless y' = 0); in fact \y'\ is increasing as 
x — > as mentioned above. So we cannot integrate down to x = 0. Instead we fix 
an xq > and use the formula 

u k (x,y',Y") = (^-)-^-^ Uk (x^{^-)-^y' p Y") 



(6.21) 



Xq x 
+ix- il - i «* f X f k (t,(j)^y' v Y")t fb+lKk T- 

Jxo 1 1 



Notice that u k (x^, yL, Ya') depends only on f k evaluated at points [x, y',Y") with 
W\ < \yi\- Thus, (|6.17p can be used to deduce properties of u k , hence of u, in 

\y'\ < s . 

If s' < -1/2 + d, then (jOT|) gives <j>{y')u e 1^-^(0, M) for every e > 0, 
with <j) as in (JHHTJ). If s' > -1/2 + d, then <j>{y')u E I^ 1/2+d ^^ (O, M) for every 
e > 0. However, this is actually a sum of terms solving the homogeneous equa- 
tion, plus a function in ~ £ )(0,.M) for every e > 0. For simplicity we show 
this only in the case that —1/2 + d < s' < —1/2 + d + 1 7*^ _ 1 1 - Then we observe 
that (x/x )- lb - iK *ti(x ,Q,Y") is a solution of the homogeneous equation, while the 
difference 

{ ±yib-i Kh ~^ (JL)-r' jy > Y ") - (-)-^u(xo, 0, Y") 
Xo Xo J Xq 

= E(£)"^ / ; »:<>, (u(xo,r(^r r h' r Y")) dr 

■ Xq Jo Xo 

has decay at least x~ Ts ~ 1 better, hence yields a term in ~ e \0,A4) for every 
e > 0. Similarly, if we replace f k (t, (j)~ rj y'j,Y") in the integral by f k (t, 0, Y") then 
we get a homogeneous term, while the difference gives a term in /( S '- £ )(0,.M) for 
every e > 0. The argument can be repeated, removing more and more terms in the 
Taylor series for u and /, for larger values of s' . Since any solution is the sum of the 
particular solution above and the general solution, the first half of the proposition 
follows with O replaced by a smaller neighbourhood O" of q. However, we recover 
the original statement by using the real principal type parametrix construction of 
Duistermaat and Hormander [3] . 

The second half can be proved as in the previous proposition. Fix some xo > 0, 
and let u k be given by the second term on the right hand side of (I6.2ip . and let 
u = J2 k u k(x, Y")v k {Y"'). Then Poii = /, and as shown above, u has the form 

(6.22) u = ^'-^ (y'f^, k (Y")v k (Y"') + u', 

f3',k 
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with u' G 7 s ( c' e) {0,M) for all e > 0. Then with 

u = Y,u k {x,Y")v k {Y'") +J2x^'- iK "(w , }k (Y") -w p ,, k (Y"))v k (Y'"), 

k f3',k 

u — <fiu, <j> G C°°(X) identically 1 near (0, . . . , 0), u satisfies all requirements. □ 

These results on the explicit normal form Po then allow us to parameterize 
microlocally outgoing solutions for every effectively nonresonant critical point. 

Theorem 6.7. Suppose that P(o~) is effectively nonresonant at q, with normal form 
Pnorm near q as in Lemma \6.1[ and (|6.4[) holds. 

(i) If in addition q is a source/sink of W , then any microlocally outgoing 
solution u of P n0 rm has the form (|6.7p . and conversely given any Schwartz 
sequence of Schwartz functions w k there is a microlocally outgoing solution 
u of P n0 rm which has the form (|6.7[) . Thus, microlocal solutions at a 
source/sink ofW are parameterized by Schwartz functions of the variables 
(Y",Y"'). 

(ii) If q is a saddle point of W , then all microlocally outgoing solutions are 
in x~ 1 / 2+€ L 2 for some e > 0. For each monomial (y')P in the variables 
y' , each igN and each Schwartz function w(Y") there is a microlocally 
outgoing solution of the form 

(6.23) u = J2x a "'~ lKk (y'f w(Y")v k (Y'") +u', 

k 

where v! is in a strictly smaller weighted L 2 space than u, and every mi- 
crolocally outgoing solution is a sum of such solutions, with the w = w k> p> 
rapidly decreasing as k — ► oo in every seminorm. 

Proof. First, P nor m = A(P + R), Re x e M j , e > 0. Thus, if O is a neighbourhood 
of q as above, WF 8C (P norm u) n O = 0, then u S 1^(0, M) for all s < -1/2, so 
Ru G lt' ] (0,M) for some s' > 1/2. Hence P u = A _1 P norm w — Ru G l£\0,M). 

If q is a source/sink of W, then Proposition ^. 31 is applicable, and we deduce that 
u is microlocally of the form (|6.7|) . Moreover, if q is a source/sink of W, then given 
any Schwartz sequence of Schwartz functions w k , let uq € C^ s <~i/2lsc\o, A4) be 
of the form (|6J]) with P Q u G C°°(X). We construct u k G n r< _ 1/2 - ke Isc } (O, M), 
k > 1, inductively so that Pou k + Ru k -i G C°°(X) for k > 1; this can be done 
by the second half of Proposition 16.31 Asymptotically summing J2 k u k to some 
u G n s< _i/ 2 /ic'(0,7W) gives a microlocally outgoing solution with the prescribed 
asymptotics, completing the proof of the theorem in this case. 

If q is a saddle point of W, we apply Proposition 16.51 with s' > —1/2 as in 
the first paragraph of the proof. If e' > is sufficiently small, all of the terms 
in (|6.9[) are in /i c 1 / 2+e ^(0,A4) proving the first claim. To show the next, let 
u Q = x a ^'- iKk (y')^' w(Y")v k (Y'"), so P u = and u a G I${0,M) for any s < 
— 1/2 + d. We construct inductively as above, using Proposition 16. 5[ to obtain 
u. □ 

Remark 6.8. From (|6.7j) or (|6.23|) it is not hard to derive the asymptotic expan- 
sion of eigenfunctions of the original operator A + V — a\ we need only apply the 
Fourier integral operator F^ 1 arising by composing any Fourier integral operators 



48 



ANDREW HASSELL, RICHARD MELROSE, AND ANDRAS VASY 



with canonical relation given by the contact maps in Lemma 12.71 and Theorem l3.11l 
to these expansions. In fact, as mentioned in Remark I3.14( this Fourier integral 
operator can be taken to be a composition of a change of coordinates with multi- 
plication by an oscillatory function if q is either a source/sink (so q € Min + (<7)) or 
the linearization of W has no non-real eigenvalues (so there are no y'" variables). 

In the case of a radial point q G Min + (cr), in appropriate coordinates y on dX, 
the expansion takes the form 
(6.24) 

u = e^/ x ^2x- ll ~ tKk w k {Y")v k {Y"') +u\ u' S r^ +t {0,M) for some e > 

k 

where $ is a smooth function (it parameterizes the Legendrian submanifold which 
is the image of the zero section under the canonical relation of F^ 1 ). For a given 
a, only a finite number of terms in the Taylor series for $ are relevant. Similarly 
in the case of radial points q € RP+(er) \ Min+(er), the expansion (|6.25p takes the 
form 

(6.25) u = e**^)/* x a P'~ lKk (y'fw(Y")v k (Y"') + u', 

k 

with $ smooth. Again it parameterizes the image of the zero section under the 
canonical relation of F~ x . In this case, the value of $ on the unstable manifold 
{y" = y'" = 0} is essential, but only a finite number of terms in the Taylor series 
for $ about this unstable manifold are relevant. 

These expansions were obtained directly in Part I (i.e. without going via a normal 
form) in the two dimensional case. 



7. Effectively resonant operators 

If P is effectively resonant, the simple expressions (|6.7[) and (|6.9p need to be 
replaced by slightly more complicated ones in which positive integral powers of 
log a; also appear. Essentially, instead of powers, or Schwartz functions, of -~rj, 
factors of log x also arise in the expressions for the Y\ . 

First define a change of coordinates inductively that simplifies the vector field 

m— 1 

(7-1) V = (xD x ) + J2 (rjVj + VjiVs, • • • . Vj-i))D Vi 

3=1 

that appears in (|6.ip as the combinations of the linear terms ^TjyjDy. and the 
effectively resonant vector fields in i? or . (Note that rjyj and Vj(y s , ■ ■ ■ iUj-i) are 
both homogeneous of degree rj.) We do this in two steps to clarify the argument, 
first only dealing with the y" terms, i.e. j = s, . . . , m — 1. 

The coordinates Yi, j — s. . . . , m — 1, are a modification of the coordinates 
that appear in (|6.5p . so that Yj — -^j are polynomials pj in Y s , . . . , Yj_x,t = \ogx. 
Thus, we let 

Y s =^r,Vl = 0, Y s (Y s ,\ogx) = Y s +V^ogx) 

T ! s 
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and provided that Y s , . . . , Yj-i, V\, . . . ,P\-\ have been defined, we let 

p]{Y s , Yj_ u t) = f Pj(Y s (Y s , t'), . . . , Yj-^Yg, Yj.ut')) dt', 
Jo 

Y j = ^--V i j (Y s ,...,Y j . 1 ,logx), 

Y 3 = y j + ^O^. Yj-ulogx), j = s,...,m-l. 

The point of the construction is that V annihilates Yj for all j. This can be seen 
iteratively: for Y s this is straightforward, and if VY S = . . . = VYj_i = then (with 
d t Vj denoting the derivative with respect to the last variable, t = logx) 

VY 3 

= ~ T 3~ + (r 3 y 3 +V,(ys, . . .,yj-i))x- r > - (d t p]){Y s , . . . .^-i.logaO 

= Vj{y.x- T -,. . . , yj -ix- r ^) -VjiZiYsAogx), . . .,Yj-x{Y„ . . . , Yj-t, logx)) 
= 

in view of the definition of Y s , . . . , Yj_i and Y s , . . . , Yj_\. 

One can deal with the j = 1 , . . . , s — 1 terms similarly. We define Pj , Yj and Yj 
inductively as above, starting with Y s -\. Thus, we let 

y -i = ^T' V s-i =0, Y s _ 1 (Y s _ 1 ,logx)=Y s _ 1 +P t s _ 1 (logx) 

and provided that Yj + \ , . . . , F s _i, Vj +1 , . . . , P\_\ have been defined, we let 
(7.2) 

p]{Y j+u Y s _ u t) = f Pj(Y j+1 (Y j+1 , . . . , Y s _ u t'), Y s -i(Y s - U t'j) dt> , 
Jo 

Y 3 = §--'P i 3(y3+i,---,ys-^ogx), 

Yj = Yj + ^(Yj+u ■ • • , ^-1, logx), ./ I a 1. 

With these definitions, in the coordinates X = x, Y\, . . . , Y m _i,y m , . . . , y n -i, i.e. 
(X,Y',Y",y"'), which correspond to a blow-up of x = y s = . . . = y m -i = 0, 
V = X 2 D X . 

The zeroth order term is a polynomial Pq in y s , . . . , y m -\ which is homogeneous 
of degree 1 (where yj has degree Tj). Thus, 

x^Poiys, ■ ■ -,y m -i) = P (Y s (Y s ,iogx), . . . ,y m _i(y s , . . . ,y m _i,iogx)). 

Let 

p*{Y a , . . . , Yj_ u t) = [ p (Y s (Y s ,t'), . . . , y,-i(y , . . . , y,-i, f )) < 

Jo 

which is thus a polynomial in y, . . . , Y,_i, i. Then e ,p o( y '— ^j-^ 1 ^ 1 ) can be used 
as an integrating factor, conjugating P, to remove the zeroth order term in R CI . 
Finally, to put the quadratic terms in a convenient form, we let 

Yj = J^, j = m,...,n-l 

as before. 



50 



ANDREW HASSELL, RICHARD MELROSE, AND ANDRAS VASY 



Suppose first that Vo = 0. With our definition of the Yj, (I6TT2]) . resp. ([67T8]) . 
holds if q is a source/sink, resp. saddle point, of Vq. Thus, the statement and the 
proof of Proposition 16.31 holds without any changes, while the statement and the 
proof of Proposition 16 . 51 carry over provided x a C {y'Y is replaced by x~ lb (Y')P . A 
minor difference is that slightly more effort is required to show that \y'\ decreases on 
the integral curves of the vector field (|7.1|) inside \y'\ < 8\ for 8\ > small. Namely 
we need to use that, as Vj, j = 1, . . . , s — 1 have no linear or constant terms by 
LemmaEU V\yf = r ]V ) + 0(\y'\ 3 ) < r s ^\y'\ 2 + 0(\y'\ 3 ), r 5 _i < 0, to 

conclude that V\y'\ 2 < for \y'\ < Si, <Si > small. 

In general, with 6 = 6 — z 2 -^ 1 as in (|6.6|) . (|6. 12[) . resp. (|6. 18[) . are replaced by 

(7.3) P = xD x \ Y + Qj (Xj"> ° Y r ) + V o + b, 

j 

respectively 

s-1 

(7.4) P = xD x \ y ,, Y „, Y ,„ + Yjir'iV'i + Vj)D y , + ]T Qj(Y-", Dyy) +V a + b. 
Thus, 

(7.5) e^'Poe-^o = xDg \ Y + £ Qj(Y"' ,Dy><>) + b, 

j 

respectively 

s-1 

(7.6) e^oPoe-^* = xD x \ v ,y„ y „, + £(rj^ + P,)P y , + ]T Q,(y/", D Y >>> ) + b. 

3=1 3 

Since multiplication by e ±l7: o preserves Isc\o, A4), the rest of the proof of the 
propositions is applicable with u replaced by e "»«, / = p)U replaced by e °f. 
We thus deduce the following analogues of Propositions 16.31 - 16.51 in the effectively 
resonant case. 

Proposition 7.1. Suppose that the radial point q is a source/sink ofW, and (|6.4p 
holds, that u G 1^(0, M), and P u S I (s,) (0, .M) w/iere s < -1/2 < s'. T/ien u 
takes the form 

(7.7) u = J2x~ il - iK "e- i ' p ow k (Y")v k (Y'")+u' 

k 

where the sum is over k G N, v k {Y) is an L? -normalized eigenfunction of the 

harmonic oscillator 

(7.8) 



n-l 



tv. n_. i a. v.\ v. 



X Q.O rlh ,. QA) ,.l)y i = Q^Dy ; - + D Yj Yj), Yj 



with eigenvalue K k , w k are Schwartz functions with each seminorm rapidly decreas- 
ing in k, and u' £ ~ e \0,A4) for every e > 0. 

Conversely, given any sequence Wk of Schwartz functions in Y" with each semi- 
norm rapidly decreasing in k, and given any f £ I^'(0,Ad), there exists u £ 
n s< _ 1/2 I^ ] (0,M) of the form ([72]) with WF sc (P u - /) n O = 0. 
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Proposition 7.2. Suppose that q is a saddle point of W , and (|6.4| holds, that 
u S 1^(0, M), and P u £ \0,M) for some s < s' < 00. Then u takes the 
form 

(7.9) u = J2 x^-^ {Y'fe- iV "wp., k (Y")v k {Y'") + u' 

where the sum is over k £ N and a finite set of multiindices f3' , v k (Y) and n k are 
as above, Wpi jk are Schwartz functions with each seminorm rapidly decreasing in k, 
and u' £ I<> s '~^{0,M) for every e > 0. 

Conversely, given any sequence of Schwartz functions Wp^kt finite in (3' with 
each seminorm rapidly decreasing in k, and any f £ J S c \Q,M) there exists u £ 
C\ s< - 1/2 I { S c{0,M) of the form (0 with WF sc (P Q u - /) n O = 0. 

We thus deduce the following analogue of Theorem 16. 7\ with a similar proof. 

Theorem 7.3. Suppose that P(o~) is effectively resonant at q, with normal form 
-Pnorm near q as in Lemma \6.1[ and (|6.4|) holds. 

(i) If in addition q is a source/sink of W , then any microlocal solution u of 
Pnorm has the form (|7.7j) . and conversely given any rapidly Schwartz se- 
quence of functions w k there is a microlocally outgoing solution u of P noTm 
which has the form (I7.7[) . Thus, microlocal eigenf unctions at a source/sink 
are parameterized by Schwartz functions of the variables (Y",Y'"). 

(ii) If q is a saddle point ofW, then all microlocal solutions are in x~ 1 / 2+e L 2 
for some e > 0. For each monomial in the variables Y' , each k £ N and 
each Schwartz function w(Y") there is a microlocally outgoing solution of 
the form 

(7.10) u = x- fb -' lKk e-<(Y'fw(Y")v k (Y'")+u', 

where u' is in a strictly faster decaying weighted I? space than u, and 
every microlocally outgoing solution is a sum of such solutions, with the 
w = w k .f3' rapidly decreasing as k — > 00 in every seminorm. 

8. From microlocal to approximate eigenfunctions 

We are interested in the structure of (global) eigenfunctions of A + V. While in 
the first half of the paper a rather general element P £ 1> S*^ 1 (X) was considered, 
from now on attention is limited to 

H = A + V £ ^* s f{X), H(a) =H-a, 

in particular the order of H at dX is 0. 

In the next section we obtain an iterative description of the 'smooth' eigen- 
functions in terms of the microlocal eigenspaces. As the first step, we show that 
if q is a radial point for H(o~) = H — a, then elements of _E m i Ci+ (g, a), which 
are the microlocally outgoing eigenfunctions near q, have representatives satisfying 
(H — a)u £ C°°(X), i.e. they extend to approximate eigenfunctions, with WF sc (w) 
a subset of the forward flow-out of q. Stated explicitly this is 

Proposition 8.1. If q £ RP + (ct) then every element of E m ic,+ (q, f) has a repre- 
sentative u such that (H — a)u £ C°°(X) ) and WF sc (u) C 
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Remark 8.2. From this result, given u as in Proposition 18. II it is easy to produce 
an exact eigenfunction v such that WF sc (w) n {v > 0} C $+({(?}): we simply take 
v = u — R(a — iO)(H — a)u. 

The key ingredient of the proof, as in the two-dimensional case studied in [5], 
is the microlocal solvability of the eigenequation through radial points. To avoid 
a microlocal construction along the lines of Hormander [13] . we introduce, as in 
Lemma 5.3], an operator H which arises from H by altering V appropriately. 
This is chosen to be equal to H near the radial point in question but to have no 
other radial points in RP+ (cr) at which v takes a smaller value. One may then 
assume, in any argument concerning q G RP+(<r), that there is no q' G RP+ (cr) 
with v(q') < v(q). 

As in [SJ Definition 11.3], we introduce a partial order on RP + (cr) corresponding 
to the flow-out under W. 

Definition 8.3. If q, q' G RP+ (cr) we say that q < q' if q' G $+({q}) and q < q' 
if q < q 1 but q 1 + 1 q. A subset T C RP + (cr) is closed under < if, for all q G T, 
{q' G RP+ (cr);q < q'} C T. We call the set {q 1 G KP + (a);q < q'} the string 
generated by q. 

Remark 8.4. This partial order relation between two radial points in RP + (er) cor- 
responds to the existence of a sequence qj G RP+(ct), j = 0, . . . , k, k > 1, with 
9o = 9, 9fe = q' and such that for every j — 0, . . . ,k — 1, there is a bicharacteristic 
7j with limt^oo 7^ = ^ and lim f ^ +00 7^ = t^- +1 . 

Lemma 8.5. Given a > minVb and v > 0, set K — Vo _1 ((— 00, a — v 2 ]) C 9A 
then there exists a potential function V G C°°(X) with Vq Morse such that 

(i) V > V , 

(ii) Vq = Vq on a neighbourhood of K , 

(iii) no critical value of V lies in the interval (cr — v 2 ,a\, 

(iv) if Yi{a) is the characteristic variety at energy a of H — A + V then 

(8.1) S(o-) n {v>v} = E(o-) n {iv > £>} and 

(v) H — o has no L 2 null space. 

Proof. Choose a smooth function / on the real line so that /' > 0, /(t) = t if 
t < cr-£> 2 and /(t) > a for t > min{V(g); dV(g) = and V(g) > a - v 2 } > a - v 2 . 
Then let V = / ° V, so the critical points of Vo and Vo are the same and are 
non-degenerate. 

On S(cr) n > £>}, i/ 2 + |^i| 2 + V Q = cr, hence V < cr - i> 2 , so V = V , and 
therefore £(er) fl > C S(cr). With the converse direction proved similarly, (JTJ) 
— (jrvj follow. Property (|vj can be arranged by a suitable perturbation of V with 
compact support in the interior. □ 

These properties of H are exploited in the proof of the following continuation 
result. 

Lemma 8.6 (Lemma 5.5 of [5]). Suppose u G C~°°(X) satisfies 

WF sc (u) C {v > ^1} and WF SC {{H - a)u) C {v > v 2 }, 

for some < v\ < 1/2, then there exists u G C~°°{X) with WF sc (ti — u) C \v > ^2} 
and {H-a)u G C°°(A). 
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Proof. We just sketch the proof here; for full details, see [5]. The obvious idea of 
subtracting R(o~ + iO)((H — cr)u) from u does not quite work, since the forward 
fiowout of other critical points q' G RP+(<r) with v{q') less than v(q) may strike 
q. To avoid this problem, choose v with V\ < v < sufficiently close to i>2 
so that there are no radial points q with v(q) G [v, 1/2), and a corresponding V 
as in Lemma [8.51 Then consider the function R(a + iO)(H — a)Au, where A is 
equal to the identity microlocally on {v < v} PI S(<x) and vanishes microlocally in 
{y > V2\. Since Vb has no critical points q with < v(q) < V2 it follows readily 
u = Au — R{a + iO)(H — a)Au satisfies the desired conditions. □ 

From this we can readily deduce 

Lemma 8.7. // q G RP+ (a) then every element of E m i Ci +(q, a) has a representative 
u such that (H — a)u G C° C (X) and WF sc (u) is contained in the union of <&+({(?}) 
and the $ + ({q'}) for those q' G RP + (er) with v{q') > v{q). 

Proof. Let O be a W^-balanced neighbourhood of q (see Definition 14. 4|) . Let A G 
^~°° ,0 (X) be microlocally equal to the identity on $ + ({g}) (~l O and supported 
in a small neighbourhood of <J>+({g}) n O. Then there exists V2 > v{q) such that 
v > V2 on$+({q})\0, and WF^ c (A)\0 C {v > v 2 }. (HereWF^ c (A) is the operator 
wavefront set of A, i.e. the complement in sc Tg X X of the set where A is microlocally 
trivial; see [TS].) Now let u be any representative. Since WF sc (u) nOc $+({q}), 
WF sc (Au - u) n O = 0. In addition, WF sc (Aw) c WF^ C (A) n WF SC («), hence 
v > v(q) on WF sc (Au). Moreover, WF sc (Au - u) n O = implies that 

WF sc ((ff - a)Au) flO = WF sc ((iJ - a)Au - (H - a)u) n O = 0, 

so WF sc ((iJ - cr)Au) C WFg C (A) \ O, hence is contained in {v > v 2 }. Then, by 
Lemma [8.61 there exists u G C~°°(X) such that ^ > ^2 on WF sc (u — Au) and 
(H — a)u G C oc (X). In particular, v > v(q) in WF sc (u). Moreover, v > V2 on 
WF sc (u — u) n O, hence by Lemma r4~5l WF sc (u — u) n O = 0, so u and u have the 
same image in E mlc ^ + (0, a). □ 

Finally, we can show that each microlocally outgoing eigenfunction is represented 
by an approximate eigenfunction. 

Proof of Proposition lff.il Let u be a representative as in Lemma 18.71 If we choose 
q' from the set 

(8.2) {q' G RP+(a) n WF SC (5); v{q') > v{q), q' $ $+({g})}, 

with v{q') minimal, then, localizing u near q', gives an element v of E m i Cy+ (q'). 
By subtracting from u a representative of v given by Lemma 18.71 we remove the 
wavefront set near q' . Inductively choosing radial points from (|8.2[) and performing 
this procedure repeatedly, all wavefront set may be removed from u except that 
contained in & + ({q}). □ 

9. Microlocal Morse decomposition 

Next we show that global smooth eigenfunctions can, in an appropriate sense, 
be decomposed into components originating, in the sense of the Introduction, at 
a single radial point. We do this by defining subspaces of E^ s (a) corresponding 
to the location of scattering wavefront set in {v > 0} and showing that suitable 
quotients of these spaces are isomorphic to the spaces of microlocal eigenfunctions 
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E™ ic + (q, a), q € RP+(cr), analyzed in sections [6] and [71 Since each of the spaces 
-^mic + (5,^), 9 £ RP+(cr), is non-trivial this shows that each such radial point gives 
rise to eigenfunctions. However, as noted previously in [9], [10j . and [5] in some 
special cases, there is a qualitative difference between the radial points correspond- 
ing to local minima of Vo and the others. This is expressed by Proposition 110.31 
where we show that the eigenfunctions u € E^ in + (<r) originating only at minimum 
radial points are dense in E® ss (a) (definitions of these spaces are given below). 
Recall from [5l Equation (3.14)] the spaces of eigenfunctions of fixed growth 

(9.1) £| s » = {ue E~^(a);WF° s r 1/2 (u) n {v = 0} = 0}. 
This condition is equivalent to requiring that 

(9.2) Bu€x s -^ 2 L 2 (X) 

for some pseudodifferential operator B 6 fy®]9(X) with boundary symbol which is 
elliptic on E(er) n \v = 0} and microsupported in < a(er)}, where 

a(cr) = mm{\v(q)\;q G RP(tr)}. 

The space E^ ss (ct) is of particular interest. Choose an operator A 6 ^^"(X) whose 
boundary symbol is for v < —a(a) and 1 for v > a(a). The space E® ss (a) is a 
Hilbert space with norm 

(9.3) \\u\\ 2 E o Be{(7) = (i[H,A]u,u). 

The positive-definiteness of this form, and its independence of the choice of operator 
A, was shown in [5], Section 12. An equivalent norm is 

(9.4) \\Bu\\ x -i/2 L 2 + 11^113.-1/2-6^2 

where e > and B is as in (|9.2|) : see [5], section 3. 

We now define subspaces of E^ ss (a) depending on the location of the scattering 
wavefront set inside {v > 0}. Given any <-closed subset T of RP+(er), we define 

(9.5) E s css (a,T) ={ue E s css (a);WF sc (u) n KP+(a) C V}. 

The set of radial points q E RP+ (a) lying above local minima of V is an example of 
a <-closed subspace and will be denoted Min + (cr). In this case we use the notation 

^ in>+ (a) ee E s eBB (a,Min+(a)) = {u E E s css (a);WF sc (u) n RP+(a) C Min + (<r)} 

to be consistent with [5]. 

Proposition 9.1. Suppose that T C RP+(cr) is <-closed and q is a <-minimal 
element of T. Then with V — T \ {q} 

E™(a, V) -U- E™(a, T) — X- E mic , + (a, q) 

is a short exact sequence, where 1 is the inclusion map and r q is the microlocal 
restriction map. 

Proof. The injectivity of l follows from the definitions. The null space of the mi- 
crolocal restriction map r q , which can be viewed as restriction to a VF-balanced 
neighbourhood of q, is precisely the subset of E^ s (a,T) with wave front set dis- 
joint from {q}, and this subset is E^ s (a,T r ). Thus it only remains to check the 
surjectivity of r q . 

We do so first for the strings generated by q € RP+ (er). For q 6 Min + (a), 
the string just consists of q itself and the result follows trivially. So consider the 
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string S(q) generated by q E RP+ (er) \ Min+ (cr) . By Proposition 18 . II any element of 
-E'mic,+ (q , a) has a representative u satisfying (H — a)u £ C°°(X) with WF sc ( , i2) C 
$+({(?}). Then u = u — R(a — iO)(H — a)u E E^ S ((T, T), which gives surjectivity in 
this case. 

For any <-closed set T and <-minimal element q, the string S(q) is contained in 
r, so the surjectivity of r q follows in general. □ 

Notice that we can always find a sequence = Tq c T± C . . . C T n = RP + (cr), 
of <-closed sets with Tj \ Tj-i consisting of a single point qj which is <-minimal 
in Tj\ we simply order the qi E RP+(cr) so that v{q\) > ^(92) > and set 
Ft = {qi, ■ ■ ■ , Then Proposition 19. II implies the following 

Theorem 9.2 (Microlocal Morse Decomposition). Suppose that = To C Ti c 
. . . C r„ = RP + (er), is as described in the previous paragraph. Then 

(9. 6) {0 } _ EZ (a,^) EZ(a,T n ^) - 

(9.7) E~ (a, T^/E™ (cr, r 3 -i) ^ ^ m ic,+ fe, tr), J = 1, 2, ... , n. 

10. L 2 -PARAMETERIZATION OF THE GENERALIZED EIGENSPACES 

Recall from Theorem 16. 7[ or Theorem 17.31 in the effectively resonant case, that 
there is a surjective map 

(10.1) M + (a) : E^ + (a) -> 8 8eMin+ ( <r )5(M ?l - 1 ), 

ae (mmyo J oo)\(Cv(V)UU, e cv(v)^Ht,*), 

given by taking it S £ , ^ , in+ (cr), microlocally restricting u to a neighbourhood of 
each q giving u q S ££^ c + (<r, q) and sending u to the sum of the leading coefficients 
J2k w k(Y")v k (Y'"), (Y",Y"') e R n_1 , of each of the u q . Since the v k are nor- 
malized eigenfunctions of a harmonic oscillator and the w k are Schwartz functions 
of Y" with seminorms rapidly decreasing in k, the sum is a Schwartz function of 
(Y",Y"'). 

Let us regard © g 5(IR' l_1 ) as a subspace of © g L 2 (R n_1 ), endowed with the norm 

(10.2) \\{w q ) qeMi n+((T )|| 2 = V / \w q {Y)\ 2 dw q ^, dw q , a = 2^0- - V{n(q)) dw q , 

where u> q is the measure induced by Riemannian measure, namely the measure 

x n-(n-rn)/2-^r'J dg 

divided by dx/x and restricted to x = 0. (It takes the form dY" dY" 1 provided 
that the y are normal coordinates, centred at the critical point, for the metric 
h(0,y,dy).) 

The next result is the main content of this section. 

Theorem 10.1. The map M + (a) in (110. ip has a unique extension to an unitary 
isomorphism 

M+(a) : E° css (a) -> e g6M ia + ( < r)£ 2 (R n ~ 1 ). 
Remark 10.2. Here, and throughout this section, we take a E (min Vb, 00) \ Cv(V). 
To prove the theorem, we establish several intermediate results. First we show 
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Proposition 10.3. The space E^ in is dense in E^ s (a) in the topology of 



Proof. The proof is by induction. We consider a sequence r C Ti C • ■ ■ C F„ — 
RP+(cr) as in the previous section, but with the additional condition that the 
radial points are ordered so that, among the points with equal values of v, those 
corresponding to local minima of V are placed last. We shall prove by induction 
that 

(10.3) E~ (<r, Ti n Min + (cr)) is dense in E™ (a, I\) in the topology of E° ss (a). 

For i = 1 there is nothing to prove. Assume that (| 10. 3[1 is true for i — k. Let 
Tk+i \ Tfe = {q}. If q arises from a local minimum of Vq, then using a microlocal 
decomposition, any u £ E^ s (a 7 Tk+i) can be written as the sum of u\ £ E^ s (a, {q}) 
and U2 £ E™ s (a, Ffc). A similar statement is true for u £ E^ s (a, Tk+i H Min + (cr)), 
which proves (|10.3|) for i = k + 1. 

Next suppose that q does not arise from a local minimum of Vq . Then we adapt 
the argument of Proposition 11.6 of [5] to prove (|10.3|) for i = k + 1. We first 
make the assumption that a is not in the point spectrum of H. Using our inductive 
assumption, it is enough to show that E^ s (a, Tk ) is dense in E^ s (a,Tk+i). Let 
u £ E^ s (<r,Tk+i). Let Q £ ^^°(X) be microlocally equal to the identity near 
F/t n Min + (cr), and microsupported sufficiently close to Tk D Min + (cr). Then away 
from Min + (cr), u £ ar 1 /^^ 2 by (ii) of Theorem and thus (H - a)Qu = 
[H, Q]u £ x 1 / 2+e i 2 for some e > 0. This is also true near Min + (cr) since Q is 
microlocally the identity there, so we have (H — o~)Qu £ x x l 2Jr<L L? everywhere. 
This implies that 



since v — u— (Qu — R(a — iO)(H—a)Qu) satisfies (H — cj)v — and v £ x~ 1 / 2+€ L 2 
microlocally for v > 0. 

Now choose a modified potential function V as in Lemma I8.5[ where we choose 
v larger than is(q) but smaller than v{q') for every q 1 £ Tk n Min + (cr). (This is 
possible because of the way we ordered the g,.) Since WF sc (Qm) lies in {v > £>}, 
we have 

(10.5) Qu = R(a + i0)(H - a)Qu. 

Now take u' 3 = ^(x/r^u, where 4> £ C°°(R), <j>(t) = 1 for t > 2, <j>(t) = for 
t < 1 and rj — > as j — > oo. Then u'- £ C°°(X), and Wj defined by 




(10.4) 



u = Qu — R(a — i0)(H — a)Qu, 



wj = R{a + i0){H - a)Qu'j 
converge to Qu in x^ 1 / 2 ^ 6 L? . Our choice of V ensures that 

WF sc (^)nRP + ((r) c F fc . 



Moreover, 



(10.6) 



(H — a)wj converges to (H — a)Qu in x 1 / 2+c i 2 . 



Now define 



R(a-i0)(H -a)wj. 
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Then Uj £ E™ s (<j, r&). We claim that Uj — > u in the topology of E® ss (a). Certainly, 
Uj — > u in x~ 1//2 ~ c L 2 . We must also show that -Bitj — > Bu in x _1 / 2 L 2 , where £> is 
as in (|9.2|) . To do this we write 

Bu.; - Bu = B{wj - i?(<7 - iO)(£T - cr)wj) - B((Id-<9)u + Qu) 
= b(R{(t + iO){H - a)Qu' J - R(a - iO)(H - a)w 3 

+R(a - iO)(H - a)Qu - R(a + iO)(H - a)Qu\ , 

using (|10.4|) and (|10.5|) . and this goes to zero in x~ x l 2 L 2 by (|10.6j) and propagation 
of singularities, Theorem 3.1 of [5], as in the proof of [5] Proposition 11.6]. 
If <t is in the point spectrum of H, then equation (|10.4|) must be replaced by 

(10.7) u = u(Qu-R{a-iO)(H-(r)Qu], 

where II is projection off the L 2 er-eigenspace. Consequently we must define Wj by 
HR(<j + iO)(H — a)Qu'j, and then the rest of the proof goes through. □ 

The second intermediate result we need is 

Proposition 10.4. The Hilbert norm (|9.3|) on the subspace E^ in + (a) C E® ss (a) 
is given by the formula 

(10.8) IMll» s >)= E ^a-V(n(q)) f \M+ (q, <r)u\ 2 dcj q . 

Proof. The proof is the same as the one dimensional case, which is proved in Propo- 
sition 12.6 of [5], so we just give a sketch here. 

Let <j) be as in the proof of Proposition 110.31 Then we can write the natural 
norm (|9.3| on E® ss (a) as a limit 

(10.9) lim i((H - a)Au, <j>{x/r)u) = lim i(Au, [H, (/>(x/r)]u). 

Since u £ x _1 / 2_£ L 2 , the only term in [H, 4>{x/r)] contributing in the limit is 
2{x 2 D x )4>{x/r)(x 2 D x ). The cutoff operator A restricts attention to {v > 0}, and 
the limit vanishes when localized to any region where u £ a; _1 / 2+e L 2 , so we can 
substitute for u a sum of expressions u q as in (|6.24p in the effectively non-resonant 
case, or its analogue in the effectively resonant setting arising from (|7.7| (namely 
e »*(z/)/z ti mes an expression as in (|7.7[) . cf. Remark l6.8p , one for each q £ Min + (cr). 
A straightforward computation then gives (|10.8|) . □ 

Proof of Theorem \10.1\ Proposition ll0.4l shows that M + (a) maps E^ in + (a) into a 
dense subspace of © g L 2 (R™ _1 ), with the Hilbert norm of M+(a)u, u £ E^ in+ (a), 
equal to that of u. By Proposition 110.31 E^ in + (a) is dense in E^ s (a), and by 
Corollary 3.13 of [5], E™ s (a) is dense in E° ss (a). The result follows. □ 

So far we have only considered the microlocal restriction of eigenfunctions near 
radial points q satisfying v(q) > 0. For each critical point of Vq, there are two cor- 
responding radial points with opposite signs of is, and we can equally well consider 
microlocal restriction near radial points with v(q) < 0. This leads to an operator 

M_(<r) : E° css (a) - © ffeM in_(a)i a (R n - 1 ) 
and the analogue of Theorem 110.11 holds also for M_(cr). 
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Definition 10.5. The inverses of M±(a), P±{a) : ©geMini^ 2 ^™ -1 ) -> E° ss (a) 
of M±(a) are called the Poisson operators at energy a. 

We can identify ® q& Min + (a)L 2 (K"" 1 ) and ® ge Min_(o-)£ 2 (K n_1 ) in the obvious 
way, and may therefore assume that the M±(a) have the same range, identified 
with the domain of P± (<r) . 

Corollary 10.6. For a ^ Cv(V), the S-matrix may be identified as the unitary 
operator S(a) = M + (a)P^(a) on ® zeWm L 2 (W 1 - 1 ). 

Remark 10.7. For n — 2, the structure of S(o~) was described rather precisely in [3] 
as an anisotropic Fourier integral operator. 

Theorem 110.11 is essentially a pointwise version of asymptotic completeness in 
a. Integrating gives a version of the usual statement, but some uniformity in a is 
required for this. So we proceed to discuss an extension of part (i) of Theorem 16.71 
that is valid in an interval rather than just at one value. For this purpose, let 
/ C (minVo,oo) be a compact interval disjoint from the set of effectively resonant 
energies, the set of Hessian thresholds and Cv(V). Then for each a € I, the 
sets Min+(<r) C RP+(cr) can be identified; we write Min + (7) for this set. Each 
element of Min+(7) is thus a continuous family q{a) of minimal radial points, with 
q(a) S Mm+(cr). 

Proposition 10.8. Let I C (minVo,oo) be as above, and let the q(a) £ Min + (7) 
be an outgoing radial point associated to a minimum point z ofVo, with Y",Y"' 
the associated coordinates given by (|5.23|) . For any h{o~, •) S C°°(/; 5(K" -1 )) there 
is (f> £ C°°(X) orthogonal to E pp (I) such that for every a £ I , 

F{a)- 1 R{a + *O)0 = ^x^-^w^", o-)v 3 {Y"' , a) + u', 

(10.10) J 

h(o-,Y",Y'") = J2 w i( Y "> a >^ Y '">^ 
o 

where Wj, Vj, Kj and b are as in Provosition 1 6'. 31 F(a) is as in Theorem \3.17\ and 
where u' £ C°°{I;I^\x,M)) for some I > -\. 

Remark 10.9. The statement u' £ C°°(I;Isc(X,M.)) is meant to underline that 
this is a global claim, namely v! £ C°°(I; Is l J(0,M.)) and that it is C°° with values 
in C°°(X) microlocally away from {q(<j);<r £ I}, i.e. for all A £ ^ SC (X) with 
WF' SC (A) n {q(a); a £ 1} = 0, Au' £ C"°(I;d°°(X)). 

Proof. By the construction of Section O for each a £ I there is an approximate 
microlocally outgoing solution u a with f a = (H — a)u a £ C°°(X) and F(a)^ 1 u a 
of the same form as the right hand side of (jlO.lOp . Indeed, the construction is 
smooth in a, in the sense that (d/da) k u £ I s (0, M) for each k and each s < —1/2, 
so that with f(a, .) = / CT (.), we have / £ C°°(I;C°°(X)). Notice that there is no 
need to 'globalize' using Proposition l8.ll since microlocally outgoing solutions over 
sources/sinks (i.e. minima of Vo) are localized at q(o~). 

Let / £ C^°(C x X) be an almost analytic extension of / with compact support, 
so d a f vanishes to infinite order at R x X, and let 

(h = — - [ R(a)d a f da A da. 
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Thus, <j) £ C°°{X) since d a f vanishes to infinite order on the real axis. 

We also claim that (|10.10p holds. Indeed, let a € R, x & x identically 

1 near ao, let x be an almost analytic extension of x of compact support. Thus, 

/(c, •) = f(cro, -)x(o-) + ( CT - o- )5 (<r, •)> /(^ •) = /(°o, -)x(o) + (<7 - (T )5(<7, .) 
with g G C C °°(R X X), <? G C C °°(C x X). Then, writing ct - a = (H - a ) - (H - a), 

4> =^ (J^R(a)d„xdaAda\ f(a 0> .) 

(H — <jq) I R(a)d a g da A da H f d a gda Ada, 

2m J c 2m J c 

where in the last term the identity (H — a)R(a) = Id is used. Since the last term 
vanishes (as g is smooth), and the integral in the second term is in C°°(X), while 
the integral in the first term is x{E), we deduce that 

<\> - fa + (H- a )f ao = (H- a )(u ao + f'J 

for some /' G C°°(IxX). Then if v G E pp (I), (H-a )v = 0, we have v G C°°(X), so 
(4>,v) = {u aa +f ao , (H-ao)v) = 0. Also R(a Q +iO)cf>-R{a +iO)U = f' G C°°(X), 
so R{a$ + i0)4> and R(ao + i0)f ao indeed have the same asymptotics. In particular, 
(|10.10P holds for every a G R. □ 

Now we state asymptotic completeness in a more standard form. 

Theorem 10.10 (Asymptotic completeness). Let I C (minVb,oo) be a compact 
interval as above. Then 

M+(.) o Sp(-) : Ran(n 7 ) E pp (I) -» ® qeM in + (i)L 2 {I x R^ 1 ; 2ndadLu q . a ) 

is unitary. Here, as before, dko q ^ a = 2\J a — V(Tr{q)) duj q . 

Proof. For / G C°°(X) orthogonal to E pp (I), let 

u = u(a) = {2m)- 1 (R(a + i0)f - R{a - i0)f) = Sp(cr)/, 

where Sp(cr) = {2m)- 1 {R{a + iO) - R(a - iO)) 

is the spectral measure. The norm of u in E® ss (a) is given by (i(H — a)Au,u), 
where A is as in (|9.3p . Notice that 

2m{H - a)Au - f = (H - a)A(R(a + iO) - R(a - iO))f - (H - a)R(a + i0)f 
= (H-a)((A-Id)R(a + iO)f- AR(er - iO)/) = (H - a)v, v G C°°(X), 
since 

WF^ C (A) n WF sc (R(a - iO)f) = and WF' SC (A - Id) n WF sc (R(a + iO)f) = 0. 
Hence 

27r|M||c, s(ff) = 2ni{(H-a)Au,u) = (f+(H-a)v,u) = (/,Sp(a)/). 

The right hand side is continuous, hence so is the left hand side. 

Integrating over a in /, denoting the spectral projection of H to I by Hi, and 
using Proposition rXUHl we deduce that M + {a) Sp(a)f is continuous with values in 
L 2 and 

(10.11) IIIWH 2 = 2n J \\M + {a) Sp(a)/|| 2 da, 
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so M+(-) o Sp(-) is an isometry on the orthocomplement of the finite dimensional 
space Epp(I) in the range of 11/. 

It remains to prove that the range is dense in © gS Min£ 2 (-f x R™ -1 ). It suf- 
fices to show that if h £ e ge MinC°°(7 x R"- 1 ), then there is a / £ C°°(X) with 
M + (a) Sp(cr)/ = h(a,.). But this was proved in Proposition 110.81 so the proof of 
the theorem is complete. □ 

Remark 10.11. The results of this section can be related more closely with The- 
orem 19.21 by considering the closure of E^ in + (a) as a subset of E™ s (a) in the 
topology of E^ ss (a) for varying values of s. We have seen in Proposition 1 1 . 31 that 
^Min +( CT ) i s dense, in the topology of E® ss (a). In fact the proof of Proposition ll0.3l 
shows that this is true in the topology of E^ ss (a) for < s < sq, where so is the 
smallest number such that every u £ E™ ic (q), for every q £ RP + (cr) \ Min + (a), 
is in a;~ 1 / 2+s °L 2 locally near n(q); that Sq is strictly positive follows from (ii) of 
Theorem 16. 71 By contrast, E^ in + (a) is closed in the E™ s (a) topology. What hap- 
pens as s increases is that the closure of E^ in +( CT ) m tne ^'oss( cr ) topology changes 
discretely, as s crosses certain values determined by the structure of eigenfunctions 
at the non-minimal critical points. 

One way to understand this is in terms of microlocally incoming eigenfunctions at 
the outgoing radial points, i.e. microlocal eigenfunctions u with scattering wavefront 
set near q is contained in as opposed to $+((7). In Part I we showed (in all 

dimensions) that there are nondegenerate pairings 

£mic,+ (<7, o) x Enac -(q, a) -* C, 

£ c s ss (a)x£- ss s (a)^C 

(Lemma 12.2 and Proposition 12.3 of [5]). The closure of E^ in + (ct), in the topology 
of E^ ss (a), may be identified with the annihilator, in E^ s (a), of the eigenfunctions 
which are in (a) and have scattering wavefront set contained in 

(J $_(g)U{^<0}. 

9 eRP + 0)\Min + (o-) 

This set is trivial for s < sq, and nontrivial for s > sq. The fact that this set 
of eigenfunctions jumps discretely with s in shown in the two dimensional case in 
Section 10 of Part I. 

11. Time-dependent Schrodinger equation 

11.1. Long-time asymptotics. In this final section we apply the earlier results 
to deduce the long-time asymptotics for solutions of the initial value problem 

(11.1) (A + H)u = 0, u\ t=0 = «o, u £ C°°{X), 

for a dense set (in L 2 E pp (H)) of initial data. 

Our approach is to use the spectral resolution of uq and the functional calculus. 
In this way, we deduce the long-time asymptotics of u from the asymptotics of 
generalized eigenfunctions of H using the stationary phase lemma. 

We first define the space X^ch on which the asymptotics of the solution u of 
(|11.1[) will be described. Let us first choose a global boundary defining function 
x satisfying (| 1 . 1 p ; we can specify, for example, that x = 1 outside a collar neigh- 
bourhood of dX. We then introduce the variable r = tx, where t is time. Let 
us compactify the real r-line R to an interval R using t _1 as a boundary defining 
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function near r = oo, and — r 1 as a boundary defining function near r = — oo. 
Then we define 

(11.2) X Sch = X x I T 

Thus Xg c h is a compact manifold with corners, with boundary hypersurfaces if (the 
'infinity face') at t = ±oo (or t = ±oo), naturally diffcomorphic to two copies of X 
(one at t = +oo, one at t = — oo), and a boundary hypersurface af (the 'asymptotic 
face') diffcomorphic to dX x R T . At af, every point with t > corresponds to 
t = +oo and every point with t < corresponds to t = — oo, so this is the place to 
look for long-time (and large-distance) asymptotics of the Schrodinger wave u. The 
variable r has an interpretation of inverse speed; a particle travelling asymptotically 
radially at speed t^ 1 will end up at af after infinite time at t = tq. 

We now specify a good subset of L 2 initial data uq, for which the asymptotics as 
t — > +oo of the solution, u, to (| 1 1 . 1[) are particularly simple. Let L c (minVo,oo) 
be a compact interval disjoint from Cv(V) and from the set of effectively resonant 
energies and Hessian thresholds. Let (h(cr,-))q £ C°°(I; <S(K™ -1 )) be a collection 
of smooth functions from / into Schwartz functions of n — 1 variables, one for 
each q e Min+(J), and let <j) = (j>(I,h) = Y, q <i>(I > K) e C°°{X) be the function 
constructed in Proposition ll0.8l Let 

Ai = {<l>(I,hy,h(a,-) GC°°(/;5(R"- 1 ))} and .4 = ^.4/ 

i 

be the (algebraic) vector space sum of Ai over all such I as above. It is clear from 
Theorem 110. 101 that Aj is dense in Ran 11/ (if) Q E pp (I), and hence that A is dense 
in L 2 Q E pp (H) — H ac (H). To give the asymptotics of (|11.1|) with initial data from 
A it suffices to give the asymptotics starting from uq = <f)(I, h) for some h as above. 

Theorem 11.1. Suppose that I is as above and that <fi — <j)(I,h) £ Aj. Let u(-,t) 
be the solution of (|11.1|) with initial data uq = (f>, regarded as a function on X^ch- 
Then u has trivial asymptotics (i.e. u and all its derivatives are 0(t~°°)) at if. 
Also, if w E dX is not a local minimum of Vq, and r > 0, then u has trivial 
asymptotics in a neighbourhood of {w, r) € af. 

Let z be a local minimum ofVo, and let (Y" ,Y'") be the coordinates given by 
(|5.23| , where a is determined in terms of r by (|11.4|) . Then, in a neighbourhood of 
(z,t) S af u takes the form 
(11.3) 

u(x,Y",Y'",t) = CT-^ 2 Y,^ lb ^ +1/2 e^^^/ x w 1 (Y'' l( 7(T))v J (Y''',a(T))+u' 1 

j 

where h(a(r),Y" ,Y"') = £ Wj (Y", a{r)) Vj (F'", a(r)), c=^j=e-^'\ 

j 

(11.4) a(T ) = y o(z ) + _L 

b is as in (|6.6p and (|6 .4[) . Kj is as in (|6 . T[) . $ is a smooth function of y and t, h 
is decomposed as in Proposition \10.8l and vl decays faster than the leading term. 

Proof. Let v(a) = Sp(a)(/) = {2ni)- 1 (R(a + iO) - R(a - i0))<f>. Then 

u(t, ■) = — [ e - ita (R(a + iO) - R(a - iO)Uda. 
2tti Jj 
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Shifting the contour of integration shows that, as t — > oo, R(a — iO)cj> has trivial 
asymptotics. Hence it is enough to consider 

(11.5) u(t,-) = — [ e- lta R(a + iOUda. 

2m J j 

Let F(a) be the FIO constructed in Theorem I3.17[ which conjugates x~ 1 (H — 
a) to normal form microlocally near the point q £ RP+(ct) with n(q) — z. By 
construction, F(a)~ 1 R(a + iO)cp has asymptotics (jlO.lQI) for every a. Since F(a) is 
a smooth family of FIOs, and (|10.10p is a Legendre distribution associated to the 
zero section (i.e. it is conormal at x — with no oscillatory factor), it follows that 
R(a + iO)(j> itself has asymptotics 

(11.6) R{<7 + lO)(t> = x -ib-iK je my,a)/x a( y", Y'" , X, a) + v' , 

where ■!>(•, c) is a smooth function, parametrizing the image of the zero section 
under the canonical relation of F(a) (as in (|6.24[0 . By assumption, a is smooth in 
a , conormal in x and Schwartz in (Y", Y'"). At the critical point z we have 

(11.7) $(z,(t) = y/cr-Vo(z), Z = n{q), q £ Min + (<j). 
We may substitute (|11.6p into (|11.5p and compute 

(11.8) u(t, -) = ^lj z~ lt ° (e l ^y^l x a{Y", Y'", x, a) + w') da, 

exploiting the smoothness of $ and a in a. 

Let p £ X be an interior point. Then (R(a ± i0)4>)(p) is a smooth function of a 
by Proposition ll0.8l It follows that for a fixed interior point p the integral (|11.8p is 
rapidly decreasing as t — > oo, being the Fourier transform of a smooth, compactly 
supported function. Hence the asymptotics of u are trivial at if. 

To investigate asymptotics at af, where x — * 0, we rewrite (|11.8p as 
(11.9) 

u(r, x, Y", Y'") = ^- J e K-™+*{y,*))/x ^ Y ", Y'" , x, a) + v'(Y", Y"',x, a)) da, 

and apply stationary phase to the integral. We first note that for any w E dX 
which is not a local minimum of Vo, the integrand is rapidly decreasing as x — > 
in a neighbourhood of (w, a) £ OX x /, uniformly in a, so u is rapidly decreasing 
as x — > in a neighbourhood of (w, r) £ af. So we may restrict attention to a 
neighbourhood of [z, r) £ af, where z is a local minimum of Vo- 

To do this we apply stationary phase to (|11.9p . The phase is critical when 
t = d<j <&(?/, a). Since $ is smooth in y, this gives 

^critical = d a $(z, a) + 0{YiX Ti ) 

and, since a is Schwartz in Y, to compute the expansion of u to leading order we 
may drop the 0(YiX Ti ) terms when we substitute t = r cr itical into a in (|11.8p . Since 
$(z,a) is given by (|11.7p . we may therefore take r in the argument of a to be 
given by r = d a &(z,a) which implies (|11.4p . Moreover, the Hessian of the phase 
function at the critical point is {Ax)^ 1 {a — Vo(z))~ 3//2 = 2x~ 1 t 3 . The stationary 
phase lemma then gives (|11.3p . with $(y, r) = — t<t(t) + $(y,a(r)). □ 

Remark 11.2. Equation (|11.4p is just the energy equation 'total energy = potential 
energy+ kinetic energy' at infinity, since 1/r is the asymptotic speed. The factor 
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1/4 comes from the fact that in writing our Hamiltonian as A + V , we have taken 
the value of mass to be 1/2 in our units. 

Remark 11.3. We may not replace ^(y, r) with ^(z, t) in (|11.3p . due to the singular 
factor l/x in the phase. In fact, if we expand ^(y, t) in a Taylor series about y = z, 
written in terms of the variables Yi — yi jx Ti , then we get an asymptotic expansion 
involving polynomials y a in the variables yi multiplied by nonncgative powers x r , 
where r — ^2 a^i. We may discard all terms in the Taylor series of \P with r > 1 
since these will only contribute to the term u' decaying faster than the leading 
term, but we must keep all terms with r < 1. The number of such terms is always 
finite, but depends on a — Vq(z) and the eigenvalues of the Hessian of Vo at z. 

11.2. Asymptotic completeness: time dependent formulation. We see that 
solutions of the time dependent Schrodinger equation (at least those with initial 
data in A) have expansions at af which are equivalent to first spectrally resolving the 
initial data and looking at the expansion of the corresponding family of generalized 
eigenfunctions; the variable a in the time-independent setting, and r in the time- 
dependent setting, play equivalent roles and are linked by (|11.4|) . In view of this, 
we can recast Theorem 1 1 . 1 01 in time-dependent terms as follows: 

Theorem 11.4. Let I and Ai be as in Theorem \11.1[ let uq £ Ai and let u be 

the solution of the time dependent Schrodinger equation (jTTTTj) with initial data uq . 
For a given local minimum zofVo, let Min+(J) be the associated family of outgoing 
radial points, and let Q = ^ Qj , b and Kj be as in Proposition \6.3[ The map 



whose existence is guaranteed by Theorem \11.1\ extends uniquely by linearity and 
continuity to a unitary isomorphism 



Here io qiT is the measure in (|10.2|) and r = r(q, a) is given by (lll.4[) . 

Remark 11.5. The operator e llo s x Q simply removes the factors of x~~ lak in the 
expansion pi.3[) . so that we can take a limit as x — > 0. 

Remark 11.6. The measure in (|11.11[> should be thought of as the product of 
T~ 1 ujq tT , which is the measure in Proposition 110.41 tensored with the measure 
da = T- 3 dr/2. 

Proof. Let uq — <f)(I, h) be as in Theorem lll.il We may take the L 2 norm of (|11.3|) 
for a fixed t, and take the limit as t — > oo. To do this, we write x = r/t and 
integrate with respect to the measure on X which is given by tr~ 2 x~ 1 ~ 2lnlb dYdT. 
If we just look at the principal term in (|11.3[) then the powers of t cancel exactly 
and we get 



(11.10) 




(11.11) 




(11.12) 




Since to, 



t x dY, and da = r 3 dr/2 using (|11.4|) . this is given by 
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The expression Y]^ WjVj is equal to M+(q, a)((R(o~ + iO) — R(a — iO))uo, or equiva- 
lently 2mM + (q, a) Sp(a)u Q . Also, the norm on e 9 L 2 (M ,l_1 ) is given by (|TIP)) . So 
we get, using Theorems 110.11 and 110. 101 

(mi =2nJ2 I \\M + {q,o-)S V (o-)u Q \\ 2 da = \\u Q \\ 2 L2 . 

But (|11.12p is precisely the square of the norm of the right hand side of (] 1 1 . 10|) . 
So we have established the conclusion of the theorem for the principal term in the 
asymptotic expansion in (| 11 . Since the remainder term u' decays faster than 
the principal term, the L 2 norm of u'{-,t) goes to zero as t — * oo, so the proof is 
complete. □ 

From Theorem 111.41 we can deduce the following result first proved by Herbst 
and Skibsted, using a direct method involving the uncertainty principle, rather than 
proceeding via the structure of generalized eigenfunctions as here. 

Corollary 11.7 (Absence of L 2 channels at non-minimal critical points). Let x € 
C°°(X) vanish in a neighbourhood of the local minima of on dX. Let u be the 
solution of (jll.ip on X x K with initial value uq € L 2 (X) Q E pp (H). Then 

(11.13) lim \\ X u(t,-)\\ L 3 ix) ^0. 

t— >oo v ' 

Proof. We may assume that < x < 1 without loss of generality. 

Let e > be given. Then by density of A in L 2 QE pp (H), we can find 4> £ A, with 
4> equal to a sum of a finite number of 4>j{Ij, hj) € Ai , such that ||uo — 4>\\l 2 < e - 
Without loss of generality we may assume that all the Ij are disjoint. Let u 1 be the 
solution with initial condition (j). By direct calculation from (|11.3|) we find that 

t Um ||(1 -X)u'(*, Ollk=£ / V^-V a (n(q))\\h,\\ 2 L HR^)d'y, 

which by Theorem 110. 101 is equal to ||0||? 2 . But by unitarity of e~ ltH , we have 

= U\\ 2 L2 for eachi. 
Since < x < L an elementary calculation shows that 

ll(i-xy||| 2 + ||x^||i 3 <lK||| 2 , 

which implies that 

lim || xu '(i,.)|| 2 L2 =0. 

So (|11.13|) is true for u 1 . On the other hand, 

limsupj^ \\x(u(t, •) - u'(t, -))\\ L 2 < e, 

so limsup^^ ||x u (^i ')IU 2 — e - Since this is true for every e > 0, the result follows. 

□ 
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11.3. Comparison with results of Herbst-Skibsted. We first show that our 
results on the asymptotics of the solutions to the time-dependent equation (jTTTTJ) 
are consistent with the comparison dynamics of Herbst-Skibsted |11) . Herbst and 
Skibsted define comparison dynamics, i.e. a family of unitary operators Uo(t) for 
a given local minimum of Vo and for either a 'low energy' range or a 'high energy' 
range which depends on the behaviour of the r, (a) from Lemma 12.71 It has the 
property that the strong limit 

lim e ltH U Q (t) 

t — >oo 

exists in L 2 (X) and defines a unitary wave operator. 

Let us compare their results on long-time asymptotics with ours. For simplicity, 
we consider the 'very low energy' energy interval in which all of the exponents 
Ti are complex, with real part 1/2 (this is 'below the Hessian threshold', in our 
terminology). For simplicity we also assume, as do Herbst and Skibsted, that 
Vq(z) = 0. In this case, the exponent — ib in (| 1 1 .3|) is equal to (n — l)/4, and there 
are no Y" variables. Moreover, the function $(y, ct) is equal to — M 2 /4) (see 

[5], section 7, particularly (7.22) and (7.23) for the case n — 2), which implies that 
vff(y 5 r) = (1- |y| 2 /4)/r. If we substitute x = r/t into (TO} then we get 

(11.14) c^t-M^-V^^M/w-^eaii-lylV^^^^'^), 

3 

To compare this with Herbst and Skibsted's comparison dynamics, we adopt their 
notation: we decompose the variable x £ M. n as x — (xi, x^) where (1,0, ... , 0) is 
the point on the sphere at infinity where Vo has a local minimum, and x 1 - are n—1 
orthogonal linear coordinates. We can identify our boundary defining function x 
with 1/xi. Thus r = tjx\ and y — x /x%, and we can write the expression above 
as 

(11.15) c^t-^-^^-^^T^-^/^^e^e^^l'/^^WwjfarL/VST.r). 

3 

In this very low-energy case the Herbst-Skibsted comparison dynamics is given 
explicitly by 

(11.16) U (t) = 5 f - 1 /2e 4 l^l 2 / 4 e-^/ 2 e- l(lost)H2 f7 , 
where the operator pi stands for D Xi — —id Xi , S t -i/2 is the scaling 

S t - 1/2 f(x 1 ,x x )=t-^- i y 4 f(x 1 ,t-^ 2 x x ), 
the operator H 2 is given by 

ff 2 = ^ + i<^,(p r v 2 )-id/4)^) 

(where V^ 2 ' is the Hessian of Vo at the critical point), and finally Uq is an arbitrary 
unitary operator. 

To compare this to our long-time asymptotic expansion (|11.15| . it is convenient 
to take Uq to be inverse Fourier transform mapping functions of p\ to functions of 
x\. Then H2 is a family of harmonic oscillators parametrized by p\. The operator 
e -itp 1 /2 ac tj n g on W(pi : x ± ) then takes the form 

(27T)- 1 / 2 / e^e-^^WipuX^dpt 
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and by stationary phase we see that the large t asymptotics of this operation is 
given by 

W( Pl ,x X ) Hr^e^^'^f^i 1 ). 

Let us expand W(pi, x x ) in eigenfunctions of the operator H2 — H^ipi) as 

W(p 1 ,x ± ) =^2u j (p l )xj{x ± ,pi), 
j 

and write r for t/x\. A computation shows that S t -i/2H2S t i/2 is equal to p{ X Q 
where Q is the operator from (|6.8p . The comparison dynamics therefore maps W 
to 

(11.17) r( n - 1 )/^ 1 / 2 e it / aTS e- i l <BX l J, / 4 *^;t««*w i (r- 1 )xi(ar L /Vt,T- 1 ). 

j 

This agrees with (|11.15|) . if we identify Wj(r) with T _1 o;j(r _1 ), and Vj(Y'",T) with 
Xj(y"T _1 / 2 , t^ 1 ). (The imaginary powers of r simply amount to a different choice 
of normalized eigenfunction Vj. Also there are some discrepancies of factors of 
1/2 since Hcrbst-Skibsted's operator is half the Laplacian plus V.) Thus, the two 
expansions are consistent. 

In the high energy regime, it is easier to check the agreement of the two expan- 
sions. In this case, there are no Y'" variables. The Herbst-Skibsted comparison 
dynamics takes the form 

where S(t, x) is a solution to the eikonal equation 

(11.18) dtS(t,x) + ^\VxS(t,x)\ 2 + V(x) = 
and k(t,x) is the energy function 

y = l\VxS(t,x)f + V(x). 

To make the link with our long time expansion (|11.3[) , we begin by showing that 
S corresponds to our phase function ty/x. Indeed, "J/ is obtained from the phase 
function in (|11.6[) by performing stationary phase as in (|11.9[) . The phase function 
<5>/x parametrizes a Legendrian submanifold which is the image of the zero section 
under the FIO F in Remark 16.81 Since the zero section is the flowout from the 
critical point in the eigendirections (of the linearized flow) with eigenvalues Ar.; (as 
opposed to A(l — rj)), as can be computed easily from (|2.13[) . the same is true of 
the Legendrian submanifold parametrized by in particular, it corresponds pre- 
cisely to Herbst-Skibsted's Lagrangian submanifold M.^ parametrized by S (using 
the correspondence between conic Lagrangian submanifolds and Legendre mani- 
folds 'at infinity'); see [11], Theorem 2.1. Then the way /x is obtained from $/x 
is exactly the same as the Legendre transform by which Herbst-Skibsted obtain S 
from S (see [IT], p559), with k 2 corresponding to our a and S corresponds to our 
^. Moreover, from [TTJ, p561, we have 

w . = t -n j {k) k x-0 j {k)^ 1 + 2^(k)) Uj + 0( Uj \u\), 

In our notation, [3j(k) = —Vj, Uj = yj and t = t/x. Setting Y" = yjX~ Tj as above, 
we get 

Wj = 9j (k)Y!' + 0(x^). 
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Thus, up to an energy-dependent factor gj(k), the coordinate w in Herbst-Skibsted 
is equivalent to our Y" . The asymptotics ()1 1 .3[) in this regime (where now there 
are no Y'" variables, hence no sum over j) thus take the form 

which is consistent with the Herbst-Skibsted comparison dynamics. 

It is a little more difficult to make the link between our asymptotics and Hcrbst- 
Skibsted's in the low-energy regime where not all the r t have real part equal to 1/2, 
but the real parts are all at least 1/3. We can, however, offer some explanation 
as to why the low-energy comparison dynamics fails to work above this energy 
level. Referring to Remark 1 1 1 .31 above this energy we cannot approximate the 
function ^/(y, cr(r)) by its quadratic approximation; we need to include at least 
cubic terms in the Taylor series of ^ at y = z. These in turn depend on the cubic 
terms in the Taylor series for Vq at z. The Herbst-Skibsted low energy comparison 
dynamics neglects these terms. It cannot therefore be expected to provide an 
accurate approximation to the long-time asymptotics of solutions to (|11.1[) , since we 
have seen that in (|11.3p that one cannot replace 'J by its quadratic approximation. 

We emphasize that our long time asymptotic formula (|11.3p works for all energies 
(except for the discrete set of eigenvalues, effectively resonant energies and Hessian 
thresholds), whereas in Herbst and Skibsted's results there is a gap of 'intermediate 
energies' in which they do not give any comparison dynamics. The formula (|11.3[) 
correctly interpolates between low energies, below the Hessian threshold, and high 
energies, where all the exponents r$ are real. 



Appendix A. Errata for [5] 

A.l. Correction to the proof of (5[ Proposition 6.7]. There is an error in the 
proof of Proposition 6.7]: with the stated assumptions, the proof of Proposition 
6.7 needs to be two-step, and the conclusion is slightly modified, although this does 
not affect any of its applications, in particular (5j Proposition 6.9], which is its only 
use in [5]. Below equation numbers of the form (6.xx) refer to (5J> while equation 
numbers of the form (A.xx) refer to this appendix. 

The error in the proof arises from the microlocalizers Q £ 4 , ~ oo,0 (A) considered 
there, in (6.27), so we recall the assumptions on it. With O m a neighborhood of q 
as (6.24) or (6.25), we assume that WF^ C (Q) C O m q & WF^Qd-Q), 

i[Q*Q, P-a]= x 1/2 (B*B + G)x l/2 + x 1/2 Fx 1/2 , where 

' A '" B, Fef^(O), Ge^(i), ^WF;(F), 

and in addition, F satisfies WFg C (F) C {v < v(q)}. (This condition on F ensures 
that WF' SC (F) n WF sc (u) = for the application in Section 9].) 

In fact, due to the two step nature of the proof below, we also need another mi- 
crolocalizer Q' £ \E , ~°°'°(A) satisfying analogous assumptions with B, etc., replaced 
by B', etc., 

(A.2) i[(Q')*Q', P-<t] = x^UB'YB 1 + G^x 1 ' 2 + x 1 ' 2 F l x 1 ' 2 , 

with properties analogous to (|A.1|I . except that WFg C (Q') C 0' m , etc., where 0' m 
is the elliptic set of Q. 
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The following is a slightly modified version of Proposition IA.1| in that we need 
to assume the existence of Q' as above, and that the conclusion is on the elliptic 
set of Q' rather than that of Q. 

Proposition A.l. (modified version of 15, Proposition 6.7],) Suppose that m > 0, 
s < —1/2, q £ RP+(<t), a $ Cv(V), either (6.14) or (6.15) hold, and let O m be as 
in (6.24) (or (6.25)). Suppose thatu£ 7 s ( c ),m_1 (O m , M), WF sc ((P-o»nO m = 
and that there exists Q, Q' £ H!~ c °° ,0 (O m ) elliptic at q that satisfies (|A.1[) - (|A.2[) with 
WF' sc (F)nWF sc (a) = 0, WF' SC (F') n WF sc (u) = 0. Then u £ I$' m (0", M) where 
O" is the elliptic set of Q' . 

The issue with the argument presented in [5J Proof of Proposition 6.7]: is that it 
gains a whole extra factor in the module at once: u £ 7 S c 1 {O m , Ai), is assumed, 
and u £ Isc m (0' ,A4) is concluded. Now, the novel part of such a statement, 
corresponding to the terms arising from factors from the module M. C V E'~ 00 '~ 1 (X), 
is properly dealt with in the (erroneous) proof presented in [5|. However, there is a 
problem with the microlocalizer Q unless (6.27) is strengthened to make the error 
term G have two orders higher decay than the main term, i.e. to make it order 
(0,2). This is of course the same issue as what makes one gain 1/2 order at a 
time usually in positive commutator proofs for the propagation of singularities for 
operators of real principal type. Factors from the module M. are fine because they 
essentially get reproduced by the commutator with P — a. The problem is that 
G cannot be written as a multiple of Q, in general. Technically, this shows up in 
(6.29) where e|| A ajS M'|| 2 cannot be absorbed in the left hand side for it does not 
have a factor of Q. (One needs to remember that Au' is the vector of QA a ^ s u', 
so all terms arising by commutators with the module generators are OK, the only 
issue is the microlocalizer Q.) 

This error is easily remedied by a two-step argument. The cost of this is that 
the open set on which we conclude regularity is shrunk slightly from the elliptic 
set of Q to that of Q', although in relevant situations one can usually recover the 
original statement of [5J Proposition 6.7] easily as in [51 Proposition 6.9]. First, the 
argument given in [5J Proof of Proposition 6.7] proves the following lemma. 

Lemma A. 2. Suppose that m > 0, r < —1/2, q £ RP + (er), a £ Cv(V), ei- 
ther (6.14) or (6.15) hold, and let O m be as in (6.24) ( or (6.25)). Suppose that 
u £ 7 s ( c r),m_1 (O m ,7W), WF SC ((P - a)u) n O m = and that there exists Q £ 
*- m '°(O m ) elliptic at q that satisfies (TaTT]) with WF' SC (P) n WF sc (u) = 0. Then 
u £ /ic~ 1 ^' m (0', M) where O' is the elliptic set of Q. 

Notice that under the same hypothesis as Proposition IA.ll this lemma proves 
regularity under A4 m (as Proposition I A . 1 1) . but does so at the cost of losing half an 
order of decay: u £ l£~ 1/2) ' m (0' ,M) rather than u £ I^' m (0',M). 

Proof of Lemma \ A. 2[ With the notation of [5J Proof of Proposition 6.7], let s = 
r — 1/2 (so in particular s < —1/2), let A a>s , etc., be as there. Then the pairing 
(A atS u' , GA ajS u') (where v! will be regularizations of u) is controlled by the a priori 
control of v! in 7 s ( c +1/2) ' m_1 (O m , M) = I^ m ~ l {O m ,M)- Indeed, x 1 ' 2 A a , s and 
x~ l l 2 GA ayS are both the product of an element of ^sc' s+1 ^ 2 - ) (O m ) and m factors in 
the module M C ^^(Om), hence in particular can be thought of (by combining 
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the factor from ^ic' s+1 ^ 2 ' (O rn ) with a factor from M.) as the product of an element 
of l4c'" s " 1/2) (O m ) with m - 1 factors in M. So this gives u € I^ } ' m (0',M) = 
proving the lemma. □ 

We can now prove Proposition lA.il 

Proof of Provosition \A.l[ Lemma lA.21 shows that u G /ic _1 ^ 2 ^ m (0', M.) with O' 
as in Lemma [A. 21 With this additional knowledge, the argument stated in [5l Proof 
of Proposition 6.7], applied with the same s, goes through. (But now we apply it 
with Q replaced by Q', etc!) Indeed, the pairing (A a ^ s u' , G 1 A a _ s v!) is controlled 
by the a priori information, as x 1 / 2 A a ^ s u' = A a s _ 1 / 2 U \ so it is controled in L 2 if 
u 1 is a priori controlled in I$c 1/>2 ' > ' m (O' , A4) (which we just have proved), and a 
similar conclusion holds for x~ x / 2 G' A a s u' as x~ 

-vac" e *°c 1/2 PO just like x 1 ! 2 is. 
Thus, u £ I^ hrn (0",M), with O" the elliptic set of Q', as desired. This finishes 
the proof. □ 

A. 2. Correction to Proposition 9.4 of |5j. The proof of Proposition 9.4 contains 
the statement " Since r\ < 0, the vector field xd x + r\yd y is nonresonant" , which 
is false. To correct the proof, that statement should be deleted and the sentence 
following it replaced by: "By a change of coordinates x' = a(y)x, y' = b(y)y, where 
a, b e C°° near y — satisfy the ODEs 

a'(y) = ny) = --^kry a(0)= m = l 

ri+yG(y) ri+yG(y) 

the F and G terms are eliminated and the vector field becomes 

2V -Mx') 2 D xl )+r iy {x'D y ,)), 



a(y) 

modulo terms in x 2 Ai 2 and subprincipal terms." This proves the proposition apart 
from the prefactor of a(y) _1 in front of Pq which is irrelevant for the application of 
this proposition. 

Of course, Proposition 9.4 also follows by applying the results of the present 
paper, noting that the case considered there is effectively nonresonant. 
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